• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    二分之五次方拋物線形明渠設(shè)計(jì)及提高水力特性效果

    2017-10-13 16:40:52韓延成徐征和高學(xué)平SaidEasa
    關(guān)鍵詞:次方過流水深

    韓延成,徐征和,高學(xué)平,Said M. Easa

    ?

    二分之五次方拋物線形明渠設(shè)計(jì)及提高水力特性效果

    韓延成1,徐征和1,高學(xué)平2,Said M. Easa3

    (1. 濟(jì)南大學(xué)資源與環(huán)境學(xué)院,濟(jì)南 250022; 2. 天津大學(xué)水利工程仿真與安全國家重點(diǎn)實(shí)驗(yàn)室,天津 300072;3. Dept. of Civil Engineering, Ryerson Univ., Toronto, ON, Canada, M5B 2K3)

    為提高拋物線形斷面的水力特性,增加輸水能力,該文提出了一種二分之五次(以下簡(jiǎn)稱2.5次)方拋物線形渠道斷面,推導(dǎo)其水力斷面特性。將濕周用高斯超幾何函數(shù)表示后,將水力最優(yōu)斷面的最優(yōu)化模型轉(zhuǎn)換為關(guān)于寬深比的一元方程,得到2.5次方拋物線形渠道水力最優(yōu)斷面的解析解,其最優(yōu)寬深比為2.088 3。比較結(jié)果表明,2.5次方拋物線形斷面較常規(guī)拋物線形斷面具有更好的水力學(xué)特性。與平方、半立方拋物線形斷面比較,在相同水深條件下,2.5次方拋物線形水力最優(yōu)斷面的過流能力更大。相反,在相同流量下,2.5次方拋物線形水力最優(yōu)斷面的過流面積、濕周、水深更小。2.5次方拋物線形水力最優(yōu)斷面的建造成本與其他2種斷面相比是最小的。進(jìn)一步地,為便于工程應(yīng)用,基于高斯勒讓德算法,提出2.5次方拋物線形斷面的三點(diǎn)和四點(diǎn)格式近似濕周算法。結(jié)果表明,四點(diǎn)格式近似算法具有較高精度。研究可為明渠設(shè)計(jì)提供理論依據(jù)。

    渠道;水力學(xué);設(shè)計(jì);明渠;二分之五次方拋物線形斷面;水力最優(yōu)斷面;濕周

    0 引 言

    渠道輸水?dāng)嗝鎸?duì)渠道過流能力、水深、建造成本等均有很大的影響。常見的渠道輸水?dāng)嗝鏋樘菪?、矩形斷面。隨著渠道建造工藝的改進(jìn),大型襯砌機(jī)器的應(yīng)用,曲線形斷面的建造越來越容易,也越來越受到歡迎。例如巴基斯坦的High Level渠,西班牙Genil-Cabra渠等[1]均采用了平方拋物線形斷面。學(xué)者們普遍認(rèn)為曲線形斷面有以下優(yōu)點(diǎn)[2-4]:曲線形渠道斷面沒有或拐角點(diǎn)少、應(yīng)力集中點(diǎn)少,因而由應(yīng)力集中可能產(chǎn)生的裂縫少,滲漏量少;自然形成的河道或非襯砌渠道更多呈現(xiàn)曲線形斷面的形狀;曲線形渠道斷面從渠底到渠堤邊坡是漸變的,因而穩(wěn)定性更好;曲線形渠道具有更好的水力學(xué)特性。

    常見的拋物線形渠道斷面主要包括平方拋物線形斷面和半立方拋物線形斷面。學(xué)者們對(duì)常規(guī)的拋物線形斷面已經(jīng)進(jìn)行了大量研究,不僅研究了平方拋物線形渠道的水力最優(yōu)斷面[2,5]、經(jīng)濟(jì)斷面[6],也研究了具有拋物線形底和三角形邊坡的復(fù)合斷面[7],平底的平方拋物線形斷面[3-4],以及平底的半立方拋物線形斷面[8]等。為方便工程應(yīng)用,張新燕等[9]進(jìn)行了拋物線形斷面渠道正常水深的顯式計(jì)算研究,文輝等[10]進(jìn)行了平底的馬蹄形斷面的水力計(jì)算,王正中等[11]進(jìn)行了收縮斷面的顯式計(jì)算研究,文輝等[12]進(jìn)一步開展了拋物線形渠道恒定漸變流水面線的計(jì)算方法的研究。然而,在寬度和深度相同的情況下,平方拋物線形和半立方拋物線形斷面的底部較尖,可能仍然會(huì)影響過流能力。為提高拋物線形斷面的水力學(xué)特性,提高輸水能力,本文提出了一種二分之五次(以下簡(jiǎn)稱2.5次)方拋物線形斷面,并對(duì)該拋物線形的水力斷面特性、最優(yōu)斷面的解析解、近似濕周算法等進(jìn)行研究,以便于工程應(yīng)用。

    1 2.5次方拋物線形渠道斷面及其特性

    1.1 斷面形狀

    2.5次拋物線形斷面形狀如圖1所示。與平方拋物線和半立方拋物線相比較,2.5次拋物線形斷面具有更加平坦的底部,這有利于施工、底部壓實(shí)及增加過流能力。

    以拋物線斷面底部中心為原點(diǎn)建立-坐標(biāo)系,2.5次拋物線形斷面可表示為

    式中為形狀系數(shù)。

    由圖1可知,當(dāng)為水面寬度的1/2(=/2)時(shí),存在=(為水深),可用和表示為

    由式(2),可用和表示為

    1.2 面積與濕周

    2.5次方拋物線渠道過水?dāng)嗝娴拿娣e可用積分方法得到,為

    2.5次方拋物線渠道的濕周可用弧長法求得,為

    式中為濕周,m。

    1.3 濕周的顯式近似解

    式(5)中的濕周為積分形式,需要用數(shù)值方法求得。為便于工程應(yīng)用,求解簡(jiǎn)單的顯式求解算法是必要的。本文經(jīng)過推導(dǎo),將式(5)用高斯勒讓德三點(diǎn)格式表示為[13]

    利用式(2),也可用和表示為

    用同樣的方法,得到基于高斯勒讓德四點(diǎn)格式的濕周表達(dá)式為

    以為0.3為例,比較三點(diǎn)法和四點(diǎn)法計(jì)算2.5次方拋物線形渠道為1.0~3.5 m時(shí)濕周的精度。將=0.3和值代入式(3)得到,然后代入式(5),用數(shù)值積分法求得濕周理論值P。同樣,將和分別代入式(6)和式(8)得到基于三點(diǎn)和四點(diǎn)高斯勒讓德的濕周近似值3和4(表1)。表中絕對(duì)誤差為近似值與理論值差值的絕對(duì)值。由表1可知,三點(diǎn)法計(jì)算的濕周的最大絕對(duì)誤差為0.004 01 m,四點(diǎn)法濕周的最大絕對(duì)誤差0.000 97 m??梢钥闯?,四點(diǎn)高斯勒讓德近似算法具有較高的精度,更接近理論值。

    表1 基于高斯勒讓德近似算法的濕周值

    2 2.5次方拋物線形水力最優(yōu)斷面的解析解

    2.1 明渠均勻流

    渠道設(shè)計(jì)流量可按謝才公式表示為[10,13-15]

    式中為設(shè)計(jì)流量,m3/s;為渠底縱坡;為水力半徑,m;C為謝才系數(shù)。根據(jù)曼寧公式[13,16-19],(為糙率)。因此式(9)可以表示為

    2.2 水力最優(yōu)斷面求解模型

    水力最優(yōu)斷面是面積一定的情況下,使過流能力最大的斷面,或過流能力一定的情況下,使過流面積最小的斷面[17, 20-23]。兩者得到的最終結(jié)果是相同的。因此,求解水力最優(yōu)斷面的模型表示為

    目標(biāo)函數(shù):

    約束條件:

    式中為等式約束函數(shù)。

    上式用拉格朗日乘子法可表示為

    式中為拉格朗日乘子。

    將式(14)代入式(13),消去,并化簡(jiǎn)后可得到

    式(15)即求解2.5次方拋物線形渠道最優(yōu)水力斷面的微分方程。

    2.3 水力最優(yōu)斷面的解析解

    由式(4)得到對(duì)和的偏導(dǎo)數(shù)

    式(5)中濕周是一個(gè)積分形式,難以求導(dǎo)。式(6)或式(8)是近似算法,不適合求解水力最優(yōu)斷面的理論解。此處可以用高斯超幾何函數(shù)表示為[24-25]

    其中S=Hypergeom([1,2],3,4) ,為高斯超幾何函數(shù),14為高斯超幾何函數(shù)的參數(shù)。

    由式(10),可以求得對(duì)和的偏導(dǎo)數(shù)分別為

    將式(16)、(17)、(20)、(21)代入式(15),簡(jiǎn)化后得到

    (22)

    設(shè)無量綱參數(shù)=/,并將其代入式(22),經(jīng)過進(jìn)一步簡(jiǎn)化可以得到

    式(23)只有一個(gè)變量,用數(shù)值方法求解式(23),可以得到唯一的可行解

    因此,2.5次方拋物線形最優(yōu)斷面的寬深比(水面寬與水深比)是一個(gè)常數(shù),/=2.088 3。由可以得到其他參數(shù),將代入式(2)可以得到最優(yōu)斷面的形狀系數(shù)

    將/=2.088 3和式(25)代入式(4)可以得到水力最優(yōu)斷面的過流面積的直接計(jì)算公式為

    =1.491 62(26)

    將/=2.088 3和式(25)代入式(18)可以得到水力最優(yōu)斷面的濕周的直接計(jì)算公式為

    =3.096 33(27)

    2.4 水深、水面寬度、過流面積、濕周與流量的關(guān)系

    將式(26)、(27)代入式(10),可以得到流量的顯式計(jì)算公式

    對(duì)上述方程求解,可以得到依據(jù)流量直接計(jì)算的顯式公式

    將式(29)代入式(24)~式(27),可以得到

    2.5 正常水深和臨界水深的計(jì)算

    正常水深和臨界水深是明渠水力學(xué)計(jì)算的重要參數(shù)[26-29]。顯然式(29)也是最優(yōu)斷面條件下正常水深的計(jì)算公式。

    臨界水深的通用計(jì)算公式為

    式中為能量修正系數(shù);為重力加速度,m/s2。

    將/=2.088 3和式(26)代入式(35),可以得到臨界水深h的顯式表達(dá)式

    2.6 案例分析

    2.6.1案例描述

    一個(gè)2.5次方拋物線形渠道,流量=25 m3/s,=0.014,=1/12 000,=1。要求:1)按2.5次方最優(yōu)拋物線形斷面設(shè)計(jì)渠道。2)驗(yàn)證最優(yōu)斷面寬深比為2.088 3。

    2.6.2結(jié)果與分析

    1)最優(yōu)斷面設(shè)計(jì)

    將已知條件代入式(29)可得到水深4.055 5 m,代入式(25)得最優(yōu)形狀系數(shù)0.109 91。由/=2.088 3得到=8.469 1 m。代入式(36)得臨界水深h=2.092 2 m。代入式(26)~(27)或式(32)~(33)得24.533 m2、=12.557 m。

    2)最優(yōu)斷面寬深比驗(yàn)證

    取=24.533 m2,取0.5~10 m內(nèi)步長0.000 5 m的不同值,根據(jù)式(4)計(jì)算得到,將和代入式(2)得到。利用式(5)得到,式(10)得到,繪制/與的關(guān)系曲線。如圖2所示,面積一定時(shí),=2.09時(shí)最大,其結(jié)果與式(24)(2.088 3)近似。

    3 2.5次方和常規(guī)拋物線形最優(yōu)斷面的比較

    3.1 過流能力的比較

    與平方拋物線形的比較。對(duì)于平方拋物線形斷面的水力最優(yōu)斷面,、、和分別為

    2.055 5,=1.370 32,=2.998,=0.946 7-1(37)

    將式(37)代入式(10),可以得到平方拋物線形斷面流量的顯式計(jì)算公式為

    求解上式,可以得到平方拋物線形斷面正常水深關(guān)于的表達(dá)式為

    比較式(38)和式(28)可知,相同水深情況下,2.5次方拋物線形最優(yōu)水力斷面的流量大于平方拋物線形最優(yōu)斷面的流量;比較式(39)與式(34)也可以得到,相同流量情況下,2.5次方拋物線形最優(yōu)水力斷面的水深小于平方拋物線形最優(yōu)斷面的水深。

    將式(39)的代入式(37),可以得到平方拋物線形斷面的,和的顯式計(jì)算公式

    用同樣的方法,可以得到半立方拋物線形的特性參數(shù)如表2。2.5次方拋物線形最優(yōu)水力斷面的過流面積、濕周較半立方拋物線形、平方拋物線形斷面的小。

    3.2 建設(shè)成本比較

    渠道的主要建設(shè)成本包括土方、襯砌和征地費(fèi)用。忽略超高的影響后,通用的單位渠長上建造成本可表示為

    =W·+W·+W·(43)

    式中為單位渠長上建造成本,元;W為單位渠長上,沿橫斷面上單位面積挖土方成本,元/m2;W為單位渠長上,沿橫斷面單位長度襯砌成本,元/m;W為單位渠長上,單位渠寬的征地費(fèi),元/m。

    表2 不同拋物線形渠道水力最優(yōu)斷面的水力特性

    注:,為流量,為糙度,為渠底縱坡。

    Note:,is discharge,is roughness,is slope of canal base。

    對(duì)2.5次方拋物線形斷面,如果不考慮安全超高的影響,將式(31)~式(33)代入式(43),可以得到2.5次方拋物線形斷面單位長度渠道造價(jià)為

    式中2.5為2.5次方拋物線形最優(yōu)斷面單位渠長建造成本,元;。

    同樣,將式(40)~式(42)代入式(43),可以得到單位渠長平方拋物線形最優(yōu)斷面的建造成本

    式中2為平方拋物線形斷面單位渠長建造成本,元。

    根據(jù)表2和式(43),也可以得到半立方拋物線形斷面建造成本(單位渠長)

    式中1.5為半立方拋物線形斷面單位渠長建造成本,元。

    比較式(44)~式(46),可以得到

    因此,2.5次方拋物線形最優(yōu)斷面的造價(jià)是最低的。

    實(shí)際上,考慮安全超高后,用類似的方法,可得到2.5次方拋物線形最優(yōu)斷面的建造成本小于平方拋物線和半立方拋物線形斷面。推導(dǎo)過程不再贅述。

    4 結(jié)論與討論

    斷面設(shè)計(jì)是輸水渠道設(shè)計(jì)的最重要內(nèi)容之一,其對(duì)渠道輸水效率、成本均有顯著影響。本文提出了一種二分之五次(以下簡(jiǎn)稱2.5次)方拋物線形明渠輸水?dāng)嗝?,并?duì)其水力特性、水力最優(yōu)斷面等進(jìn)行了研究。

    2.5次方拋物線形明渠輸水?dāng)嗝婢哂休^平坦的渠底,具有良好的水力特性。推導(dǎo)了2.5次方拋物線形渠道水力斷面特性。為便于工程應(yīng)用,提出了基于高斯勒讓德算法的三點(diǎn)和四點(diǎn)近似濕周算法,結(jié)果表明,四點(diǎn)近似算法具有很高的精度。

    建立了2.5次方拋物線形斷面的水力最優(yōu)斷面模型,用拉格朗日法推導(dǎo)了最優(yōu)斷面的微分方程。將濕周用高斯超幾何函數(shù)表示后,將水力最優(yōu)斷面模型轉(zhuǎn)化成了關(guān)于寬深比的方程,并最終求得了2.5次方拋物線形斷面的水力最優(yōu)斷面。結(jié)果表明其最佳寬深比為2.088 3。

    利用水力最佳寬深比,進(jìn)一步得到了水力最優(yōu)斷面條件下的正常水深、臨界水深、水面寬、過流面積、濕周的顯式計(jì)算公式。對(duì)比結(jié)果表明,水深相同的情況下,2.5次方拋物線形最優(yōu)水力斷面的流量較半立方拋物線形、常規(guī)(平方)拋物線形斷面的大;相反,在流量相同的情況下,2.5次方拋物線形最優(yōu)水力斷面的過流面積、濕周較半立方拋物線形、平方拋物線形斷面的小。建造成本對(duì)比結(jié)果表明,相同情況下,2.5次方拋物線形最優(yōu)水力斷面的建造成本較半立方拋物線形、平方拋物線形斷面的小。

    到目前為止,包括梯形、矩形、拋物線形、半圓形、蛋形等在內(nèi)的多種斷面已被應(yīng)用于渠道或排水管道。這使設(shè)計(jì)者可以根據(jù)地質(zhì)條件、輸水規(guī)模、滲漏、經(jīng)濟(jì)條件等選擇最合適的斷面形式。本文提出的2.5次方拋物線形斷面為工程提供了一種新的選擇。但在實(shí)際工程中,仍然需要根據(jù)過流能力、水深、地質(zhì)條件、經(jīng)濟(jì)性等選擇平方拋物線形斷面還是2.5次方拋物線形斷面或其他拋物線形斷面。不同次方的拋物線形斷面均有不同的水力學(xué)特點(diǎn)和適宜條件。

    本文選擇2.5次方拋物線形而不是更高次拋物線形,主要是其相對(duì)于平方拋物線形斷面,其水力條件有改善,但邊坡等方面變化不是很大,因此有較好的實(shí)用性,其更多性質(zhì)及與其他拋物線形斷面的比較需要進(jìn)一步研究。

    [1] Anwar A A, de Vries T T. Hydraulically efficient power-law channels[J]. Journal of Irrigation and Drainage Engineering, 2003, 129(1): 18-26.

    [2] Mironenko A P, Willardson L S, Jenab S A J. Parabolic canal design and analysis[J]. Journal of Irrigation and Drainage Engineering, 1984, 110(2): 241-246.

    [3] Das J A. Optimal design of channel having horizontal bottom and parabolic sides[J]. Journal of Irrigation and Drainage Engineering, 2007; 133(2): 192-197.

    [4] Easa S M. Improved channel cross section with two-segment parabolic sides and horizontal bottom[J]. Journal of Irrigation and Drainage Engineering, 2009, 135(3): 357-365.

    [5] Loganathan G V. Optimal design of parabolic canals[J]. Journal of Irrigation and Drainage Engineering, 1991, 117(5): 716-735.

    [6] Chahar J, Bhagu R. Optimal design of a special class of curvilinear bottomed channel section[J]. Journal of Hydraulic Engineering, 2007, 133(5): 571-576.

    [7] Babaeyan-Koopaei K, Valentine E M, Swailes, D C J. Optimal design of parabolic-bottomed triangle canals[J]. Journal of Irrigation and Drainage Engineering, 2000,126 (6): 408-411.

    [8] Han Yancheng. Horizontal bottomed semi-cubic parabolic channel and best hydraulic section[J]. Flow Measurement and Instrumentation, 2015, 45: 56-61.

    [9] 張新燕,呂宏興. 拋物線形斷面渠道正常水深的顯式計(jì)算[J]. 農(nóng)業(yè)工程學(xué)報(bào),2012,28(21):121-125. Zhang Xinyan, Lü Hongxing. Explicit solution for normal depth in parabolic-shape channel[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2012, 28(21): 121-125. (in Chinese with English abstracts)

    [10] 文輝,李風(fēng)玲. 平底馬蹄形斷面的水力計(jì)算[J]. 農(nóng)業(yè)工程學(xué)報(bào),2013,29(10):130-135.Wen Hui, Li Fengling. Hydraulic calculation of horseshoe cross-section with flat-bottom[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2013, 29(10): 130-135. (in Chinese with English abstract)

    [11] 王正中,王羿,趙延風(fēng),等. 拋物線形斷面河渠收縮水深的直接計(jì)算公式[J]. 武漢大學(xué)學(xué)報(bào):工學(xué)版,2011,44(2):175-177,191. Wang Zhengzhong, Wang Yi, Zhao Yanfeng, et al. Formula for direct calculation of contracted depth of parabola shaped canal[J]. Engineering Journal of Wuhan University, 2011, 44(2): 175-177, 191. (in Chinese with English abstracts)

    [12] 文輝,李風(fēng)玲. 數(shù)值積分法計(jì)算拋物線形渠道恒定漸變流水面線[J],農(nóng)業(yè)工程學(xué)報(bào),2014,30(24):82-86. Wen Hui, Li Fengling. Numerical integration method for calculating water surface profile of gradually varied steady flow in parabola shaped channel[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2014, 30(24): 82-86. (in Chinese with English abstracts)

    [13] Chow V T. Open Channel Hydraulics[M]. New York, N Y, United States: McGraw-Hill, 1959.

    [14] Chanson H. The Hydraulics of Open Channel Flow: An Introduction[M]. London, UK: Edward Arnold, 1999.

    [15] Easa S M. Versatile general elliptic open channel cross section[J]. KSCE Journal of Civil Engineering, 2016, 20(4): 1572-1581.

    [16] Easa S M. New and improved channel cross section with piece-wise linear or smooth sides[J]. Canadian Journal of Civil Engineering, 2011, 38(6): 690-697.

    [17] Flynn Lawrence E, Marino Miguel A. Canal design: Optimal

    cross sections[J]. Journal of Irrigation and Drainage Engineering, 1987, 113(3): 335-355.

    [18] Hussein A S. Simplified design of hydraulically efficient power-law channels with freeboard[J]. Journal of Irrigation and Drainage Engineering, 2008, 134(3): 380-386.

    [19] Bhattacharjya R K, Satish M G. Optimal design of a stable trapezoidal cross section using hybrid optimization techniques[J]. Journal of Irrigation and Drainage Engineering, 2007, 133(4): 323-329.

    [20] Aksoy B, Altan-Sakarya A. Optimal lined channel design[J]. Canadian Journal of Civil Engineering, 2006, 33(5): 535-545.

    [21] Guo C, Hughes W. Optimal channel cross section with freeboard[J]. Journal of Irrigation and Drainage Engineering, 1984,110(3): 304-314.

    [22] Han Yancheng, Said M Easa. Superior cubic channel section and analytical solution of best hydraulic properties[J]. Flow Measurement and Instrumentation, 2015, 50:169-177.

    [23] Abdulrahman A. Best hydraulic section of a composite channel[J]. Journal of Hydraulic Engineering, 2007, 133(6): 695-697.

    [24] Abramowitz M, Stegun I A. Handbook of Mathematical Functions[M]. New York: Dover Publications, 1964.

    [25] 裘松良. Gauss超幾何函數(shù)的導(dǎo)數(shù)和廣義Legendre關(guān)系[J]. 杭州電子科技大學(xué)學(xué)報(bào),1998,18(2):1-8. Qiu Songliang. Derivatives and generalized Legend’s relation of the Gaussian hypergeometric function[J]. Journal of the Hanzhou Institute of Electronic Engineering, 1998, 18(2): 1-8. (in Chinese with English abstract)

    [26] Vatankhaha Ali R, Said M Easa. Explicit solutions for critica and normal depths in channels with different shapes[J]. Flow Measurement and Instrumentation, 2011, 22: 43-49.

    [27] 張寬地,呂宏興,趙延風(fēng). 明流條件下圓形隧洞正常水深與臨界水深的直接計(jì)算[J], 農(nóng)業(yè)工程學(xué)報(bào),2009,25(3):1-5. Zhang Kuandi, Lü Hongxing, Zhao Yanfeng. Direct calculation for normal depth and critical depth of circular section tunnel under free flow[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2009, 25(3): 1-5. (in Chinese with English abstract)

    [28] 張新燕,呂宏興,朱德蘭. U 形渠道正常水深的直接水力計(jì)算公式[J]. 農(nóng)業(yè)工程學(xué)報(bào),2013,29(14):115-119. Zhang Xinyan, Lü Hongxing, Zhu Delan. Direct calculation formula for normal depth of U-shaped channel[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2013, 29(14): 115-119. (in Chinese with English abstract)

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

    Design of two and a half parabola-shaped canal and its effect in improving hydraulic property

    Han Yancheng1, Xu Zhenghe1, Gao Xueping2, Said M. Easa3

    (1.250022,; 2.300072,; 3.)

    Shapes of canal cross sections affect their discharge capacity, water depth and construction cost. Researches have shown that the curve-shaped canal such as quadratic parabolic and semi-cubic parabolic shape has good hydraulic property. However, the less smooth base of the quadratic parabolic and semi-cubic parabolic shape canal can affect the discharge capacity. To improve the hydraulic property and increase the discharge of the quadratic parabolic sections of canals, a section with two and a half parabola shape was proposed in this paper. Formulas for flow area, shape factor and water surface width for this new section were derived. The theoretical formula for the wetted perimeter was deduced using Gauss super-geometric functions. A model of the optimum hydraulic section that minimized the flow area for a given discharge was developed based on the Manning formula. The partial differential equation for the optimum hydraulic section was deduced using Lagrange’s multiplier method. After substituting the derivatives of wetted perimeter and flow area with respect to water depth and water surface width into this partial differential equation, the optimum model was successfully converted into an equation about the water surface width-depth ratio. Various explicit formulas to compute the characteristic’s parameters such as wetted perimeter, shape factor, flow area, normal water depth and critical water depth for the best hydraulic section were obtained. Using these formulas, the hydraulic design could be achieved easily. The results showed that the best ratio of water surface width-depth ratio for the optimum hydraulic section of the two and a half parabola-shaped canal was a constant (2.0 883). The two and a half parabola-shaped canal had better hydraulic properties than that with quadratic or semi-cubic parabolic sections. Comparisons with quadratic and semi-cubic parabolic sections showed that the flow discharge of the two and a half parabola-shaped section was the largest under the same water depth, which means it is an economical section. Under the same discharge, the water depth of the two and a half parabola-shaped section was smaller than the quadratic parabolic and semi-cubic parabolic sections. The flow area, wetted perimeter and water surface width of the two and a half parabola-shaped section was the least under the same discharge among the three sections. Minimum wetted perimeter and flow area implied that the cost of construction (excavation and lining cost) was minimized. In theory, the comparisons with quadratic and semi-cubic parabolic sections also showed that the construction cost of the proposed best hydraulic section was the lowest under the same discharge. To aid practical use, the 3- point and 4-point method of Gauss-Legendre approximate algorithm were presented for the wetted perimeter calculation. The application example with the water depth of 1.0-3.5 m showed that the approximate algorithm was highly accurate. The 3-point approximate format formula could meet the practical use and design with the maximum absolute error of 0.004 01 m. The results from the 4-point format formula almost equaled to those of the theoretical results with the maximum absolute error of 0.000 97 m. This research provides a theoretical basis for the design of the two and a half parabola-shaped canals with improved hydraulic properties.

    canals; hydraulics; design; open canal; two and a half parabola-shaped section; optimum hydraulic section; wetted perimeter

    10.11975/j.issn.1002-6819.2017.04.019

    TV131.4

    A

    1002-6819(2017)-04-0131-06

    2016-04-27

    2016-09-10

    國家“十二五”科技支撐計(jì)劃(2015BAB07B02-6);山東省重點(diǎn)研發(fā)計(jì)劃(2016GSF117038);南市科技發(fā)展計(jì)劃(201302052)。

    韓延成,男,甘肅武威人,副教授,博士,碩士生導(dǎo)師,主要從事水力學(xué)及河流動(dòng)力學(xué)方面的研究。濟(jì)南 濟(jì)南大學(xué)資源與環(huán)境學(xué)院,250022。Email:stu_hanyc@ujn.edu.cn

    韓延成,徐征和,高學(xué)平,Said M. Easa. 二分之五次方拋物線形明渠設(shè)計(jì)及提高水力特性效果[J]. 農(nóng)業(yè)工程學(xué)報(bào),2017,33(4):131-136. doi:10.11975/j.issn.1002-6819.2017.04.019 http://www.tcsae.org

    Han Yancheng, Xu Zhenghe, Gao Xueping, Said M. Easa. Design of two and a half parabola-shaped canal and its effect in improving hydraulic property[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2017, 33(4): 131-136. (in Chinese with English abstract) doi:10.11975/j.issn.1002-6819.2017.04.019 http://www.tcsae.org

    猜你喜歡
    次方過流水深
    書法靜水深流
    河北水利(2022年10期)2022-12-29 11:48:12
    基于水深分段選擇因子的多光譜影像反演水深
    變壓器相間過流保護(hù)不正確動(dòng)作原因的探討
    電子制作(2017年19期)2017-02-02 07:08:58
    國產(chǎn)500kW短波發(fā)射機(jī)過流保護(hù)電路的應(yīng)用缺陷與改進(jìn)
    電子制作(2016年11期)2016-11-07 08:43:36
    手表+手鏈+戒指 N次方組合
    Coco薇(2016年7期)2016-06-28 02:09:09
    水輪機(jī)過流部件改造與節(jié)能增效
    一組計(jì)算題的啟示
    GPS RTK技術(shù)在水深測(cè)量中的應(yīng)用
    配網(wǎng)雙端閉鎖式過流保護(hù)的實(shí)現(xiàn)方法
    浸入式水深監(jiān)測(cè)儀器的設(shè)計(jì)
    国产精品亚洲一级av第二区| 窝窝影院91人妻| 制服诱惑二区| a级毛片在线看网站| 亚洲九九香蕉| a级毛片a级免费在线| 中文字幕人成人乱码亚洲影| 韩国精品一区二区三区| 精品高清国产在线一区| 亚洲欧洲精品一区二区精品久久久| 欧洲精品卡2卡3卡4卡5卡区| 男人的好看免费观看在线视频 | 天堂影院成人在线观看| 亚洲人成电影免费在线| 国产精品av久久久久免费| 91成年电影在线观看| 黄色视频不卡| 国产av不卡久久| 香蕉av资源在线| 午夜福利18| 老汉色∧v一级毛片| 在线观看免费日韩欧美大片| 国产91精品成人一区二区三区| 欧美在线黄色| 日韩精品青青久久久久久| 在线观看66精品国产| 男男h啪啪无遮挡| 精品一区二区三区av网在线观看| 欧美乱妇无乱码| 亚洲av熟女| 精品久久久久久久毛片微露脸| 国产精品久久久久久精品电影 | 午夜精品久久久久久毛片777| 亚洲一区高清亚洲精品| 免费一级毛片在线播放高清视频| 国产成人欧美在线观看| 久久久久久国产a免费观看| 国产亚洲精品第一综合不卡| 国产三级在线视频| 村上凉子中文字幕在线| 亚洲成人精品中文字幕电影| 亚洲av成人不卡在线观看播放网| 日韩免费av在线播放| 日韩免费av在线播放| av中文乱码字幕在线| 国产精品影院久久| 婷婷精品国产亚洲av在线| 国产精品一区二区精品视频观看| 午夜激情av网站| 亚洲午夜精品一区,二区,三区| 中文资源天堂在线| 精品卡一卡二卡四卡免费| 欧美另类亚洲清纯唯美| 久久伊人香网站| 高清在线国产一区| 亚洲欧洲精品一区二区精品久久久| 国产单亲对白刺激| 真人一进一出gif抽搐免费| 91av网站免费观看| 久久香蕉国产精品| 亚洲av美国av| 久久精品aⅴ一区二区三区四区| avwww免费| 午夜福利欧美成人| 人人妻人人澡欧美一区二区| 久久久久久免费高清国产稀缺| 午夜两性在线视频| 麻豆久久精品国产亚洲av| 国产一级毛片七仙女欲春2 | 精品久久久久久成人av| 亚洲午夜精品一区,二区,三区| 99国产极品粉嫩在线观看| 亚洲男人天堂网一区| 欧美成人性av电影在线观看| 性欧美人与动物交配| 怎么达到女性高潮| 熟女电影av网| 午夜免费观看网址| 中亚洲国语对白在线视频| e午夜精品久久久久久久| 亚洲一区高清亚洲精品| 91九色精品人成在线观看| 黑人操中国人逼视频| 中文字幕高清在线视频| 男人舔女人的私密视频| 国产一区二区激情短视频| 久久天躁狠狠躁夜夜2o2o| 国产成+人综合+亚洲专区| 欧美亚洲日本最大视频资源| 欧美黄色片欧美黄色片| 这个男人来自地球电影免费观看| 日韩欧美在线二视频| 一本综合久久免费| 色在线成人网| 欧美日韩中文字幕国产精品一区二区三区| 国产欧美日韩一区二区精品| 久久午夜综合久久蜜桃| 757午夜福利合集在线观看| 男人舔奶头视频| 他把我摸到了高潮在线观看| 中国美女看黄片| 成熟少妇高潮喷水视频| 又黄又粗又硬又大视频| 好男人电影高清在线观看| 欧美黄色片欧美黄色片| 精品久久久久久久末码| 99国产精品99久久久久| 嫩草影视91久久| a级毛片a级免费在线| 免费在线观看亚洲国产| 亚洲国产日韩欧美精品在线观看 | 草草在线视频免费看| 久久久水蜜桃国产精品网| 在线观看午夜福利视频| 亚洲国产中文字幕在线视频| 欧美日韩瑟瑟在线播放| 亚洲片人在线观看| 丝袜人妻中文字幕| 成人18禁高潮啪啪吃奶动态图| 精品国内亚洲2022精品成人| 成人av一区二区三区在线看| 久久精品影院6| 久久久久国内视频| 热99re8久久精品国产| 窝窝影院91人妻| 动漫黄色视频在线观看| 1024香蕉在线观看| 高清毛片免费观看视频网站| 最新美女视频免费是黄的| 午夜激情av网站| 人人妻人人澡欧美一区二区| 听说在线观看完整版免费高清| 操出白浆在线播放| 黄色a级毛片大全视频| 国产精品98久久久久久宅男小说| av在线播放免费不卡| a级毛片a级免费在线| 中文字幕人妻熟女乱码| 免费一级毛片在线播放高清视频| 天天躁夜夜躁狠狠躁躁| 一二三四社区在线视频社区8| 18美女黄网站色大片免费观看| 99久久综合精品五月天人人| 岛国在线观看网站| 婷婷亚洲欧美| av电影中文网址| 欧美久久黑人一区二区| 法律面前人人平等表现在哪些方面| 午夜久久久久精精品| 九色国产91popny在线| 观看免费一级毛片| 久久 成人 亚洲| 成人国产综合亚洲| 国产黄片美女视频| 在线观看66精品国产| 欧美不卡视频在线免费观看 | 一个人观看的视频www高清免费观看 | 国内精品久久久久久久电影| 国内精品久久久久精免费| 亚洲中文字幕一区二区三区有码在线看 | 69av精品久久久久久| 国产亚洲欧美98| 99久久99久久久精品蜜桃| 丝袜在线中文字幕| 久久精品影院6| 日韩大码丰满熟妇| 国产又黄又爽又无遮挡在线| 黄网站色视频无遮挡免费观看| 日韩欧美一区视频在线观看| 99国产精品一区二区蜜桃av| 亚洲精品国产一区二区精华液| 国产av一区在线观看免费| 国产精品国产高清国产av| 久久国产乱子伦精品免费另类| 国产一区二区三区在线臀色熟女| 人妻久久中文字幕网| 日韩国内少妇激情av| 男女午夜视频在线观看| 一个人免费在线观看的高清视频| 免费在线观看黄色视频的| 午夜日韩欧美国产| 午夜福利欧美成人| 两性夫妻黄色片| 啦啦啦观看免费观看视频高清| av电影中文网址| 美女国产高潮福利片在线看| 男女午夜视频在线观看| 久久久久久久精品吃奶| 成人18禁高潮啪啪吃奶动态图| 亚洲欧美精品综合一区二区三区| 亚洲熟妇中文字幕五十中出| 午夜精品久久久久久毛片777| 级片在线观看| 国产一级毛片七仙女欲春2 | 免费在线观看日本一区| 国产一区二区三区视频了| 长腿黑丝高跟| 在线av久久热| 亚洲午夜理论影院| 国产精品野战在线观看| 日韩高清综合在线| 黑人欧美特级aaaaaa片| 国产亚洲av高清不卡| avwww免费| 亚洲国产精品成人综合色| 99在线视频只有这里精品首页| 国产人伦9x9x在线观看| 成人国产综合亚洲| 久久国产乱子伦精品免费另类| 成人三级黄色视频| 欧美精品啪啪一区二区三区| 成人特级黄色片久久久久久久| 老汉色av国产亚洲站长工具| 99精品久久久久人妻精品| 一区二区三区国产精品乱码| 精华霜和精华液先用哪个| 久久久久久久精品吃奶| 成人18禁高潮啪啪吃奶动态图| 久久婷婷成人综合色麻豆| 美女免费视频网站| 亚洲精品美女久久久久99蜜臀| 一个人免费在线观看的高清视频| 亚洲av电影不卡..在线观看| 亚洲国产高清在线一区二区三 | 久久国产精品男人的天堂亚洲| 日韩av在线大香蕉| 高清在线国产一区| 久久久久久久久中文| 精品欧美国产一区二区三| 欧美丝袜亚洲另类 | 精品福利观看| 欧美午夜高清在线| 大香蕉久久成人网| 一级毛片女人18水好多| 999精品在线视频| 老汉色av国产亚洲站长工具| 国产视频内射| 在线播放国产精品三级| 成人精品一区二区免费| 韩国av一区二区三区四区| 亚洲午夜精品一区,二区,三区| 亚洲中文字幕日韩| 999久久久国产精品视频| 老司机靠b影院| 一个人免费在线观看的高清视频| 亚洲欧美激情综合另类| 一级毛片女人18水好多| 男人舔奶头视频| 欧美激情极品国产一区二区三区| videosex国产| 日韩欧美国产在线观看| 亚洲专区中文字幕在线| 成在线人永久免费视频| 亚洲国产中文字幕在线视频| 夜夜看夜夜爽夜夜摸| www日本在线高清视频| 亚洲精品久久成人aⅴ小说| 一个人观看的视频www高清免费观看 | 精品久久久久久成人av| 熟女电影av网| 国产真人三级小视频在线观看| 久久久精品欧美日韩精品| 啪啪无遮挡十八禁网站| 国产亚洲精品一区二区www| 老熟妇仑乱视频hdxx| 90打野战视频偷拍视频| 人成视频在线观看免费观看| 国产亚洲欧美精品永久| 看免费av毛片| 波多野结衣高清无吗| 他把我摸到了高潮在线观看| 免费搜索国产男女视频| 高潮久久久久久久久久久不卡| 妹子高潮喷水视频| 男女之事视频高清在线观看| 我的亚洲天堂| av中文乱码字幕在线| 精品免费久久久久久久清纯| 欧美日韩瑟瑟在线播放| 久久国产乱子伦精品免费另类| 国产91精品成人一区二区三区| 岛国视频午夜一区免费看| 国产亚洲欧美在线一区二区| 亚洲精品久久国产高清桃花| 中文在线观看免费www的网站 | 一区二区三区精品91| 中文资源天堂在线| 亚洲中文av在线| 久久精品国产亚洲av香蕉五月| 无人区码免费观看不卡| 亚洲国产欧美日韩在线播放| 男人的好看免费观看在线视频 | 亚洲国产精品999在线| 欧美色视频一区免费| 午夜激情av网站| 亚洲性夜色夜夜综合| 亚洲人成网站高清观看| 一区二区日韩欧美中文字幕| 99在线视频只有这里精品首页| 免费无遮挡裸体视频| 极品教师在线免费播放| 人人妻人人看人人澡| 中文字幕人妻丝袜一区二区| 天堂影院成人在线观看| 久热这里只有精品99| 国产97色在线日韩免费| 亚洲av成人不卡在线观看播放网| 国产真实乱freesex| 一级a爱视频在线免费观看| 俄罗斯特黄特色一大片| x7x7x7水蜜桃| 亚洲男人的天堂狠狠| 这个男人来自地球电影免费观看| 十分钟在线观看高清视频www| 亚洲中文日韩欧美视频| 成人三级黄色视频| 久久精品91无色码中文字幕| 亚洲欧洲精品一区二区精品久久久| 老司机靠b影院| 久久久久久大精品| 特大巨黑吊av在线直播 | 97人妻精品一区二区三区麻豆 | av片东京热男人的天堂| 国产伦一二天堂av在线观看| 午夜福利在线在线| 中文字幕另类日韩欧美亚洲嫩草| 在线观看免费午夜福利视频| 欧洲精品卡2卡3卡4卡5卡区| www日本黄色视频网| 成年免费大片在线观看| 男女午夜视频在线观看| 老司机午夜十八禁免费视频| 男人舔奶头视频| 国产欧美日韩一区二区精品| 12—13女人毛片做爰片一| 十八禁人妻一区二区| 777久久人妻少妇嫩草av网站| 亚洲欧美一区二区三区黑人| 超碰成人久久| 51午夜福利影视在线观看| 色播在线永久视频| or卡值多少钱| 亚洲一区二区三区不卡视频| 成在线人永久免费视频| 国产亚洲精品av在线| 国产精品九九99| 人人澡人人妻人| 免费av毛片视频| 国内久久婷婷六月综合欲色啪| 久久国产精品男人的天堂亚洲| 午夜激情福利司机影院| 日韩欧美国产一区二区入口| 99热6这里只有精品| 少妇粗大呻吟视频| 国产真实乱freesex| 国产亚洲av高清不卡| 母亲3免费完整高清在线观看| 欧美性长视频在线观看| 成人午夜高清在线视频 | av视频在线观看入口| 女人爽到高潮嗷嗷叫在线视频| 亚洲av日韩精品久久久久久密| 丁香六月欧美| 国产一区二区三区在线臀色熟女| 9191精品国产免费久久| 亚洲九九香蕉| 国产免费男女视频| 性色av乱码一区二区三区2| 国产成+人综合+亚洲专区| 国产一区二区三区在线臀色熟女| 国产精品亚洲美女久久久| 嫩草影院精品99| 热re99久久国产66热| 亚洲精品国产精品久久久不卡| 国产99久久九九免费精品| 国产精品一区二区三区四区久久 | 黄色女人牲交| 久久国产精品人妻蜜桃| 日本精品一区二区三区蜜桃| 日日爽夜夜爽网站| 欧美日韩中文字幕国产精品一区二区三区| 一区二区三区国产精品乱码| 成人免费观看视频高清| 美国免费a级毛片| 极品教师在线免费播放| 欧美精品啪啪一区二区三区| 人成视频在线观看免费观看| 精品少妇一区二区三区视频日本电影| 高清毛片免费观看视频网站| 久久久久久免费高清国产稀缺| 一本一本综合久久| 成人av一区二区三区在线看| 一进一出抽搐动态| 2021天堂中文幕一二区在线观 | 中文字幕人成人乱码亚洲影| 99热这里只有精品一区 | 白带黄色成豆腐渣| 国产精品久久视频播放| 精品国产亚洲在线| 国产精品美女特级片免费视频播放器 | 欧美乱妇无乱码| 妹子高潮喷水视频| 99久久国产精品久久久| 精品久久久久久,| 亚洲第一欧美日韩一区二区三区| 久久狼人影院| 亚洲第一电影网av| 欧美黄色片欧美黄色片| 成人一区二区视频在线观看| 午夜福利免费观看在线| 色播在线永久视频| 51午夜福利影视在线观看| 曰老女人黄片| 禁无遮挡网站| 好男人电影高清在线观看| 三级毛片av免费| av天堂在线播放| 久久久久久亚洲精品国产蜜桃av| 久久国产精品男人的天堂亚洲| 国产激情久久老熟女| 久久精品国产亚洲av香蕉五月| 中文字幕最新亚洲高清| 19禁男女啪啪无遮挡网站| 午夜日韩欧美国产| 黄频高清免费视频| 国语自产精品视频在线第100页| 最近最新免费中文字幕在线| 亚洲精品国产区一区二| 美女 人体艺术 gogo| 人人妻,人人澡人人爽秒播| 欧美黄色淫秽网站| 在线观看一区二区三区| 制服人妻中文乱码| 在线天堂中文资源库| 男男h啪啪无遮挡| 成在线人永久免费视频| 国产欧美日韩一区二区精品| 黄色 视频免费看| 久久性视频一级片| 妹子高潮喷水视频| 三级毛片av免费| 亚洲欧美日韩无卡精品| 日本在线视频免费播放| 久久精品影院6| 真人做人爱边吃奶动态| 久久中文字幕一级| 婷婷亚洲欧美| 黄片大片在线免费观看| 国产爱豆传媒在线观看 | 亚洲国产日韩欧美精品在线观看 | 午夜免费观看网址| 中文字幕久久专区| 听说在线观看完整版免费高清| 中文字幕精品免费在线观看视频| 久久亚洲真实| 悠悠久久av| 中文字幕精品亚洲无线码一区 | 黄色片一级片一级黄色片| 久久人人精品亚洲av| 欧美成人免费av一区二区三区| 久久精品aⅴ一区二区三区四区| 中文在线观看免费www的网站 | 黄色片一级片一级黄色片| 亚洲自拍偷在线| 亚洲精品国产一区二区精华液| 在线免费观看的www视频| 脱女人内裤的视频| 国产亚洲精品久久久久5区| 精品福利观看| 一级作爱视频免费观看| 亚洲国产精品sss在线观看| 国产激情久久老熟女| 欧美黑人欧美精品刺激| 在线国产一区二区在线| 日日干狠狠操夜夜爽| 日本成人三级电影网站| 国产真实乱freesex| 欧美日本亚洲视频在线播放| 亚洲国产中文字幕在线视频| 男女床上黄色一级片免费看| 女人爽到高潮嗷嗷叫在线视频| 一个人观看的视频www高清免费观看 | 1024手机看黄色片| av天堂在线播放| 欧洲精品卡2卡3卡4卡5卡区| 变态另类丝袜制服| 又黄又粗又硬又大视频| 午夜免费激情av| 999久久久国产精品视频| 久久久久久亚洲精品国产蜜桃av| 精品国产美女av久久久久小说| 成人欧美大片| 日本a在线网址| 夜夜看夜夜爽夜夜摸| 午夜激情福利司机影院| 草草在线视频免费看| 午夜免费鲁丝| 日本免费a在线| 一级黄色大片毛片| 99热6这里只有精品| 国产主播在线观看一区二区| 日日爽夜夜爽网站| 亚洲aⅴ乱码一区二区在线播放 | АⅤ资源中文在线天堂| 一本大道久久a久久精品| a级毛片在线看网站| 久久久久久人人人人人| 中出人妻视频一区二区| 亚洲久久久国产精品| 国产亚洲精品av在线| 麻豆av在线久日| 成人国产一区最新在线观看| √禁漫天堂资源中文www| 亚洲av成人av| 熟妇人妻久久中文字幕3abv| 俺也久久电影网| 欧美又色又爽又黄视频| www.熟女人妻精品国产| 一个人免费在线观看的高清视频| 欧美不卡视频在线免费观看 | 国产激情欧美一区二区| 51午夜福利影视在线观看| 妹子高潮喷水视频| 精品久久蜜臀av无| 日本精品一区二区三区蜜桃| 99国产精品一区二区三区| 黄片大片在线免费观看| 亚洲成人精品中文字幕电影| 日本一本二区三区精品| 国内毛片毛片毛片毛片毛片| 久久香蕉精品热| 亚洲人成伊人成综合网2020| 欧美午夜高清在线| 亚洲午夜理论影院| 国产精品久久视频播放| 婷婷精品国产亚洲av在线| 国产久久久一区二区三区| 视频区欧美日本亚洲| 亚洲专区国产一区二区| 国产日本99.免费观看| 国产精品久久久人人做人人爽| 午夜日韩欧美国产| 99riav亚洲国产免费| 精品人妻1区二区| 国产aⅴ精品一区二区三区波| 一个人免费在线观看的高清视频| 黑人巨大精品欧美一区二区mp4| 中文字幕精品亚洲无线码一区 | 欧美成人免费av一区二区三区| 老汉色av国产亚洲站长工具| 又黄又爽又免费观看的视频| 精品熟女少妇八av免费久了| 男人舔女人下体高潮全视频| 丁香六月欧美| 麻豆av在线久日| 最好的美女福利视频网| 91大片在线观看| 午夜亚洲福利在线播放| 亚洲中文字幕一区二区三区有码在线看 | 免费观看人在逋| 亚洲中文字幕日韩| 亚洲一区高清亚洲精品| 亚洲av成人av| 精品福利观看| 国产日本99.免费观看| 欧美激情 高清一区二区三区| 中文字幕精品免费在线观看视频| 麻豆国产av国片精品| 久久狼人影院| 亚洲国产欧美网| 18禁裸乳无遮挡免费网站照片 | av福利片在线| 久久精品aⅴ一区二区三区四区| 99久久99久久久精品蜜桃| 一区二区三区精品91| 巨乳人妻的诱惑在线观看| 精品国内亚洲2022精品成人| 啪啪无遮挡十八禁网站| 精品人妻1区二区| 91老司机精品| 亚洲黑人精品在线| 成人国产一区最新在线观看| 91大片在线观看| 亚洲欧美一区二区三区黑人| 男人舔女人的私密视频| 中文在线观看免费www的网站 | 亚洲精品国产一区二区精华液| 淫秽高清视频在线观看| 精品国产亚洲在线| 欧美日韩亚洲国产一区二区在线观看| 女生性感内裤真人,穿戴方法视频| 精品福利观看| 国产亚洲欧美98| 国产精品二区激情视频| 精品福利观看| 十分钟在线观看高清视频www| 哪里可以看免费的av片| 777久久人妻少妇嫩草av网站| xxxwww97欧美| 男女床上黄色一级片免费看| 亚洲色图av天堂| 夜夜夜夜夜久久久久| 中文字幕人妻丝袜一区二区| 国语自产精品视频在线第100页| 国产亚洲欧美98| 19禁男女啪啪无遮挡网站| 黑人巨大精品欧美一区二区mp4| 久久草成人影院| 国产一区二区激情短视频| 侵犯人妻中文字幕一二三四区| 狂野欧美激情性xxxx| 久久香蕉国产精品| 精品不卡国产一区二区三区| 欧美黄色淫秽网站| 女性生殖器流出的白浆|