張志昌，賈 斌，李若冰

(西安理工大學(xué) 水利水電學(xué)院，陜西 西安710048)

六圓弧蛋形斷面正常水深直接計算方法

張志昌，賈 斌，李若冰

(西安理工大學(xué) 水利水電學(xué)院，陜西 西安710048)

【目的】 研究六圓弧蛋形斷面正常水深的簡單計算方法，為該斷面的工程設(shè)計提供技術(shù)支持。【方法】 根據(jù)面積分割法計算六圓弧蛋形斷面的面積、濕周和水力半徑，根據(jù)明渠均勻流理論計算不同流量時渠道的正常水深，通過優(yōu)化擬合研究蛋形斷面正常水深的簡單計算方法。【結(jié)果】 優(yōu)化擬合得到了蛋形斷面水深處于不同位置時相對斷面面積、相對濕周、相對水力半徑的計算公式，確定了正常水深與流量的關(guān)系及水深的迭代計算公式，通過優(yōu)化擬合研究，最后給出了正常水深的簡化計算公式?！窘Y(jié)論】 所建立的六圓弧蛋形斷面正常水深的理論計算公式和近似計算公式計算過程簡單、精度高，平均誤差為0.204%，最大誤差為1.129%。

六圓弧蛋形斷面；正常水深；計算方法；參數(shù)修改

正常水深的計算在明渠設(shè)計、運行、維護中具有重要作用。明渠斷面有矩形、梯形、U形、馬蹄型、蛋形等多種形式，對于矩形和梯形斷面正常水深的研究已有很多成果[1-3]，對于U形斷面文獻[4]給出了直接計算公式，馬蹄型斷面分為標準Ⅰ型和標準Ⅱ型，對標準Ⅰ型馬蹄型斷面文獻[5]給出了正常水深的迭代公式，文獻[6]給出了標準Ⅱ型馬蹄型斷面正常水深的直接計算公式。

蛋形斷面由于受力條件好，水流條件優(yōu)越，常用于城市的排水工程和水工建筑物的隧洞工程。例如美國曼哈頓的排水工程、印度卡納塔克邦Varahi水電站的輸水工程[7]，以及我國湖南省臨澧縣青山水輪泵站灌區(qū)的冉鋪灣隧洞[8]和衡東縣白蓮灌區(qū)武家坳隧洞[9]均采用蛋形斷面。

蛋形斷面形狀比較復(fù)雜，形式多樣。斷面形式主要有下部為半圓形、上部為半橢圓形的蛋形斷面[10]；底部為半圓形、上部為拋物線形的蛋形斷面[9]和四圓弧蛋形斷面[11]及六圓弧蛋形斷面[7-8]等。目前有關(guān)蛋形斷面水力特性的研究主要集中在正常水深和臨界水深方面。就正常水深而言，文獻[10]研究了下部為半圓、上部為半橢圓的蛋形斷面正常水深的近似計算方法；文獻[12]運用CAD技術(shù), 編制了倒立式四圓弧蛋形斷面諾謨圖；文獻[13]也對倒立式四圓弧蛋形斷面的水力要素進行了分析，給出了正常水深的計算方法；文獻[14]將倒立式四圓弧蛋形斷面分成3個區(qū)域，分別給出了正常水深的迭代公式，與理論公式相比該迭代公式的誤差小于1%；文獻[15]研究了四圓弧蛋形斷面隧洞正常水深的簡易算法，但從公式結(jié)構(gòu)來看其計算過程比較復(fù)雜；文獻[11]研究了四圓弧蛋形斷面的正常水深，但在計算時需借助表格查算，計算過程比較麻煩；文獻[16]依據(jù)優(yōu)化擬合理論，給出了我國目前應(yīng)用的六圓弧蛋形斷面正常水深的簡化計算方法，但在計算濕周時誤將底部弓形斷面的半角當作全角，公式存在著明顯的錯誤；文獻[7]研究了另一種形式的六圓弧蛋形斷面正常水深和臨界水深的計算方法，與我國目前應(yīng)用的六圓弧蛋形斷面相比，其圓弧的圓心位置不同。

綜上可以看出，蛋形斷面正常水深的計算主要是針對四圓弧蛋形斷面或者圓弧和橢圓組合形式。六圓弧蛋形斷面由6個圓心和3個不同半徑的圓弧組成，幾何形狀復(fù)雜，計算過程繁瑣。因此，本研究在詳細推導(dǎo)六圓弧蛋形斷面正常水深計算公式的基礎(chǔ)上，分析了蛋形斷面相對面積、相對水力半徑、相對濕周、相對流量與相對水深的關(guān)系，以期確定該斷面正常水深的簡化計算方法，從而為六圓弧蛋形斷面的工程設(shè)計提供技術(shù)支撐。

1 蛋形斷面的形狀及參數(shù)

六圓弧蛋形斷面是在文獻[17]的基礎(chǔ)上改進的。文獻[17]在設(shè)計中取r1=280 cm，r2=120 cm，r3=85 cm，α1=25.375°(0.442 88 rad)，α2=122.2°(2.132 79 rad)，α3=27.066 7°(0.472 4 rad)，α4沒有給出角度，渠道最大寬度為r1，最寬處距渠底的距離為82 cm，但使用計算機繪圖時發(fā)現(xiàn)，在每個圓弧與圓弧相切點都找不到切點。因此作者對原設(shè)計進行了修改，在修改中仍保持r1、r3、α1，由圖1的幾何關(guān)系計算得：hab=27.013 74 cm，h2=54.782 39 cm，h3=78.912 02 cm，r2=119.291 86 cm，α2=121.183 8°(2.115 06 rad)，α3=29.408 1°(0.513 27 rad)，α4=40.128 11°(0.700 37 rad)。為了使蛋形斷面具有通用性，文獻[18]將hab、h2、h3、r2、r3均換成r1的函數(shù)，得hab=0.096 478r1，h2=0.195 651 4r1，h3=0.281 829r1，r2=0.426 042r1，r3=(17/56)r1。

現(xiàn)以圖1中ab、cd、ef作為分界線，對應(yīng)水深分別為hab、hc d、hef，則可以將斷面水深h分為4種情況，即h≤hab、hab≤h≤hc d、hc d≤h≤hef和h≥hef。

圖1 蛋形斷面設(shè)計示意圖Fig.1 Diagrammatic design of egg-shaped cross section

2 不同水深時蛋形斷面面積、濕周、水力半徑和水深的計算

2.1 水深處于底部弓形斷面(h≤hab)

如圖2所示，當斷面實際水深h處于底部弓形，即h≤hab時，斷面相對面積、相對濕周、相對水力半徑及相對水深按照下列公式計算。

斷面相對面積為：

(1)

式中：A0為水深h對應(yīng)的斷面面積；φ為水深h處于底部弓形斷面時對應(yīng)的半角。

相對濕周為：

χ/r1=2φ。

(2)

式中：χ為濕周。

相對水力半徑為：

R/r1=A0/(χr1)=(2φ-sin 2φ)/(4φ)。

(3)

式中：R為水力半徑。

相對水深為：

h/r1=1-cosφ。

(4)

以上各式中，0<φ≤25.375°(0.442 88 rad)。

2.2 水深處于最大寬度以下(含最大寬度)斷面(hab≤h≤hcd)

如圖3所示，當斷面實際水深h處于最大寬度以下時，可將計算區(qū)域分為底部、兩側(cè)弓形和中央梯形4個部分，其面積分別為A1、A2、A3、A4，文獻[18]給出的相對斷面面積計算公式為：

(5)

式中：α4=40.128 11°(0.700 37 rad)，sin (α4-β)=sinα4cosβ-cosα4sinβ，代入公式(5)有：

(6)

同理，相對濕周、相對水力半徑和相對水深的計算公式為：

χ/r1=2α1+(34/56)(α4-β)=1.310 978-(34/56)β，

(7)

R/r1=[0.242 5-(17/56)2(β+sinβcosβ)-0.119 26sinβ]/[1.310 978-(34/56)β]，

(8)

h/r1=(hab+h2-r3sinβ)/r1=0.292 129-(17/56)sinβ。

(9)

式中：0≤β≤40.128 11°(0.700 37 rad)。

圖2 水深處于底部弓形斷面Fig.2 Water depth at the bottom bow of cross section

2.3 水深處于ef線(含ef線)以下斷面(hcd≤h≤hef)

(10)

χ/r1=2α1+(34/56)α4+2β1=1.310 978+2β1，

(11)

R/r1=(0.242 5+β1-sinβ1+sinβ1cosβ1)/(1.310 978+2β1)，

(12)

h/r1=0.292 129+sinβ1。

(13)

式中：0≤β1≤29.408 1°(0.513 27 rad)。

2.4 水深處于頂部弓形(ef線以上)斷面(h≥hef)

當斷面水深h處于圖5的ef線以下時，在ef線以下(含ef線)的相對斷面面積可以由公式(10)計算，此時的β1=α3=29.408 1°(0.513 27rad)。當斷面水深h處于圖5的ef線以上時，在ef線以上面積可以分為兩側(cè)弓形面積A8、A9及中間梯形面積A10。

相對斷面面積、相對濕周、相對水力半徑和相對水深的計算公式為：

(14)

χ/r1=2.337 514+0.852 084β2，

(15)

R/r1={0.614 854+0.181 512[β2+0.5sin (α-2β2)]}/(2.337 514+0.852 084β2)，

(16)

h/r1=0.573 958+0.426 042cos (α2/2-β2)。

(17)

式中：0≤β2<α2/2=60.591 9°(1.057 53 rad)。

圖4 水深處于ef線(含ef線)以下斷面 Fig.4 Water depth under the ef line of cross section

3 蛋形斷面正常水深的計算

3.1 不同斷面水深條件下正常水深的理論計算

蛋形斷面正常水深的計算公式為明渠均勻流公式[19]，即：

(18)

式中：Q為流量，A為過水斷面面積，C為謝才系數(shù)，R為水力半徑，i為明渠的比降，n為糙率，χ為濕周。

由前面研究可以看出，六圓弧蛋形斷面可以分為h≤hab、hab≤h≤hcd、hcd≤h≤hef和h≥hef4種斷面水深條件，表1給出了4種水力條件下已知水深求流量和已知流量求水深的計算公式及迭代初值。

迭代公式經(jīng)驗算均收斂，由于證明過程較繁，此處從略。

表1 六圓弧蛋形斷面正常水深的計算公式Table 1 Formula for normal water depth of Six arc egg-shaped section

3.2 蛋形斷面正常水深的近似計算

圖6 當0

當0

(19)

(19)式的平均誤差為0.336%，最大誤差為1.074%。

當0.15≤h/r1<0.90時，有：

(20)

(20)式的平均誤差為0.204%，最大誤差為1.129%。

(21)

(22)

在迭代時，h/r1初值可取0.15～0.90。

公式(22)需要迭代求解，為了簡化計算，給出的顯函數(shù)關(guān)系式為：

(23)

公式(23)的平均誤差為0.229%，最大誤差為1.057%。

4 算 例

已知某輸水隧洞為六圓弧蛋形斷面，最大半徑r1=1.5 m，洞底設(shè)計坡降i=1/2 500，洞內(nèi)壁糙率n=0.014，試求當設(shè)計流量Q=0.1～1.2 m3/s時隧洞內(nèi)相應(yīng)的正常水深h值。

根據(jù)設(shè)計流量及已知條件，由表1迭代法及公式(22)、(23)得到的正常水深h的計算結(jié)果及誤差分析見表2。

表2 某六圓弧蛋形斷面輸水隧洞不同流量下正常水深的計算結(jié)果及誤差分析Table 2 Calculation results and error analysis of normal water depth of egg-shaped cross section with six arcs in water-conveyance channel

由表2可以看出，h理論計算值與用簡化公式計算的h值非常接近，最大誤差為0.723%。而文獻[16]由于濕周公式有誤，在流量Q=0.1和1.2 m3/s時，計算得到的h=0.19和0.934 3 m，與h理論計算結(jié)果分別相差18.78%和13.77%。由此可見，水深越小，濕周對流量計算的影響越大。

5 結(jié) 語

本研究根據(jù)面積分割法和明渠均勻流理論，推導(dǎo)了六圓弧蛋形相對斷面面積、相對濕周和相對水力半徑的計算公式，給出了水深處于不同位置時正常水深與流量的關(guān)系，在此基礎(chǔ)上，擬合了相對水深與相對流量的簡單計算公式，經(jīng)驗證，該公式不僅計算方便而且精度較高，完全可以滿足工程設(shè)計要求。

[1] 張志昌，肖宏武，毛兆民.明渠測流的理論和方法 [M].西安：陜西人民出版社，2004：292-293.

Zhang Z C,Xiao H W,Mao Z M.The theory and method of the open channel flow [M].Xi’an：Shaanxi People’s Publishing House,2004：292-293.(in Chinese)

[2] 趙延風，祝晗英，王正中，等.梯形明渠正常水深的直接計算方法 [J].西北農(nóng)林科技大學(xué)學(xué)報：自然科學(xué)版,2009,37(4)：220-223.

Zhao Y F,Zhu H Y,Wang Z Z,et al.A direct formula for normal depth in trapezoidal open channel [J].Journal of Northwest A&F University：Natural Science Edition,2009,37(4)：220-223.(in Chinese)

[3] 王正中，席跟戰(zhàn)，宋松柏，等.梯形明渠正常水深直接計算公式 [J].長江科學(xué)院院報,1998,15(6)：1-3.

Wang Z Z,Xi G Z,Song S B,et al.A direct calculation formula for normal depth in open trapezoidal channel [J].Journal of Yangtze River Scientific Research Institute,1998,15(6)：1-3.(in Chinese)

[4] 張新燕，呂宏興，朱德蘭.U形渠道正常水深的直接水力計算公式 [J].農(nóng)業(yè)工程學(xué)報，2013，29(14)：115-119.

Zhang X Y,Lü H X,Zhu D L.Direct calculation formula for normal depth of U-shaped channel [J].Transactions of the CSAE,2013,29(14)：115-119.(in Chinese)

[5] 張志昌，李若冰.標準Ⅰ型馬蹄形斷面水力特性的研究 [J].長江科學(xué)院院報，2013,30(5)：55-59.

Zhang Z C，Li R B.Hydraulic characteristics of standard Ⅰ-type horseshoe cross section [J].Journal of Yangtze River Scientific Research Institute，2013,30(5)：55-59.(in Chinese)

[6] 趙延風，王正中，蘆 琴.馬蹄形斷面正常水深的直接計算公式 [J].水力發(fā)電學(xué)報，2012,31(1)：173-177.

Zhao Y F，Wang Z Z，Lu Q．Direct calculation formulas for normal depth of horseshoe section [J].Journal of Hydroelectric Engineering，2012，31(1)：173-177．(in Chinese)

[7] Raikar R V,Shiva Reddy M S,Vishwanadh G K.Normal and critical depth computations for egg-shaped conduit sections [J].Flow Measurement and Instrumentation,2010,21(3):367-372.

[8] 傅功年，唐耿紅.冉鋪灣隧洞襯砌斷面型式選擇 [J].人民長江，2002，33(5)：14-16.

Fu G N，Tang G H.Type selection of lined cross-section for Ranpuwan tunnel [J].Yangtze River,2002,33(5)：14-16.(in Chinese)

[9] 王子宜.蛋形襯砌在白蓮灌區(qū)武家坳隧洞中的應(yīng)用 [J].湖南水利水電，2006(6)：70-71.

Wang Z Y.Application of egg-shaped lining in Bailian irrigation area Wujiaao tunnel [J].Hunan Hydro & Power,2006(6):70-71.(in Chinese)

[10] 李風玲,文 輝,涂寧宇.蛋形斷面管道正常水深近似算法 [J].人民長江，2008，39(18)：77-78.

Li F L，Wen H，Tu N Y.Approximation algorithm of normal depth in egg-shaped cross section pipeline [J].Yangtze River，2008，39(18):77-78.(in Chinese)

[11] Gill M A.Hydraulics of partially filled egg sewers [J].Journal of Environmental Engineering,1987,113(2)：407-425.

[12] 寧希南.蛋形斷面管道簡明水力計算及其諾謨圖[J].林業(yè)建設(shè)，2007(5):22-25.

Ning X N.Concise hydraulic calculation and nomogram of egg-shaped cross section pipeline [J].Forestry Construction，2007(5):22-25.(in Chinese)

[13] 陶 冶，賈 悅.蛋形斷面排水無壓管道水力計算及其水力要素分析 [J].給水排水，2009，35(S1)：428-430.

Tao Y,Jia Y.Analysis of the hydraulic calculation and hydraulic elements of egg-shaped cross section in open channel [J].Water & Wastewater Engineering,2009，35(S1)：428-430.(in Chinese)

[14] Bijankhan M,Kouchakzadeh S.Egg-shaped cross section:Uniform flow direct solution and stability identification [J].Flow Measurement and Instrumentation，2011,22(6):511-516.

[15] 滕 凱.蛋形斷面隧洞正常水深的簡易算法 [J].長江科學(xué)院院報，2013,30(12):39-42.

Teng K.The simple algorithm of normal depth in egg-shaped cross section tunnel [J].Journal of Yangtze River Scitentific Research Institute，2013,30(12):39-42.(in Chinese)

[16] 滕 凱，李新宇.六圓弧蛋形斷面無壓隧洞水面線解析計算模型 [J].水資源與水工程學(xué)報，2013,24(4)：174-179.

Teng K,Li X Y.Calculation model of water surface line analysis in tunnel of six circular egg-shaped section without pressure [J].Journal of Water Resources&Water Engineering,2013,24(4):174-179.(in Chinese)

[17] 武漢水利電力學(xué)院.水工建筑物 [M].北京：水利出版社，1980：188-190.

Wuhan Institute of Hydraulic and Electric Engineering.Hydraulic structure [M].Beijing:Water Conservancy Press,1980:188-190.(in Chinese)

[18] 李若冰，張志昌.明渠六圓弧蛋形斷面臨界水深和收縮斷面水深的計算 [J].武漢大學(xué)學(xué)報，2012，45(4)：463-467.

Li R B,Zhang Z C.Calculation of critical water depth and contraction section water depth in egg-shaped cross-section with six arcs of open channel [J].Engineering Journal of Wuhan University,2012,45(4):463-467.(in Chinese)

[19] 張志昌.水力學(xué) [M].北京：中國水利水電出版社，2011.

Zhang Z C.Hydraulics [M].Beijing:China Water Power Press,2011.(in Chinese)

Direction calculation of normal water depth of egg-shaped cross section with six arcs in open channel

ZHANG Zhi-chang，JIA Bin，LI Ruo-bing

(InstituteofWaterResourcesandHydro-electricEngineering,Xi’anUniversityofTechnology，Xi’an，Shaanxi710048，China)

【Objective】 This study proposed a simple method to calculate the normal water depth of egg-shaped cross section with six arcs to provide basis for engineering design of such section.【Method】 Based on area segmentation method,the area,wetted perimeter and hydraulic radius were computed.Using open channel uniform flow theory,the normal depths were studied under different flows and the approximate fitting was used to obtain simple calculation method of normal water depth.【Result】 Approximate fitting was used to obtain the simple calculation methods of section area, wetted perimeter and hydraulic radius at different positions.The relationship between normal depth and rate of flow,the iterative calculation formula of the water depth and the relationship between relative water depth and relative flow were also determined.At last,the simplified calculation method of normal water depth was proposed by optimizing fitting.【Conclusion】 The proposed theoretical calculation formula and approximate calculation formula of normal water depth of egg-shaped cross section were simple and precise with average error of 0.204% and maximum error of 1.129%.

egg-shaped cross section with six arcs;normal depth;computing method;parametric modification

2014-01-20

張志昌(1954-)，男，陜西西安人，教授級高級工程師，主要從事水工水力學(xué)研究。 E-mail：zhangzhichang1954@163.com

時間:2015-06-30 13:47

10.13207/j.cnki.jnwafu.2015.08.022

TV133

A

1671-9387(2015)08-0229-06

網(wǎng)絡(luò)出版地址:http://www.cnki.net/kcms/detail/61.1390.S.20150630.1347.022.html