三次拋物線形渠道斷面收縮水深的簡化計算公式

2012-07-11 11:47王志云

水利技術(shù)監(jiān)督 2012年2期

王志云滕凱

（1.齊齊哈爾市水利工程質(zhì)量與安全監(jiān)督站，黑龍江齊齊哈爾 161006；2.齊齊哈爾市水務局，黑龍江齊齊哈爾 161006）

由于拋物線形斷面具有良好的水流條件及力學性能，因此在城市供排水及水利水電工程中，工程設計部門常根據(jù)輸水流量及地質(zhì)條件，將渠道斷面設計為半立方拋物線、二次拋物線或三次拋物線形式之一，因此相關(guān)拋物線形式的水力計算問題也逐漸引起了有關(guān)學者的重視[1～5]，其中關(guān)于半立方拋物線及二次拋物線形渠道斷面收縮水深的計算已有相關(guān)研究成果[6～9]，而三次拋物線形斷面渠道收縮水深的計算研究則相對較少。

文獻[10]針對三次拋物線形渠道斷面收縮水深計算涉及高次方程的求解問題，棄用常規(guī)的計算過程繁復且成果精度不高的試算法及圖解法，從三次拋物線形斷面渠道收縮水深的基本方程入手，通過對初始迭代值的優(yōu)化擬合，獲得了可直接完成求解的近似計算公式，具有一定的實際意義。但因該公式需分別完成中間參數(shù)及初值計算，求解過程仍顯繁瑣。

為了進一步簡化三次拋物線形斷面收縮水深的計算過程，本文采用優(yōu)化擬合的方法，以標準剩余差最小為目標函數(shù)，獲得了一種表達式較為簡捷、計算精度較高的近似公式。

1 收縮水深的基本計算公式

收縮水深的基本方程為[11]：

式中：E0為以收縮斷面底部為基準面的過水建筑物上游總水頭，m；hc為收縮斷面處的水深，m；Q為過水流量，m3/s；g為重力加速度，m/s2；φ為流速系數(shù)；Ac為三次拋物線形斷面面積，m2。

三次拋物線形斷面的曲線方程為：

其過水斷面面積為：

將式（3）、（4）代人式（1），并設：

經(jīng)整理即可獲得計算三次拋物線形斷面收縮水深的計算公式為：

2 近似計算公式的建立及精度分析

2.1 公式建立

式（6）為高次函數(shù)的超越方程，無法直接獲解。為避免利用式（6）的超越方程求解問題，現(xiàn)通過以下方法尋求其替代函數(shù)：

(1) 根據(jù)式（6）函數(shù)，展繪α～k關(guān)系曲線。

(2) 經(jīng)對α～k關(guān)系曲線的線形分析，初步擬定以下函數(shù)為備選替代函數(shù)。

式中：A、B、C、x、β分別為待定系數(shù)及指數(shù)。

(3) 在工程實用范圍內(nèi)[10]（即 0.01≤α≤0.5），采用優(yōu)化擬合的方法，以標準剩余差最小為目標函數(shù)[12]即：

經(jīng)對式（7）～（10）逐次逼近擬合[13]可知，式（10）的標準剩余差S最小，其替代函數(shù)式為

2.2 精度分析及比較

為比較式（11）與式（6）的擬合精度，考慮在工程實用范圍內(nèi)（即 0.01≤α≤0.5， 2.59≤k≤100.38），取不同的iα值即可由式（6）分別計算出與之相對應的ki，再將ki代入式（11）求得與之相對應的αi′，并由下式完成式（11）替代式（6）的擬合相對誤差，結(jié)果見表1。