无压六圆弧蛋形断面隧洞正常水深的简化计算

2013-08-29滕凯

中国水能及电气化 2013年3期

滕凯

(齐齐哈尔市水务局，黑龙江省齐齐哈尔市 161006)

由于蛋形断面与其他断面形式 (马蹄形、门洞形、圆形及矩形等)相比具有受力条件好、过流能力强、适应地质环境能力突出、节省空间及材料显著等优点，因此，随着水利机械化施工技术及工艺的不断提高，该种断面正在越来越广泛地应用于水利水电供排水工程［1，2］。随着工程综合技术的进一步发展，该种断面必将以其突出的优越性在未来的输配水工程中发挥越来越重要的作用。

由于该种断面由六个圆心三个不同半径的六段连续的圆弧曲线构成，几何图形复杂，正常水深计算不但涉及四个区间的分段函数而且需完成超越方程求解，常规的解析法不能直接获解，试算法［3，4］由于公式分段、表达形式为繁复的超越方程，大量的重复计算不但过程繁索，且易出错、精度差、效率低，不便实际应用，利用微机编程求解又不便于实际工作。因此，研究提出一种方便实际工程应用的简化近似计算方法具有重要的意义，而将六圆弧蛋形断面正常水深计算的分段函数统一由一个近似计算公式表示并直接解算，到目前为止尚未见到相关研究成果。

为了解决常规算法及目前现有研究成果所存在的问题，本文依据优化拟合原理，取目标函数为标准剩余差最小，经逐次逼近拟合计算，获得了在工程实用范围内用一个计算通式直接完成求解计算的简化解析式，使求解计算过程更加简捷、直观，便于实际应用。

1 蛋形断面正常水深的基本计算公式

根据明渠均匀流基本方程［5］，六圆弧蛋形过水断面正常水深应按下式计算:

式中 Q——渠道通过流量，m3/s;

A——过水断面面积，m2;

χ——湿周，m;

i——隧洞底坡降;

n——隧洞内壁糙率系数。

1.1 蛋形断面水力要素

六圆弧蛋形过水断面 (见图1)是由三种不同长度的半径分别在六个圆心位置通过弧线构成的闭合曲线，其半径分别为r1、r2和r3。各圆弧段所对应的圆心角分别为 α1、α2、α3和 α4。其数值关系为:

图1 六圆弧蛋形断面

图2 六圆弧蛋形断面计算图

由图1的几何关系可以求得:h1=0.096478r1，h2=0.195651r1，其水力要素为［6］

以上式中 h——洞内水深，m;

φ、β、β1、β2——与 α1、α4、α3、α2同圆心的不同区段水深所对应的圆心角 (见图 2)，rad。

1.2 蛋形断面正常水深基本计算公式

将蛋形断面四个区段所对应的水力要素式式(2)～式 (4)分别代人式 (1)，并设

式中 k——已知综合参数;

h0——蛋形断面的正常水深，m;

x——无量纲正常水深，m。

经进一步整理即可获得求解六圆弧蛋形过水断面正常水深的基本公式为

式 (6)为对应自变量x四个区间的四个超越方程，利用式 (6)即可根据已知的k值采用相应区间的计算公式通过非解析法完成无量纲正常水深x的求解，进而由下式求得正常水深，即

由式 (6)可见，x的取值范围为四个连续区间，值域为［0，1］，而在实际工程中，x＜0.09情况极为少见;同时，考虑为避免洞内产生不稳定的半有压流现象，《水工隧洞设计规范》［7］规定，洞内水面线以上的富余空间面积不宜小于隧洞断面积的15%，且净空高度不应小于0.4m，由此可以求得x＜0.775，因此，本文将x的取值范围确定为［0.09，0.8］，相应参数k的取值范围为［0.012，0.360］。

2 蛋形断面正常水深的近似公式及精度分析

2.1 近似公式的建立

由式 (6)可见，该方程为含未知量x的不连续分段超越方程，无法直接获解。为解决这一复杂方程的求解问题，在工程实用参数范围内 (即0.09≤x≤0.8，0.012≤k≤0.360)，假定函数x=f(k)可以替代式 (6)，依据式 (6)即可通过设定不同的x求得k，并展绘x－k关系曲线 (限于篇幅，本文曲线图略)，根据曲线图形关系，经多组备选函数的数值拟合分析及回归计算，以标准剩余差最小为最终希望目标函数［8-15］，即

n——拟合计算的数组数。经逐次逼近拟合［16，17］即可获得如下替代函数，即

利用式 (8)即可根据已知的k值直接求解无量纲正常水深x，进而完成正常水深h0计算。

2.2 近似公式的精度分析

式 (8)与式 (7)拟合精度的高低，将直接关系到利用式 (8)计算正常水深的成果精度。为比较近似公式 (8)与由四个分段函数构成的公式 (7)的拟合精度，在工程实用参数范围内，选取不同的xi，并根据xi所在区间代入式 (7)与其相对应的区段函数即可分别计算出与之相对应的ki，再将ki代入式(8)求得与之相对应的，并由式 (9)完成拟合相对误差计算。计算结果见表1。

式中 zi——拟合相对误差，%;

i——在拟合计算中所选取的第i个数据。

表1 式 (8)拟合精度比较

由表1可见，在工程实用参数范围内，用本文式(8)替代由四个区段函数组成的式 (7)的最大拟合替代相对误差仅为0.887%，且位于常用区域边缘，平均相对误差为－0.019%，且拟合相对误差小于0.5%的点占总计算点数的83.3%。可见，本文式(8)具有较高的拟合替代精度，完全可以满足实际工程的设计精度要求。