悬链线形断面渠道正常水深的简化算法

2012-11-25滕凯

水科学与工程技术 2012年6期

滕凯

（齐齐哈尔市水务局，黑龙江齐齐哈尔161006）

TENG Kai

（Qiqihar City Water Authority，Qiqihar 161006，China）

随着现代化施工机械连续施工作业技术的发展，具有受力条件好、水流条件优越的悬链线形断面渠道将越来越广泛地应用于水利水电供排水工程。由于悬链线形断面渠道正常水深计算涉及高次方程求解问题，常规的试算法、图表解法［1-3］、迭代法［4-5］或是计算过程繁复、依赖图表，或是求解成果精度不高，不便应用；而由计算机编程求解又不便实际工作，因此，寻求一种更加简单实用的简化计算方法具有一定的实际意义，但到目前为止，尚没有比较理想的有关悬链线形断面渠道正常水深计算方法的研究成果。为有效简化该种断面正常水深的求解计算方法，本文在对该断面正常水深计算公式变形整理的基础上，通过引入无量纲水深参数，采用优化拟合的方法，经逐次拟合逼近，获得了计算过程简捷、求解精度高、使用范围广的简化近似计算式，可在实际工程设计中推广应用。

1 正常水深的基本计算方程

根据水力学原理［6］，正常水深的基本计算方程为：

式中 Q为渠道通过流量（m3/s）；X为湿周（m）；A为过水断面面积（m2）；n为渠道糙率系数；i为渠道设计坡降。

悬链线形断面渠道如图1所示。

图1 悬链线形渠道过水断面

断面的曲线方程为：

设水面宽度为B，则其过水断面水力要素为：

式中 a为悬链线形断面形状参数。

由式（3）可得：

将式（4）～（7）分别代入式（1），并设

式中 x为无量纲收缩水深；h0为渠道正常水深（m）；β 为渠道宽深比系数；k为无量纲中间变量参数。

经对式（1）方程组变形整理，即可获得计算悬链线形断面渠道正常水深的基本公式：

由式（9）可见，在理论上x 的值域范围为（0，∞）。而在实际工程中，根据文献［7］的研究结果，悬链线形断面的最优断面形式为B/a=3.21223，将其代入式（6），可求得x=1.5921，因此，在实际工程中，x的值域范围不会超出［0.06，29.0］，相应的宽深比系数β 的值域范围为［0.282，11.49］，相应k的值域范围为［0.0209，23.8211］。

2 正常水深简化公式及精度分析

2.1 简化公式的建立

式（9）为超越方程求解问题，无法通过解析法获解。据此，本文采用拟合替代形式提出简化计算方法。

依据上述值域范围，利用式（9）完成k～x关系曲线绘制（曲线图略），并假定函数x=f（k）可以拟合k～x关系曲线并替代式（9），经数值替代分析及优化拟合计算，以标准剩余差最小为目标函数［8］，即

式中 x′为第i（i=1，2，3，…，n）拟合替代值；n为拟合计算的数组数。

经逐次逼近拟合［9］，即可获得式（9）的最优拟合替代函数式：

由式（10）求出x后，即可用式（11）求得正常水深h0为:

2.2 简化公式的精度分析

为比较式（10）与式（9）的拟合精度，在给定的工程实用参数范围内（即0.06≤x≤29.0），选取不同的xi代入式（9）即可计算出与之相对应的ki，再将ki代入式（10），求得与之相对应的xi′（xi′为近似计算值），并由式（12）进行拟合相对误差计算。

式中 zi为拟合相对误差（%）；i为拟合计算的第i个数据比较（i=1，2，3，…，n）。

计算结果见表1。

由表1精度比较可见，在工程实用范围内，即0.06≤x≤29.0，0.0209≤k≤23.8211情况下，式（10）的最大拟合替代相对误差仅为0.257%，且有79%的拟合点的相对误差小于0.2%，可见，式（10）具有较好的拟合替代精度，完全可以满足实际工程的设计精度要求，且该公式在表达形式上较文献［5］更直观，具体比较见表2。