河流模拟课程设计—水库一维泥沙_淤积计算

《河流模拟课程设计—水库一维泥沙_淤积计算》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、一维泥沙淤积计算水库一维泥沙淤积计算课程设计武汉大学水利水电学院2013-3-15一维泥沙淤积计算目录一、目的与要求1二、基本原理11、基本方程12、方程离散13、公式补充2三、计算步骤3四、计算框图4五、计算结果51、历年输沙量特征值52、各年淤积总量53、各年水位库容关系64、水面线的变化75、深泓变化86、坝前断面变化9六、结果分析121、剖面形态分析122、库容损失合理性分析12七、计算程序13一维泥沙淤积计算一、目的与要求通过课程设计，初步掌握一维数学模型建立数学模型的基本过程和计算方法，具备一定的解决实际问题的能力。以水流、泥沙方程为基础，构建恒定流条件下的河道一维

2、水沙数学模型，并编制出完整的计算程序，并以某个水库为实例，进行水库泥沙淤积计算。水流条件：恒定非均匀流。泥沙条件：包括悬移质，推移质的均匀沙模型，推移质计算模式为饱和输沙，悬移质计算模式为不饱和输沙，水流泥沙方程采用非耦合解。二、基本原理1、基本方程水流连续方程：①水流运动方程②或③泥沙连续方程④河床变形方程⑤推移质平衡输沙方程G=G*⑥水流挟沙力公式采用张瑞瑾公式，推移质输沙率公式采用Mayer-_Peter公式，MAYER-PETER公式中的能坡J按均匀流曼宁公式近似计算（每个断面不同）。2、方程离散方程①在恒定流情况下有，离散为:Q=const方程③变形为~27~一维泥沙

3、淤积计算或上式离散为方程（4）去掉时间项得到该方程的解析解为：由方程（4-5）可得对2号断面以下，上式可以离散为：对于进口断面，推移质不考虑，悬移质采用单点离散方程（5）可离散为：1、公式补充K取0.124，m取1.05，干密度取1.3恢复饱和系数均匀沙粒径为d=0.041mm（悬移质），d=2mm（推移质）~27~一维泥沙淤积计算一、计算步骤1、输入河床地形糙率等数据求得断面面积与水位的关系（AZ），进而求得断面平均流速，水力学半径2、读入一个时段的水沙数据　（特别注意，不要一次性将数据全部读入）读入第一时段（Q,S）值3、计算水面线，同时得到各断面的水力要素求得各个断面的河

4、宽、断面面积、水深、平均流速等值计算前要注意在坝前输入水位，各断面均应对流量赋值4、计算悬移质水流挟沙力K取0.124，m取1.05。5、计算推移质输沙率（采用mayer-peter公式）6、计算各断面含沙量公式7、计算各断面冲淤厚度进口断面8、修改水各断面水下河床高程9、重新进入（2）进行下一循环10、计算10年河床变形，计算时段为一天，单位为秒（s）11、淤积总量年输出一次，其余每两年输出一次计算结果~27~一维泥沙淤积计算一、计算框图开始读入地形资料和糙率输出初始库容、深泓等i=i+1读入一个QYQ=0计算水流子程序计算S*和Gbj=j+1计算S(j)和dy(j)YDy(

5、i)=0,重新计算S(j)判断库尾冲刷与否YJ

6、990438565.7818001991597289.5818701992198229.7312901993645696.841900最大值645696.841900最小值198229.731290多年平均值428764.3116402、各年淤积总量表2历年淤积总量的变化年悬移质输沙量(万吨)年均流量(万吨)总淤积量(亿立方米)000011821.490.1425291.2050.387316993.540.879420027.341.085521668.331.204623853.981.368730155.051.701834879.591.987936754.182.11

7、71039289.662.243~27~一维泥沙淤积计算图1历年淤积总量的变化3、各年水位库容关系表3各年水位库容关系水位（m）0年2年4年6年8年10年2732.5742.2071.5721.3130.7650.5692712.3852.0181.3841.1320.6160.4382692.2021.8351.2020.980.4930.3332672.0241.6611.0440.8590.3980.262651.8521.510.9280.7730.3460.232631.6881.