抛物线形断面渠道收缩水深简易算法.docVIP

抛物线形断面渠道收缩水深的简易算法 滕凯 摘要:针对目前求解抛物线形断面渠道收缩水深普遍存在的计算繁复问题，通过引入无量纲相对水深参数，对抛物线形断面渠道收缩水深基本计算公式进行变形整理。在此基础上，以保证求解成果精度满足工程设计要求为前提，以获得最简化的计算公式为目的，对函数高次方程优化拟合，获得了表达形式简单、计算简捷、适用范围广的近似计算公式。具有实际应用推广价值。 关键词:抛物线形断面; 收缩水深; 优化拟合; 简化计算 中图法分类号: TV131．4 文献标志码: A 文章编号: 1001-4179(2013) 09-0097-03 A simple algorithm of contraction depth of parabola channel TENG Kai Abstract: The calculation of contraction depth of parabola channel is complicated．By introducing dimensionless relative depth parameter，and on the basis of reference and transformation of the basic contraction depth formula of parabola channel，we optimized and fitted the higher order equation obtaining the approximate formula with simple expression．The presented calculation method is convenient and has wide application range． Key words: parabolic cross-section; contraction depth; optimal fitting;simplified calculation 收缩水深是工程设计人员研究水工建筑物下游水跃类型及确定消能工形式的重要水力参数。因此，在水工建筑物的水力计算中，该参数的简化获解，对提高设计工作的实际效率具有重要意义。抛物线形断面渠道由于具有受力条件好、水流条件优越、断面连续、便于机械化施工作业等优点，正越来越广泛地应用于水利水电供排水工程。由于抛物线形断面渠道收缩水深计算涉及高次方程求解，常规的试算法、图表解法［1-3］、迭代法［4］，或是计算过程繁复、依赖图表，或是求解成果精度不高，不便应用; 而由计算机编程求解又不便于实际工作。为改进上述计算方法，文献［5-7］通过数学方法提出了相关的近似计算公式，有效简化了求解过程。但这些公式均存在一个共性问题，即相关公式仍显过于繁复。因此，寻求一种更加简单实用的简化计算方法具有一定的实际意义。 本文在对抛物线形断面渠道收缩水深计算公式变形整理的基础上，通过引入无量纲水深参数，在保证工程设计精度要求的前提下，采用优化拟合的方法，经逐次拟合逼近，获得了表达形式简单、计算过程简捷、使用范围广的简化近似计算式，可在实际工程设计中推广应用。 1 收缩水深基本计算公式 收缩水深的基本方程为［8］ （1） 式中，E0为以收缩断面底部为基准面的过水建筑物上游总水头; hc为收缩断面处的水深; Q 为过水流量; g为重力加速度; φ为流速系数; Ac为与收缩水深hc相对应的抛物线形断面面积。 抛物线形断面曲线见图1。 曲线方程为 y=px2　 p＞0 (2) 其过水断面面积为 （3） 式中，p 为抛物线形断面形状参数。 将式(3) 式代入式(1) 式，并设 （4） 式中，x 为无量纲收缩水深; β为综合已知参数。 经进一步整理可得 β=x(1-x)1/3 （5） 由式(5) 求解出x，再由式(4) 即可求得抛物线形断面渠道收缩水深hc。 2 收缩水深简化公式及精度分析 2．1 简化公式的建立 式(5) 为4 次方程求解问题，x 有精确解，但计算过程繁琐［9］，不便实际应用。为简化相关计算，考虑在实际工程中0≤x≤0．5［5 －7］，利用式(5) 展绘x ～β关系曲线图，并假设x ～ f(β) 函数(x 为近似值) 可以拟合x ～β关系曲线并替代式(5) ，经数值替代分析及优化拟合计算，以标准剩余差最小为目标函数［10］，即 式中，n 为拟合计算的组数。 经逐次逼近拟合［11］，即可获得式(5) 的最优拟合替代函数式为 x=β·e0.84β1.436 (6) x 求出后即可用下式求得收缩水深。 hc=xE0　 （7） 2．2 简化公式的精度分析 为比较式(6) 与式(5) 的拟合精度，在给定的工程实用参数范围内(即0≤ x≤0．5，0≤β