本程序采用了特征线法计算管道阀门关闭时的水锤
1