计算水力学课程论文.docxVIP

.定解问题 (1)基本方程：圣维南方程组 ［ 空+3=ql(连续方程) TOC \o 1-5 \h \z :t : x 2 -r —+ —(———)+gA J +gASf =0(动量方程) \o Current Document Ft ；x A ；x (2)定解条件：初始条件和边界条件 条件1:水位边界条件已知 j 上游边界：Z=4.5( m) t 下游边界：Z=4.0+1.5sin(—)(m) 12 条件2:流量边界条件已知 ［上游边界：Q=500(m3/s) t ［下游边界：Z=4.0+1.5sin(— )(m) 条件3:水位一流量关系条件已知 ‘ 上游边界：湖泊A=20km2 下游边界：Z=4.0+1.5sin(^1)(m) .差分格式 用简化的四点线性隐格式方法的差分方程对圣维南方程组进行离散： :Qj书—Qj +CjZj++gZj=Dj 、EjQj +GjQj 书+匕乙书―FjZj =4 .计算方法 追赶法：边界条件的追赶关系 条件1：水位边界条件已知 Qj = Sj 1 -Tj 1Qj 1 Zj 1 = Pj 1 q 1 条件2:流量边界条件已知 ;Zj=Sj+-Tj^Zj^ @书=耳由-Vj? + 条件3:水位一流量关系条件已知 可线性化处理后，按照条件 2方法进行计算。 4.程序框图开始? 上边界条件计算数组说明计算追赶系数数组基本数据下边界条件计算单位换算回代计算水位流量 4.程序框图 开始 ? 上边界条件计算 数组说明 计算追赶系数 数组基本数据 下边界条件计算 单位换算 回代计算水位流量 初始条件赋值 输出计算结果 .成果分析 基于教科书附录中的 FORTRAN代码，自主转换编译为 Visual C#程序进行成果分析。 程序说明 (1)红色曲线代表流量过程线(48h),黑色代表水位流量过程线(48h),通过对各参 数进行设置，自动绘出过程线，并能将数据结果输出到 TXT文档。 (2)因部分断面流量起伏较大或水位起伏较小，致使曲线无法显示或不便进行观察比 较，故可通过乘以缩放倍比对曲线进行调整(同一条件下选择最优缩放倍比) 。 (3)本程序的重要优势在于可对计算中的重要参数进行迅速调整，简单明了模拟出其 对河道水流计算的影响程度和作用情况。 (4)本程序无法绘制同一条件下不同断面流量(或水位)过程线，因而给断面情况对 比分析带来一定的不便。 编译结果 编译结果界面如下图所示（根据题意，以断面 1、6、11、16、21为例）： （1）条件1结果（水位缩放倍比 50,流量缩放倍比0.5）: 断面1: 时向 断面挪 计算d用射门 _时算量件ic 河较长nJ 时间步恸r 权重翘仃 恨皿 21 48 1 1.0- 10?0 0~75 0.02~ 绘制画面雌二 ［也m ］二出颉函而「~|水~fc丽~ 融籥放0.5 断面6: 断面 PAGE 断面 PAGE #: 断面11: 断面21: (2)条件2结果(水位缩放倍比 25,流量缩放倍比0.15): 断面：1 : 断面 PAGE 断面 PAGE #: 断面6: 断面 PAGE 断面 PAGE #: 断面16: 断面 断面6: (2)条件3结果(水位缩放倍比 50,流量缩放倍比0.15): 断面1: 断面 断面 PAGE #: 断面11: 断面21: 规律分析 (1)水位过程线分析：条件 1中越往下游，水位不断上涨，最高水位缓慢左移；随时 间推移，逐渐由水平线转为稳定的三角函数曲线。 条件2中上中游水位过程线形状与上游流 量过程线类似，下游形状与三角函数曲线相同。条件 3中越往下游，水位不断上涨，最高水 位增幅较大且缓慢左移；全过程的三角函数曲线形状由扁平逐渐趋于稳定。 (2)流量过程线分析：条件 1中越往下游，曲线缓慢左移(洪峰提前) ，各断面洪峰值 变幅较小；全过程基本呈三角函数状。条件2中上游断面开始时陡然提升， 随后流量趋于平 缓；中下游断面则呈逐渐上升趋势，且越往下游，增幅越大，三角函数曲线形状不明显。条 件3中越往下游，曲线缓慢左移(洪峰提前) ，各断面洪峰值变幅较小；全过程曲线形状与 水位类似，但变幅更大。 (3)整体分析：依据洪水的叠加效应，越往下游洪量越集中，洪峰有提前的趋向。 (4)稳定性分析：本题中的 8值取0.75, Preissmann格式的稳定条件是无条件稳定， 存在着数值弥散。而从图形中也正验证了这一点， 即全过程流量、水位曲线形状与下游初边 DT *C 值条件(或原修分方程斛的性质)有关。此外，由于 == 因而可考虑通过对时间步 DX 长DT和河段步长DX的搭配，进一步改善格式的稳定性， 保证相位误差控制在一定范围内， 也正因为如此，本题中 DT取10min , DX取1km (可取0-1.6km)。 (5)在对