计算水力学-单一河道程序设计.ppt

计算水力学教学课件计算水力学教学课件单一河道程序设计单一河道程序设计题 目：计算程序框图计算追赶系数上边界条件计算下边界条件计算回代计算水位流量输出计算结果结 束数组说明开 始输入基本数据单位换算初始条件赋值按按时时间间步步长长循循环环1.数组说明数组说明|水位、流量；河段长；|P、V、S、T；糙率。123NN-1 2.输入基本数据输入基本数据计算断面数、计算小时数、计算条件、河段长度、时步长、权重系数、计算糙率123NN-1 3.单位换算单位换算时间单位：秒长度单位：米123NN-1 N-14.初始条件赋值初始条件赋值|初始水位条件：Z（L+1）=水面线|初始流量条件：Q（L+1）=0.0123NN-1 ZZ(I)=4.0+0.025*(N-I)5.边界条件赋值边界条件赋值上边界条件上边界流量条件：123NN-1 5.边界条件赋值边界条件赋值下边界条件下边界水位条件：123NN-1 6.按时间步长循环按时间步长循环计算追赶系数上边界条件计算下边界条件计算回代计算水位流量输出计算结果按按时时间间步步长长循循环环6.按时间步长循环按时间步长循环其 中：简化形式：追追 赶：赶：|所以所以流量上边界条件流量上边界条件追赶：追赶：6.追赶系数追赶系数EXTERN子程序子程序1iN上边界条件P1=Q(t)V1=0Pi、Vi、Si、TiZi、QiAi、Bi、Ri、UiCi、Di、Ei、Fi、Gi、iY1i、Y2i、Y3i、Y4i ABR函数6.追赶系数追赶系数EXTERN子程序子程序12NN-1上边界条件P1V1P2V2S2T2iPiViSiTiPN-1VN-1SN-1TN-1PNVNSNTN6.追赶系数追赶系数ABR子程序子程序100 m31ZZD6.按时间步长循环按时间步长循环计算追赶系数上边界条件计算下边界条件计算回代计算水位流量输出计算结果按按时时间间步步长长循循环环7.回代计算回代计算BAEXT子程序子程序1iNN-1下边界条件Z1QupZNQN=PN-VNZNZN-1=SN-TNZNQN-1=PN-1-VN-1ZN-18.输出计算结果输出计算结果计算追赶系数上边界条件计算下边界条件计算回代计算水位流量输出计算结果123NN-1 按时间步长循环