河流模擬课程设计稿—水库一维泥沙淤积计算.docxVIP

 目录 一、目的与要求 1 二、基本原理 1 1、基本方程 1 2、方程离散 1 3、公式补充 1 三、计算步骤 1 四、计算框图 1 五、计算结果 1 1、历年输沙量特征值 1 2、各年淤积总量 1 3、各年水位库容关系 1 4、水面线的变化 1 5、深泓变化 1 6、坝前断面变化 1 六、结果分析 1 1、剖面形态分析 1 2、库容损失合理性分析 1 七、计算程序 1 一、 目的与要求 通过课程设计，初步掌握一维数学模型建立数学模型的基本过程和计算方法，具备一定的解决实际问题的能力。 以水流、泥沙方程为基础，构建恒定流条件下的河道一维水沙数学模型，并编制出完整的计算程序，并以某个水库为实例，进行水库泥沙淤积计算。 水流条件：恒定非均匀流。 泥沙条件：包括悬移质，推移质的均匀沙模型，推移质计算模式为饱和输沙，悬移质计算模式为不饱和输沙，水流泥沙方程采用非耦合解。 二、 基本原理 1、 基本方程 水流连续方程： ① 水流运动方程 ② 或 ③ 泥沙连续方程 ④ 河床变形方程 ⑤ 推移质平衡输沙方程 G=G* ⑥ 水流挟沙力公式采用张瑞瑾公式，推移质输沙率公式采用Mayer-_Peter公式，MAYER-PETER公式中的能坡J按均匀流曼宁公式近似计算（每个断面不同）。 2、 方程离散 方程 ①在恒定流情况下有 ，离散为:Q=const 方程 ③变形为 或 上式离散为 方程（4）去掉时间项得到 该方程的解析解为： 由方程（4-5）可得 对2 号断面以下，上式可以离散为： 对于进口断面，推移质不考虑，悬移质采用单点离散 方程（5）可离散为： 3、 公式补充 K取 0.124，m取1.05，干密度取1.3 恢复饱和系数 均匀沙粒径为d=0.041mm（悬移质），d=2 mm（推移质） 三、 计算步骤 1、输入河床地形糙率等数据 求得断面面积与水位的关系（AZ），进而求得断面平均流速 ，水力学半径 2、读入一个时段的水沙数据　（特别注意，不要一次性将数据全部读入） 读入第一时段（Q,S）值 3、计算水面线，同时得到各断面的水力要素 求得各个断面的河宽、断面面积、水深、平均流速等值 计算前要注意在坝前输入水位，各断面均应对流量赋值 4、计算悬移质水流挟沙力 K取 0.124，m取1.05。 5、计算推移质输沙率（采用mayer-peter公式） 6、计算各断面含沙量公式 7、 计算各断面冲淤厚度 进口断面 8、 修改水各断面水下河床高程 9、重新进入（2）进行下一循环 10、计算10年河床变形，计算时段为一天，单位为秒（s） 11、淤积总量年输出一次，其余每两年输出一次计算结果 四、 计算框图 开始 读入地形资料和糙率 输出初始库容、深泓等 i=i+1 读入一个Q Q=0 计算水流子程序 计算S*和Gb j=j+1 计算S(j)和dy(j) 判断库尾冲刷与否 Jnpnxt 修改河床地形 判断年份输出结果 I3653 结束 五、 计算结果 1、历年输沙量特征值 计算各年悬疑质输沙量、推移质输沙量、年均流量见表1。（其中推移质输沙量是按照悬疑质输沙量的1.5%计算的）。 表1龙开口坝址历年流量输沙量特征值 代表年 年平均悬移质 输沙量 （万t） 年平均推移质 输沙量 （万t） 年平均流量 （万t） 1984 2941 44.12 1330 1985 4120 61.80 1750 1986 3221 48.32 1360 1987 4828 72.42 1760 1988 2543 38.15 1440 1989 6426 96.39 1900 1990 4385 65.78 1800 1991 5972 89.58 1870 1992 1982 29.73 1290 1993 6456 96.84 1900 最大值 6456 96.84 1900 最小值 1982 29.73 1290 多年平均值 4287 64.31 1640 2、各年淤积总量 表2 历年淤积总量的变化 年 悬移质输沙量(万吨) 年均流量(万吨) 总淤积量(亿立方米) 0 0 0 0 1 1821.49 484556 0.14 2 5291.205 1099418 0.387 3 16993.54 1965807 0.879 4 20027.34 2469775 1.085 5 21668.33 2765617 1.204