心墙土石坝水利枢纽工程毕业设计计算书.doc

《心墙土石坝水利枢纽工程毕业设计计算书.doc》由会员分享，可在线阅读，更多相关《心墙土石坝水利枢纽工程毕业设计计算书.doc（41页珍藏版）》请在沃文网上搜索。

1、目 录第一节 调洪计算2第二节 大坝高程的确定7第三节 土石料的设计13第四节 渗流计算18第五节 稳定分析25第六节 细部结构计算32第七节 泄水建筑物的计算34第八节 施工组织设计39第一节 调洪计算主要建筑物为2级，次要建筑物为3级，临时建筑物为4级。永久建筑物洪水标准：正常运用（设计）洪水重现期100年；非常运用（校核）洪水重现期2000年。采用隧洞泄洪方案水库运用方式：洪水来临时用闸门控制下泄流量等于来流量，水库保持汛前限制水位不变，当来流量继续加大，则闸门全开，下泄流量随水位的升高而加大，流态为自由泄流。调洪演算原理：采用以峰控制的同倍比放大法对典型洪水进行放大，得出设计与校核洪水

2、过程线。拟定几组不同堰顶高程及孔口宽度B的方案。堰顶自由泄流公式Q=Bm(2g)1/2H3/2可确定设计洪水和校核洪水情况下的起调流量Q起，由Q起开始，假定三条泄洪过程线（为简便计算，假设都为直线），在洪水过程线上查出Q泄，并求出相应的蓄水库容V。根据库容水位关系曲线可得相应的库水位H，由三组（Q泄，H）绘制的QH曲线与由Q=Bm(2g)1/2H3/2绘制的QH曲线相交，所得交点即为所要求的下泄流量及相应水位。 方案一：I=2818m, B=7m起调流量=0.90.5710.8=495.22 m3/s设计洪水时：Q设 =1680 m3/s计算QH曲线列表如下：表1-1 设计水位流量关系表（方案

3、一）假设泄水流量 Q泄 （/s）水库拦蓄洪水 V（10）水库总水量V （10）水库水位H（m）56033.8579.82830.170028.9574.92829.998020.8566.82829.6QH曲线与Q=Bm(2g)1/2H3/2 的交点为：Q泄 =585 m3/s，H=2830.06 m校核洪水时： Q校=2320 m3/s计算QH曲线列表如下：表1-2 校核水位流量关系表（方案一）假设泄水流量 Q泄 （/s）水库拦蓄洪水 V（10）水库总水量V（10）水库水位H（m）582.461.8607.82831.3776.652.4598.42830.9970.744.1590.128

4、30.6QH曲线与Q=Bm(2g)1/2H3/2 的交点为：Q泄 =660 m3/s，H=2831.125 m图1-1 方案一调洪演算结果方案二：I=2819m, B=7m起调流量=0.90.579.8=428.05 m3/s设计洪水时：Q设 =1680 m3/s计算QH曲线列表如下：表1-3 设计水位流量关系表（方案二）假设泄水流量 Q泄 （/s）水库拦蓄洪水 V（10）水库总水量V （10）水库水位H（m）49038.7584.72830.456035.8581.82830.270030.4576.42830.098021.7567.72829.6QH曲线与Q=Bm(2g)1/2H3/2

5、的交点为：Q泄 =530 m3/s，H=2830.275m校核洪水时： Q校=2320 m3/s计算QH曲线列表如下：表1-4 校核水位流量关系表（方案二）假设泄水流量Q泄（/s）水库拦蓄洪水 V（10）水库总水量V（10）水库水位H（m）485.470.2616.22831.8776.654.5600.52831.0970.745.5591.52830.7QH曲线与Q=Bm(2g)1/2H3/2 的交点为：Q泄 =610 m3/s，H=2831.425m图1-2 方案二调洪演算结果方案三：I=2819m, B=8m起调流量=0.90.589.8=489.2 m3/s设计洪水时：Q设 =168

6、0 m3/s计算QH曲线列表如下：表1-5 设计水位流量关系表（方案三）假设泄水流量Q泄 （/s）水库拦蓄洪水V（10）水库总水量V（10）水库水位H（m）49036.5582.52830.356034.0580.02830.270029.1575.12829.998020.9566.92829.7QH曲线与Q=Bm(2g)1/2H3/2 的交点为：Q泄 =590 m3/s，H=2830.125m校核洪水时： Q校=2320 m3/s计算QH曲线列表如下：表1-6 校核水位流量关系表（方案三）假设泄水流量Q泄（/s）水库拦蓄洪水V（10）水库总水量V（10）水库水位H（m）582.462.06

7、08.02831.3776.652.6598.62831.0970.744.3590.32830.6QH曲线与Q=Bm(2g)1/2H3/2 的交点为：Q泄 =675 m3/s，H=2831.15m图1-3 方案三调洪演算结果方案四：I=2820m, B=8m起调流量=0.90.588.8=416.27 m3/s设计洪水时：Q设 =1680 m3/s计算QH曲线列表如下：表1-7 设计水位流量关系表（方案四）假设泄水流量Q泄（/s）水库拦蓄洪水V（10）水库总水量V（10）水库水位H（m）49039.2585.22830.456036.2582.22830.370030.7576.72830.

8、098021.8567.82829.6QH曲线与Q=Bm(2g)1/2H3/2 的交点为：Q泄 =535 m3/s，H=2830.34m校核洪水时： Q校=2320 m3/s计算QH曲线列表如下：表1-8 校核水位流量关系表（方案四）假设泄水流量Q泄（/s）水库拦蓄洪水V（10）水库总水量V（10）水库水位H（m）582.465.2611.22831.5776.654.8600.82831.0970.745.7591.72830.7QH曲线与Q=Bm(2g)1/2H3/2 的交点为：Q泄 =620 m3/s，H=2831.4m图1-4 方案四调洪演算结果第二节 大坝高程的确定坝顶在水库静水位以

9、上的超高按下式确定： 其中：y-坝顶超高；R-最大波浪在坝顶的爬高；e-最大风壅水面高度；A-安全超高。该坝为二级建筑物，设计时取A=1.0，校核时取A=0.5坝顶高程等于水库静水位与坝顶超高之和，应按以下运用条件计算，取其最大值:1. 设计水位加正常运用条件下的坝顶超高;2. 正常蓄水位加正常运用条件下的坝顶超高;3. 校核洪水位加非常运用条件下的坝顶超高； 4. 正常蓄水位加非常运用条件下的坝顶超高。波浪的平均波高和平均波周期宜采用莆田试验站公式g hm/w2 = T=4.438hmL= hm平均波高，mT平均波周期，sD风区长度km水域平均水深，mW 计算风速, m/s风壅高度可按下式计

10、算: 式中 e计算处的风壅水面高度m D风区长度km K综合摩阻系数3.610计算风向与坝轴线的夹角22.5（1）设计水位加正常运用条件下的坝顶超高设计水位 2830.275m 吹程D=15 风速W=25m/s 坝前水深Hm=80.275m =22.5解得hm=1.31423m T=4.438hm=5.0877sL= =40.41439m波浪爬高设计波浪爬高值应根据工程等级确定，2级坝采用累积频率为1%的爬高值.正向来波在单坡上的平均波浪爬高可按下式或有关规定计算:m=2.5= 0.771.0 / =2.084m式中 Rm平均波浪爬高 m单坡的坡度系数,若坡角为a,即等于cota K斜坡的糙率

11、渗透性系数，根据护面类型查规范得0.77 Kw经验系数,查规范得1.0查规范1累积频率下的爬高与平均爬高的比值为2.23 R/R=2.23R1%=2.232.084=4.648m因风向与坝轴线的法线成22.5,波浪爬高应按等于正向来波计算爬高值乘以折减系数K查规范得0.95 KR1%=0.95 4.648=4.42m坝顶在水库静水位以上的超高按下式确定：y=R+e+A=4.42+0.0198+1.00=5.4398m坝顶高程2830.275+5.4398=2835.7m坝高2835.7-2750=85.7m为预防坝体竣工后的沉陷,预留0.4%的坝高坝高+0.4%沉陷=2835.7+85.70.

12、004=2836m（2）校核洪水位加非常运用条件下的坝顶超高校核洪水位 2831.425m 吹程D=15 风速W=14m/s 坝前水深Hm=81.425m =22.5解得hm=0.68904m T=4.438hm=3.683905sL= =21.18877m波浪爬高设计波浪爬高值应根据工程等级确定，2级坝采用累积频率为1%的爬高值.正向来波在单坡上的平均波浪爬高可按下式或有关规定计算:m=2.5= 0.771.0 / =1.093m查规范不同累积频率下的爬高与平均爬高的比值为2.23 R/R=2.23R1%=2.231.093=2.44m因风向与坝轴线的法线成22.5,波浪爬高应按等于正向来波

13、计算爬高值乘以折减系数K查规范得0.95KR1%=0.95 2.44=2.315m坝顶在水库静水位以上的超高按下式确定：y=R+e+A=2.315+0.00612+0.5=2.82m坝顶高程2831.425+2.82=2834.25m坝高2834.25-2750=84.25m为预防坝体竣工后的沉陷,预留0.4%的坝高坝高+0.4%沉陷=2834.25+84.250.004=2834.6m（3）正常蓄水位加正常运用条件下的坝顶超高正常蓄水位 2828.8m 吹程D=15 风速W=25m/s 坝前水深Hm=78.8m =22.5解得hm=1.31381mT=4.438hm=5.08sL= =40.

14、4014m波浪爬高设计波浪爬高值应根据工程等级确定，2级坝采用累积频率为1%的爬高值.正向来波在单坡上的平均波浪爬高可按下式或有关规定计算:m=2.5= 0.771.0 / =2.083m查规范不同累积频率下的爬高与平均爬高的比值为2.23 R/R=2.23R1%=2.232.083=4.646m因风向与坝轴线的法线成22.5,波浪爬高应按等于正向来波计算爬高值乘以折减系数K查规范得0.95KR1%=0.95 4.646=4.41m坝顶在水库静水位以上的超高按下式确定：y=R+e+A=4.41+0.02017+1.00=5.43m坝顶高程2828.8+5.43=2834.2m坝高2834.2-

15、2750=84.2m为预防坝体竣工后的沉陷,预留0.4%的坝高坝高+0.4%沉陷=2834.2+84.20.004=2834.57m（4）正常蓄水位加非常运用条件下的坝顶超高正常蓄水位 2828.8m 吹程D=15 风速W=14m/s 坝前水深Hm=78.8m =22.5解得hm=0.689mT=4.438hm=3.6835sL= =21.18409m波浪爬高设计波浪爬高值应根据工程等级确定，2级坝采用累积频率为1%的爬高值.正向来波在单坡上的平均波浪爬高可按下式或有关规定计算:m=2.5= 0.771.0 / =1.093m查规范不同累积频率下的爬高与平均爬高的比值为2.23 R/R=2.2

16、3R1%=2.231.093=2.44m因风向与坝轴线的法线成22.5,波浪爬高应按等于正向来波计算爬高值乘以折减系数K查规范得0.95KR1%=0.95 2.44=2.315m坝顶在水库静水位以上的超高按下式确定：y=R+e+A=2.315+0.00632+0.50=2.821m坝顶高程2828.8+2.821=2831.62m坝高2831.62-2750=81.62m为预防坝体竣工后的沉陷,预留0.4%的坝高坝高+0.4%沉陷=2831.62+81.620.004=2832m坝顶高程计算成果表：表2-1 坝顶高程计算成果表计算情况计算项目正常运用情况非常运用情况上游静水位（m）正常蓄水位设

17、计洪水位正常蓄水位校核洪水位2828.82830.2752828.82831.425河底高程 （m）27502750坝前水深H（m）78.880.27578.881.425吹程D（km）1515风向与坝轴线法向的夹角 ()22.522.5风浪引起坝前壅高 e(m)0.020170.01980.006320.00612风速V（m/s）2514平均波高hm）1.313811.314230.6890.689护坡粗糙系数0.770.77上游坝面坡角cot1/2.5co1/2.5波浪沿坝坡爬高（m）4.414.422.3152.315安全超高A（m）1.00.5坝顶高程（m）2834.22835.728

18、31.622834.25坝顶高程加0.4沉陷（m）2834.57283628322834.6坝顶高程取为2836 m。第三节 土石料的设计3.1粘土料的设计3.1.1计算公式粘壤土用南京水利科学研究所标准击实试验求最大干容重、最优含水量。应该使土样最优含水量接近其塑限含水量，据此确定击数，得出多组平均最大干容重gdmax 和平均最优含水量W，设计干容重为gd设计填筑干容重gdmax 标准击实试验最大干容重m 施工条件系数，或称压实度，m值对于一二级坝或高坝采用0.960.99，三四级坝或低坝采用0.930.96。本设计的m=0.98。设计最优含水量为 = 3.1.2计算结果粘性土料设计的计算成

19、果见表3-1。表3-1 粘性土料设计成果表料场比重S最大密度gdma（g/cm3 ）设计干密度（g/cm3）塑限含水量p%流限%塑性指数Ip填筑含水量%自然含水量%最优含水量%1#下2.671.601.56823.1442.6019.4622.0724.822.072#下2.671.651.61722.2043.9021.7021.0224.221.021#上2.651.561.52925.0049.5724.5722.3025.622.302#上2.741.541.50926.3049.9023.5023.8026.323.803#下2.701.801.76420.0034.0014.001

20、6.9015.916.90料场孔隙比e内摩擦角渗透系数10-6cm/s有机含量灼热法%可溶盐含量%凝聚力kPa固结压缩系数cm2/kg1#下0.73424.674.321.730.07024.00.0212#下0.72125.504.801.900.01923.00.0201#上0.99023.171.902.200.11025.00.0262#上1.09321.503.960.250.11038.00.0333#下0.58028.003.001.900.08017.00.0103.1.3土料的选用上下游共有5个粘土料场，储量丰富。因地理位置不同，各料场的物理性质，力学性质和化学性质也存在一定

21、差异，土料采用以“近而好”为原则。规范指出粘土的渗透系数小于1010-6 cm/s,所有料场均满足要求。规范指出塑性指数大于20和液限大于40%的冲击粘土不宜作为坝的防渗体填筑材料，2#上 和2#下 料场因此不满足要求不予选用。1#上 料场土料的有机质含量为2.20%超过了规范规定的2%也不予采用。3#下与 1#下相比，3#下的渗透系数较小，更有利于防渗，设计干密度较大，压缩性能更好。 因此采用3#下料场作为主料场，1#下料场作为辅助及备用料场。3.2 坝壳砂砾料设计3.2.1计算公式 坝壳砂砾料填筑的设计指标以相对密实度表示如下：Dr =(emax-e)/(emax-emin) 或 Dr =

22、(rd-rmin)rmax/(rmax-rmin)rd e=n/(1-n)式中：emax为最大孔隙比；emin为最小孔隙比；e为填筑的砂、砂卵石或地基原状砂、砂卵石的孔隙比；rd为填筑的砂、砂卵石或原状砂、砂卵石干容重。设计相对密实度Dr要求不低于0.700.75；地震情况下，浸润线以下土体按设计烈度大小Dr不低于0.750.80。3.2.2计算结果砂砾料的计算成果见表3-24。表3-2 砂砾料的计算成果 料场不均匀系数5mm砾石含量（%）比重s天然孔隙比相对密实度Dr设计干容重gd(g/cm3)设计孔隙比e1#上48.853.02.750.4810.7321.860.4802#上44.053

23、.02.740.5310.6121.790.5303#上25.048.02.760.4490.8091.910.4504#上38.854.82.750.4600.7841.880.4601#下48.955.02.750.4810.7321.860.4802#下39.358.02.730.4750.7481.850.4753#下63.857.02.730.4810.7321.840.4804#下45.055.02.720.5110.66200.510料场保持含水量（%）湿容重(g/cm3)浮容重(g/cm3)内摩擦角粘聚力ckPa渗透系数K（10-2cm/s）1#上51.901.18351002

24、.02#上51.821.14360002.03#上51.951.21354002.04#上51.941.20363002.01#下51.901.18352002.02#下51.891.17364002.03#下51.881.17355002.04#下51.831.14371002.0各沙石料场的级配曲线如下：图3-1 各沙石料场的级配曲线3.2.3砂砾料的选用土石坝的坝壳材料主要为了保持坝体的稳定性，要求有较高的强度。根据规范要求砂砾石的相对密实度Dr要求不低于0.75，上述料场中只有3#上 和4#上 料场满足要求，其余料场不符合要求不予采用。下游坝壳水下部位和上游坝壳水位变动区宜有较高的透水

25、性，且具有抗渗和抗震稳定性，应优先选用不均匀和连续级配的砂石料。认为不均匀系数=30100时较易压实，510时则压实性能不好。而3#上 料场砂砾料的不均匀系数=2530,不满足要求，从颗粒级配曲线上也可以看出，4#上料场的砂砾料颗粒级配较好，物理力学指标也较高，可优先使用，故选择4#上 料场作为砂砾料的主料场。第四节 渗流计算4.1 渗流计算应包括以下内容:1. 确定坝体浸润线及其下游出逸点的位置，绘制坝体及其坝基内的等势线分布图或流网图；2. 确定坝体与坝基的渗流量；3. 确定坝坡出逸段与下游坝基表面的出逸比降，以及不同土层之间的渗透比降；4. 确定库水位降落时上游坝坡内的浸润线位置空隙压力

26、；5. 确定坝肩的的等势线、渗流量和渗透比降；4.2 渗流计算应包括以下水位组合情况：上游正常蓄水位与下游相应的最低水位上游设计水位与下游相应的水位3. 上游校核水位与下游相应水位4.2.1 计算方法选择水力学方法解决土坝渗流问题。根据坝体内部各部分渗流状况的特点，将坝体分为若干段，应用达西定律近似解土坝渗流问题。计算中假定任一铅直过水断面内各点的渗透坡降相等。计算简图见图6-41。通过防渗体流量： q=K(H2-H12)/2 + K2(H-H1)T/D 通过防渗体后的流量： q=K1(H12-T12)/2L1 + KT(H1-T1)T/(L+0.44T)假设：1）不考虑防渗体上游侧坝壳损耗水

27、头的作用；2）由于砂砾料渗透系数较大，防渗体又损耗了大部分水头，逸出水与下游水位相差不是很大，认为不会形成逸出高度；3）对于岸坡断面，下游水位在坝底以下，水流从上往下流时由于横向落差，此时实际上不为平面渗流，但计算仍按平面渗流计算，近似认为下游水位为零。由于河床冲积层的作用，岸坡实际不会形成逸出点，计算时假定浸润线末端即为坝趾。4.2.2 计算断面与计算情况对河床中间断面-以及左右岸坡段各一断面-、-三个典型断面进行渗流计算，计算按正常蓄水，设计洪水和校核洪水三种情情况进行。-断面：图4-1 I-I剖面图 （1）正常蓄水位与相应下游水位正常蓄水位2828.8m 下游相应的最低水位为2752.2

28、mK=3.0cm/s, K2=1.0 cm/s K1=2.010-2 cm/s KT取2.010-2 cm/s T1=2.2m T=32m D=0.9m H=78.8m =21.8mL1=189.275m L=196.875m=+ =由= 计算得： = 6.98, =2.413m（2）设计水位与下游相应的水位H=80.275m T1=4.83m L1=193.22m L=196.875m=+ =由=计算得： = 6.878, =5.03m（3）校核水位与下游相应水位H=81.425m T1=5m L1=193.475m L=196.875m=+ =由=计算得： = 7.2534, =5.20m

29、-断面：图4-2 II-II剖面图上游坝底高程取在2752m 处，下游坝底高程取在2805m 处，经简化近似等效为上下游坝底高程均在2790m，坝与地基的接触面近似为一水平面，覆盖层深28m。（1）正常蓄水位与相应下游水位K=3.0cm/s, K2=1.0 cm/s K1=2.010-2 cm/s KT取2.010-2 cm/s T1=0m T=28m D=0.9m H=38.8m =13.8mL1=108.1m L=115m=+ =由=计算得： = 2.84, =0.065m（2）设计水位与下游相应的水位H=40.275m L1=108.1m L=115m=+ =由=计算得： = 3.07,

30、 =0.0684m（3）校核水位与下游相应水位H=41.425m L1=108.1m L=115m=+ =由=计算得： = 3.15, =0.072m-断面：图4-3 III-III剖面图上游坝底高程取在2790m 处，下游坝底高程取在2765m 处，经简化近似等效为上下游坝底高程均在2775m，坝与地基的接触面近似为一水平面，覆盖层深18m。（1）正常蓄水位与相应下游水位K=3.0cm/s, K2=1.0 cm/s K1=2.010-2 cm/s KT取2.010-2 cm/s T1=0m T=18m D=0.9m H=53.8m =16.8mL1=147.35m L=155.75m=+ =

31、由= 计算得： = 3.66, =0.165m（2）设计水位与下游相应的水位H=55.275m L1=147.35m L=155.75m=+ =由= 计算得： = 3.93, =0.173m（3）校核水位与下游相应水位H=56.425m L1=147.35m L=155.75m=+ =由= 计算得： = 3.97, =0.1794m结果整理如下：表4-1 渗流计算结果汇总表计算情况计算项目正常蓄水位设计洪水位校核洪水位上游水深（m）-78.880.27581.425-38.840.27541.425-53.855.27556.425下游水深（m）-2.24.835-000-000逸出水深（m）

32、-2.4135.035.20-0.0650.0680.072-0.1650.1730.179渗流量（10-6 m3/sm）-6.986.887.25-2.843.073.15-3.663.933.97总渗流量Q（m3/d）155.52159.24165.024.3 逸出点坡降计算用流网法可准确求出渗流区任一点的渗透压力、渗透坡降和渗流量。绘制等势线和流网时，两者互相正交，并且分割的网格大致长宽相等。以I-I断面为例，各种工况下渗流逸出点坡降计算成过如表4-2表4-2 逸出点渗透坡降成果表断 面I-I工况正常蓄水位设计洪水位校核洪水位坡降 J1.781.501.60-断面（2750高程）正常蓄水

33、位情况:h=76.387/10=7.6387, l=4.3m,i=h/l=7.6387/4.3=1.78设计水位情况:h=75.245/10=7.5245, l=5m,i=h/l=7.5245/5=1.50图4-4 校核洪水位情况：h=76.225/10=7.6225，l=4.78m ，i=h/l=7.6225/4.78m =1.60第五节 稳定分析土石坝在自重、水荷载、渗透压力和地震荷载等作用下，若剖面尺寸不当或坝体、坝基土料的抗剪强度不足，坝体或坝体连同坝基有可能发生失稳。本设计中的稳定分析，主要指边坡的抗滑稳定。稳定分析计算的目的在于分析坝坡在各种不同工作条件下可能产生的失稳形式，校核其

34、稳定性。5.1 计算方法采用折线滑动法，认为折点在水位附近，并假设滑动面只在坝壳中，而防渗体不连同坝壳一起滑动，如图所示。滑动面上的抗剪强度利用充分程度应该是一样的，其安全系数表达方式为： tgi = tg(i)/KC式中：KC即为安全系数，（i）为试验得到的抗剪强度指标。图5-1 非粘性土坡稳定计算图计算中把滑动土体ABC分成两个滑块，P为滑面BD上的作用力。计算时先假设不同的安全系数KC，从滑块BCD平衡条件中求出P，然后将其作用在ABD上，看它是否与G2、R2等平衡，如不平衡则重新假定安全系数，重复上述步骤，直到平衡为止，此时的KC即为该滑动面的安全系数。在上述计算中、是任意假定的。设计

35、中为简化程序起见，还进行了以下假定：1）滑动面的折点在水位处；2）滑面BD铅直；3）将上下游变坡等效成一均匀坡。这些假定会对计算结果的精确度产生一定影响，但总体影响不是太大，近似计算中可以忽略。5.2 计算程序流程稳定计算的程序如下：Option ExplicitPrivate Sub Command1_Click()Dim g1, g2, p1, p2, , , As DoubleDim x1, y1, x2, y2, x3, y3, x4, y4, x5, y5, x6, y6, gan, hu, h1, xh1, xh2, xh3, xa, ya, xb, yb, xc, yc As S

36、ingleDim xd, yd, xd2, xb0, kk(), kc, k, ppp, temp, xa1, xa2 As SingleDim i, j, ii, jjj, m, n, iii As IntegerDim sss1, sss2, sss3, sss4, p, pp, pd1 As Double, pd2 As DoubleDim f1, f2, f3, fi1, fi2, fi3, 1, 2, 3, kh1 As Double, kh2 As DoubleDim xxa, yya, xxb, yyb, xxc, yyc As Singlex1 = 0: y1 = 0 定义坐标

37、，并输入土坝轮廓x2 = 225.2: y2 = 0x3 = 135: y3 = 45x4 = 137: y4 = 45x5 = 239.5: y5 = 86x6 = 242.4: y6 = 86gan = 19.11 = 0.63745166: 2 = 3 = 1h1 = Val(Text1.Text)If h1 = 45 Then temp = MsgBox(请输入合法水位！水位必须大于25，小于75.8)If Option2.Value = True Thenhu = 19.17Elsehu = 11.988End IfFor i = 0 To 27 选定一个 滑块终止A 点 xa =

38、x1 + 5 * i - 0.1 ya = xa / 3For j = 0 To 6 选定一个滑块起始C点 xc = x5 + 0.8 * j yc = y5 For m = 0 To 25 选定一个滑块转折B点 xb = 137.1 + 4 * m For n = 1 To 21 yb = ya + 4 * n If ya 25 Then If (yb - ya) / (xb - xa) (y4 - ya) / (x4 - xa) Or (y6 - yb) / (x6 - xb) (y6 - y2) / (x6 - x2) Or (yb = y5) Then GoTo 10 Else If

39、(yb - ya) / (xb - xa) (y5 - y4) / (x5 - x4) Or (y6 - yb) / (x6 - xb) (y6 - y2) / (x6 - x2) Or (yb = y5) Then GoTo 10 End If xd = xb If (xd = x5) Then yd = (y5 - y4) * (xd - 24.5) / (x5 - x4) Else yd = 86 End If = Atn(yb - ya) / (xb - xa) 求角度 = Atn(yc - yb) / (xc - xb) = 0 If (yb h1) And (yd = h1) Then 计算两个土块的相互作用力P1，P2 xh3 = (xc - xb) * (h1 - yb) / (yc - yb) + xb xh1 = 137 + 2.5 * (h1 - 45) xa2 = xb xa1 = 24.5 + (x4 - 24.5) * ya / (y4 - 0) sss1 = (xc - x5 + xh3 - xh1) * (y5 - h1) / 2 xd2 = (xc - xb) * (yd - yb) / (yc - yb) + xb sss2 = (xd2 - xd + xh3 - xh1) * (h1 - yd) / 2 + (xd2 - xd) * (yd