基于計算機VB編程解決明渠相關水力計算問題

高 雷

（鄭州大學水利與環境學院，河南鄭州，450001）

基于計算機VB編程解決明渠相關水力計算問題

高 雷

（鄭州大學水利與環境學院，河南鄭州，450001）

在水利工程設計中,經常會遇到一些明渠的水力計算,設計人員普遍采用Excel軟件建立水力計算模塊，借助假設分析中的單變量求解法,實現其計算目標,但是求解操作相當麻煩，計算方法不快捷，本文介紹了利用VB開發可視化計算程序,快速實現明渠不同條件下的水力計算。

VB；明渠；臨界水深；水面曲線

0 前言

水力學分析問題常采用的是一元流理論，然而實際工程問題大多屬于二元流甚至三元流，要精確地分析這些水力學復雜問題，必須從其控制微分方程出發，水力學相關軟件可以簡化這些計算，而縱觀國際上的大型水力學計算軟件,如Delft3D、SMS、MIKE等,其計算程序均使用 FORTRAN 編寫，對設計或施工人員專業知識水平要求較高，而VB編程作為一種解決水力相關計算問題的工具，是借助于計算機數值模擬的方法來實現的，它不僅把水力計算繁瑣或重復的過程程序化，大大節省了工作時間，而且VB編程簡單、面向大眾化，可以根據實際水力條件來設計不同的計算程序。

1 明渠恒定非均勻流臨界水深

1.1 臨界水深計算理論

明渠中由于渠道斷面的幾何形狀或尺寸、粗糙度或底坡沿程改變，或在明渠中修建人工建筑物等都會改變水流的均勻狀態，造成水深和流速等水利要素沿程改變，從而產生非均勻流動，明渠恒定非均勻流特點是明渠的底坡線、水面線、總水頭線彼此互不平行，故水力坡度J、水面坡度JZ、渠底坡度i互不相等，即J≠JZ≠i。在渠道的流量、斷面形狀和尺寸均確定的情況下，相應于斷面比能最小值Esmin，即斷面比能具有極小值=0的水深稱為臨界水深，以hK表示。

1.2 迭代法求解臨界水深

例1:假設某一梯形斷面渠道,設計流量為Q,邊坡系數為m，明渠底寬為b,允許誤差為e，求梯形斷面明渠恒定非均勻流臨界水深hK。

VB編程的程序代碼（迭代法）：

建立迭代法函數:Public Sub diedaifa() A =(b+m

*hk1)* hk1d=b+2 * m * hk1If=1-(q^2)*d/9.8/(A^3) F0=-(q^2)/9.8/(A^3)*(2*m-3*(d^2)/A) hk2=hk1-f/F0

建立循環關系：If Abs(hk2-hk1)<=e Then hk=hk2Else hk1=hk2Call diedaifa End If Print"hk"End Sub

以上即為迭代法計算臨界水深的程序，hK1為水深計算初值，當時，滿足精度要求，即可輸出hK值，通過計算機VB編程減少了人工迭代法計算量，節省了計算時間，并且可以準確輸出計算結果。

2 棱柱體明渠水面曲線

2.1 水面曲線計算理論

當棱柱體明渠中通過一定流量時，由于底坡、上下游進出流邊界條件差異及渠道內建筑物所形成的控制水深不同，明渠中的水流可以形成12種不同形式的水面曲線，通常可以通過經驗對這些水面曲線進行定性判斷，但在實際水利工程中，還需定量知道非均勻流斷面的水力要素變化，如水深、斷面平均流速等，即要對水面曲線進行具體計算和繪制。

2.2 二分法計算棱柱體明渠水面曲線

例2：假設某一棱柱體明渠,設計流量為Q=45m3/s,邊坡系數為m=1.5，明渠底寬為b=10m,粗糙系數n=0.022，底坡i=0.0009，明渠正常水深h0=1.96m，允許誤差為e=0.0001，渠道末端水深h=3.4m，計算該渠道的水面曲線。

VB編程的程序代碼，二分法計算不同段水力要素：

Fori=1 To9 Step1 A0(i)=(b+m*h(i))*h(i)

X(i)=b+2*h(i)*Sqr(1+m*m) J(i)=(v(i)*n/R0(i)^(2 / 3))^2 Es(i)=h(i)+Q^2/(2*g*A0(i)* A0(i))

將VB編程的計算結果輸出，如表一所示，根據表一計算出的各斷面水力要素可以在Excel中畫出想要的棱柱體水面曲線圖（以末端水頭為原點），如圖一所示，水流由右側向左側流動，棱柱體明渠水面曲線為壅水曲線，水深沿程增加，流速逐漸放緩。

表一 各斷面水力要素計算結果

圖一 棱柱體水面曲線

3 總結

水力學計算中經常會遇到一些繁瑣的計算方法，如二分法、迭代法、公式計算法等，這些方法人工計算量大，費時費力，VB開發可視化計算程序恰好巧妙解決了這一問題，把復雜、重復性的工作交給電腦。運用VB編程解決水利工程的實際問題是計算水力學的一部分，隨著計算機VB編程軟件的發展，計算水力學將進一步取代經典水力學中一些近似計算法和圖解法，例如管網和渠系的過水或輸沙（排污）能力的計算、有水輪機負荷改變時水力震蕩系統的穩定性計算研究等等。

[1] 呂宏興,裴國霞，楊玲霞·水力學[M].北京:中國農業出版社,2011.6[2] 馮民權，趙明登，鄭邦民·河渠非恒定流及其物質輸運的數值模擬[M].北京:科學出版社，2012.2

[3] 李占松，王玲玲，朱士江.二維淺水流動數學模型的多解性評價[J].河 南 科 學，2009,27（1）:64-66

[4] 吳丹鳳，丁全林，戴會超，戴凌全.計算水力學通用后處理模塊設計與快速開發[J].水電能源科學，2013,31（6）:195-197

高雷（1991.07- ），男，山東濰坊人，鄭州大學水利與環境學院 2011級本科生，水利水電工程專業

Based on computer VB programming to solve the hydraulic calculation problems related to water conservancy project

Gao Lei
（College of Water Conservancy & Environmental Engineering，Zhengzhou University，Zhengzhou，450001，Henan，China）

In the design of water conservancy projects,often encounter some hydraulic calculation of flow in open channel,designers generally Excel software is adopted to establish the hydraulic calculation module, with the help of assumption in the analysis of single variable method,achieving the calculation,but it's quite a trouble,solve the operation calculation method is not quick,visual calculation program,this paper introduces the use of VB rapidly open channel under the condition of different hydraulic calculation.

VB；Open channel；The critical depth of water；The surface curve