回顾
连续方程:
动量方程:
1.定义用于连续方程的蓄存宽度和用于运动方程的计算宽度B(B不太理解)
改写方程组(不理解为什么被省略了)
其中,C为谢才(chezy)系数,Rh为水力半径。
复习曼宁公式、谢才公式
管渠沿程水头损失常用谢才公式计算
R为过水断面水力半径,C为谢才系数,l为管渠长度
曼宁引入粗糙系数n,计算谢才系数C,适用于明渠、非满管流或较粗糙的管道水力计算。
代入谢才公式
2.链式法则
其中
引入计算宽度B,令A=Bh
代入运功方程就写成
3.离散方式
在每个节点按顺序交错算出Q或h
连续方程
连续方程的差分形式(应写作)
运动方程
运动方程的差分形式
以上D1、D2、D3、D4都取t=(n+1/2)Δt的值
4.简化方程组,写作Ax+By+Cz=D的形式
上标表示t,下标表示x,假设时间步长均匀有单位步长Δt,空间步长(节点距离)不均匀。
连续方程
动量方程
个人理解应为
5.假设Q、h在空间上顺序交错求解
代入上一步的简化方程
连续方程
运动方程
形式同假设,因此
,
,
6.边界条件
假设x=J处流量Q给定,则最后一个方程是
可得,
又
可得,
两式系数相比较,可得
,
类似的,如果边界x=J处水深给定,则
,
7.追赶法
假设x=J处边界条件Q给定,即已知
根据,,可得,
又根据4.简化方程组-连续方程,可得,,,但(不会求)
根据 , ,可得
又根据4.简化方程组-运动方程,D1、D2、D3、D4都取t=(n+1/2)Δt的值(不会求),再求得,但(不会求)
根据 , ,可得
假设另一边界为h,即已知
追的过程:→,→→→......→
根据 ,
赶的过程: →→→......→