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

计算水力学教学课件计算水力学教学课件单一河道程序设计单一河道程序设计题 目：{{{计算程序框图计算追赶系数上边界条件计算下边界条件计算回代计算水位流量输出计算结果结 束数组说明开 始输入基本数据单位换算初始条件赋值按按时时间间步步长长循循环环§1. 数组说明数组说明| 水位、流量；河段长；| P、、V、、S、、T；糙率123NN-1… … …§2. 输入基本数据输入基本数据计算断面数、计算小时数、计算条件、河段长度、时步长、权重系数、计算糙率123NN-1… … …§3. 单位换算单位换算时间单位：秒长度单位：米123NN-1… … …N-1§4. 初始条件赋值初始条件赋值| 初始水位条件：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-1PNVNSNTN……§6. 追赶系数追赶系数—ABR子程序子程序100 m31ZZD§6. 按时间步长循环按时间步长循环计算追赶系数上边界条件计算下边界条件计算回代计算水位流量输出计算结果按按时时间间步步长长循循环环§7. 回代计算回代计算—BAEXT子程序子程序1iNN-1下边界条件Z1QupZNQN=PN-VNZNZN-1=SN-TNZNQN-1=PN-1-VN-1ZN-1……§8. 输出计算结果输出计算结果计算追赶系数上边界条件计算下边界条件计算回代计算水位流量输出计算结果123NN-1… … …按时间步长循环。