河流模拟课程设计—水库一维泥沙-淤积计算(共29页).docx

《河流模拟课程设计—水库一维泥沙-淤积计算(共29页).docx》由会员分享，可在线阅读，更多相关《河流模拟课程设计—水库一维泥沙-淤积计算(共29页).docx（29页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、精选优质文档-倾情为你奉上水库一维泥沙淤积计算课程设计武汉大学水利水电学院2013-3-15目录专心-专注-专业一、 目的与要求通过课程设计，初步掌握一维数学模型建立数学模型的基本过程和计算方法，具备一定的解决实际问题的能力。以水流、泥沙方程为基础，构建恒定流条件下的河道一维水沙数学模型，并编制出完整的计算程序，并以某个水库为实例，进行水库泥沙淤积计算。水流条件：恒定非均匀流。泥沙条件：包括悬移质，推移质的均匀沙模型，推移质计算模式为饱和输沙，悬移质计算模式为不饱和输沙，水流泥沙方程采用非耦合解。二、 基本原理1、 基本方程水流连续方程： 水流运动方程 或 泥沙连续方程 河床变形方程 推移质平

2、衡输沙方程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、输入河床地形

3、糙率等数据求得断面面积与水位的关系（AZ），进而求得断面平均流速 ，水力学半径 2、读入一个时段的水沙数据（特别注意，不要一次性将数据全部读入）读入第一时段（Q,S）值3、计算水面线，同时得到各断面的水力要素求得各个断面的河宽、断面面积、水深、平均流速等值计算前要注意在坝前输入水位，各断面均应对流量赋值4、计算悬移质水流挟沙力K取 0.124，m取1.05。5、计算推移质输沙率（采用mayer-peter公式）6、计算各断面含沙量公式7、 计算各断面冲淤厚度进口断面 8、 修改水各断面水下河床高程9、重新进入（2）进行下一循环10、计算10年河床变形，计算时段为一天，单位为秒（s）11、淤积总

4、量年输出一次，其余每两年输出一次计算结果四、 计算框图开始读入地形资料和糙率输出初始库容、深泓等i=i+1读入一个QYQ=0计算水流子程序计算S*和Gbj=j+1计算S(j)和dy(j)YDy(i)=0,重新计算S(j) 判断库尾冲刷与否YJnpnxt修改河床地形判断年份输出结果YIx(m,2,i)thenalow0(i)=x(m,2,i)endifenddo enddoend!*泥沙沉降速度的计算,采用张瑞瑾公式*FUNCTION FW(T,gama,gamas,d,ndisp) IF(NDISP.EQ.-1)WRITE(*,*)INTO FWcall VISCOS(T,CMU)A1=CMU

5、*13.95/dA2=1.09*9.8*d*(gamas-gama)/gamaFW=(A1*2.0+A2)*0.5-A1 RETURN ENDSUBROUTINE VISCOS(T,CMU) X=1.775E-06 A=1+0.0337*T+0.*T*T CMU=X/A RETURN END!*水面线及各断面水力要素计算函数*subroutine level(x,rough,npxt,zlevel,dx,q,npoint,b,a,xw,nn,mm,failev,NDISP) dimension x(mm,2,nn),rough(npxt),dx(npxt),zlevel(npxt) dimen

6、sion q(npxt),npoint(npxt),b(npxt),a(npxt),xw(npxt) IF(NDISP=-1)WRITE(*,*)INTO LEVEL nc=1000 dz=0.5 dz1=0.1317 dz2=2.079 call area(npoint(npxt),x(1,1,npxt),x(1,2,npxt),b(npxt),a(npxt),xw(npxt),zlevel(npxt),ndisp) do ip=npxt-1,1,-1NR=0 zmin=zlevel(ip+1)+dz30 call area(npoint(ip),x(1,1,ip),x(1,2,ip),b(ip),a(ip),xw(npxt),zmin,ndisp) if(A(ip)0)thenfr=(q(ip)/a(ip)*2.0*b(ip)/(9.8*a(ip) if(fr0)thenzmax=zmax+dzgoto 20endifelsezmin=zmin+dz2NR=NR+1IF(NRNC)THENzmin=zmin-dz1write(*,*)the loop is death,pauseWRITE(*,*)nr,ip,zmin,fmin,nr,ip,zmin,fminread(*,*)endifgoto 30endif elsezmin=zmin-dz1