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

    堰槽組合設施水力特性試驗與數(shù)值模擬

    2024-10-23 00:00:00張文科何新林李小龍曾鈺峰王佳鑫
    排灌機械工程學報 2024年10期
    關鍵詞:數(shù)值模擬

    摘要: 為探究堰槽組合設施的水力特性及適用性,通過模型試驗與數(shù)值模擬相結合的方法,設計 3 組坡度 (0.1%,0.4%,0.8%)、9 組流量進行水力性能試驗,驗證了數(shù)值模擬對于堰槽組合的適用性,探究了不同底坡下堰槽組合裝置的測流機制.結果表明:3 組坡度閾值相對水深為0.880,0.930和0.873,測流公式精度較高,對應的平均誤差分別為1.20%,1.69%和1.74%;不同坡度下,量水槽內水深沿程降低,中垂線平均流速沿程增加,流量小于閾值流量時,綜合流量系數(shù)和中垂線平均流速隨流量的增加而增大;流量大于閾值流量時,在過渡段上游1/3附近中垂線平均流速大小一致,流量系數(shù)和上游中垂線平均流速隨流量的增大而減小,下游隨流量的增大而增大;數(shù)值模擬結果最大誤差為4.16%,與試驗相吻合,數(shù)值模擬適用于堰槽組合量水設施的研究.

    關鍵詞: 堰槽組合;流量測量;數(shù)值模擬;水力特性;測流公式

    中圖分類號: S274.4 文獻標志碼: A 文章編號: 1674-8530(2024)10-1044-08

    DOI:10.3969/j.issn.1674-8530.23.0198

    張文科,何新林,李小龍,等. 堰槽組合設施水力特性試驗與數(shù)值模擬[J]. 排灌機械工程學報,2024,42(10):1044-1051.

    ZHANG Wenke, HE Xinlin, LI Xiaolong,et al. Experimental and numerical simulation study on hydraulic performance of weir-flume combination facility[J]. Journal of drainage and irrigation machinery engineering(JDIME), 2024, 42(10): 1044-1051. (in Chinese)

    Experimental and numerical simulation study on hydraulic

    performance of weir-flume combination facility

    ZHANG Wenke1,2, HE Xinlin1,2*, LI Xiaolong1,2, ZENG Yufeng1,2, WANG Jiaxin1,2

    (1. College of Water Conservancy amp; Architectural Engineering, Shihezi University, Shihezi, Xinjiang 832000, China; 2. Key Laboratory of Cold and Arid Regions Eco-Hydraulic Engineering of Xinjiang Production amp; Construction Corps, Shihezi, Xinjiang 832000, China)

    Abstract: In order to investigate the hydraulic properties and applicability of the weir-flume combination facilities, 3 groups of slopes (0.1%, 0.4%, 0.8%) and 9 groups of flow rates were designed to conduct hydraulic performance tests by combining modeling tests with numerical simulations. The applicability of numerical simulation for weir-flume combination was verified, and the flow measurement mechanism of weir-flume combination devices under different bottom slopes was explored. The results show that the relative water depths of three groups of slope thresholds are 0.880, 0.930 and 0.873, and the accuracy of the flow measurement formulas is high, with corresponding average errors of 1.20%, 1.69% and 1.74%, respectively. Under different slopes, the water depth in the measuring flume decreases along the course, and the average longitudinal velocity increases along the course. When the flow rate is less than the threshold flow rate, the comprehensive flow coefficient and the average longitudinal velocity increases with the increase of the flow rate. When the flow rate is greater than the threshold flow rate, the size is consistent near the average longitudinal velocity in the 1/3 upstream of the transition section, the flow coefficient and the average longitudinal velocity of the upstream decrease with the increase of the flow rate, while the average longitudinal velocity of the downstream increases with the increase of the flow rate. The maximum error of numerical simulation results is 4.16%, which is consistent with the test, and the numerical simulation is suitable for the study of weir-flume combination water measuring facilities.

    Key words: weir-flume combination;flow measurement;numerical simulation;hydraulic properties;flow-measuring formula

    新疆氣候干旱,降水稀少,是典型的“綠洲灌溉農業(yè)”,農業(yè)生產長期依賴灌溉,農業(yè)用水比例高,隨著用水需求增加,水資源供需不平衡日益突出,水資源緊缺制約著新疆農業(yè)的發(fā)展.因此,對用水量進行精準測量不僅是優(yōu)化水資源配置的重要基礎,而且為推進農業(yè)水價綜合改革提供了依據,減少了供水單位和用水組織之間的矛盾.

    西北地區(qū)渠道中泥沙含量高,大面積推行膜下滴灌技術,灌溉水中的泥沙易堵塞滴灌帶.量水堰作為各級渠道中廣泛應用的工程措施,通常用于治理河道泥沙.渠道中沉積物和漂浮物等雜質多、流量變幅大等問題使測流工作變得復雜,故需改進現(xiàn)有的測量結構以解決當?shù)氐膶嶋H問題并提高精度.PITON 等[1]發(fā)現(xiàn)量水堰不僅能夠穩(wěn)定河床,還能通過緩沖作用影響推移質的運輸和泥沙傳播規(guī)律.KALITA 等[2]通過試驗模擬克倫普堰和溢流堰的流動,檢驗了數(shù)值模型的質量.凌剛等[3]對不同出口寬度的堰槽組合設施水力性能進行研究,發(fā)現(xiàn)出口寬度對槽內的水流有明顯影響.隨后,他們還研究了臂坡對堰槽組合設施泄流能力的影響,結果表明增大臂坡可有效提高小流量工況下的測流精度[4].廖偉等[5]探究了臂坡對堰槽組合設施水流紊動特性的影響,發(fā)現(xiàn)臂坡能有效提高設施的消能率.JIN 等[6]對魚形量水槽在不同工況下進行了試驗和數(shù)值模擬,研究了便攜式魚形量水槽在“U”形渠道中的水力特性和測流性能.XU 等[7]評估了巴歇爾槽不同翼墻形式對測流精度和水力特性的影響.

    前人對量水堰和量水槽多是進行獨立研究[8-10〗,對于堰槽組合應用于灌區(qū)末級渠道量水的研究相對較少,缺乏數(shù)值模擬適用性的研究,且底坡對堰槽組合的水力性能和測流機理的影響尚不明確[11].文中研究基于模型試驗與數(shù)值模擬,對比不同工況下堰槽組合設施的水力特性,建立不同過流機制下流量與測點水深的關系,為堰槽組合設施的推廣及應用提供理論參考.

    1 材料與方法

    1.1 堰槽組合設施結構

    根據堰槽組合設施的結構特點,在克倫普堰中間加設量水槽.量水槽分為矩形段、過渡段和梯形喉段,槽底為渠底,量水槽深 15 cm.上游分水墻厚度為 2 cm,迎水面呈 45° 斜坡的半圓柱形;梯形喉段長15 cm;克倫普堰迎水坡面坡度為1∶2,背水坡面坡度為1∶5,具有水面線過渡平順的特點,堰高和槽深相等.堰槽組合設施通過 3D 打印技術采用樹脂材料制成,用玻璃膠固定在試驗渠道中.體型參數(shù)詳見圖 1.

    1.2 試驗裝置與試驗方案

    試驗于石河子大學水力大廳進行,試驗裝置的主要組成部分有地下水庫、供水系統(tǒng)、電磁流量計、調節(jié)閥門、上游穩(wěn)水池、渠道、水位測量系統(tǒng)、流速測量系統(tǒng)、堰槽組合設施、變坡系統(tǒng)等,試驗體系及測點布置平面圖如圖2所示.通過水泵將地下水庫中的水抽入供水管道,電磁流量計實時監(jiān)測流量,水流經過穩(wěn)水池和平水柵后平穩(wěn)地進入渠道.試驗渠道為矩形,尺寸為30.0 m×0.6 m×0.8 m,材質為鋼化玻璃,通過DCMS型流量控制系統(tǒng)調節(jié)流量.

    試驗設計 9 組不同的流量,流量范圍5~35 L/s,設置 0.1%,0.4%,0.8% 3種坡度.以量水槽進口上游190 cm處的斷面為起始斷面,分別在上游、槽內和下游區(qū)域布置控制斷面及測點,上游沿中軸線及兩側沿軸線共布置23個測點(L0—L7,M0—M6,R0—R7);量水槽內各段各設置3個控制斷面,共10個測點(M7—M16);下游設置17個測點(M17—M33),總計設置50個測點.控制斷面測點位置如表1所示,表中L為各測點距測點M0的距離.

    依據試驗布置圖在渠道中安裝堰槽組合設施,先通過變坡控制系統(tǒng)將渠道坡度調整為設計坡度,采用流量控制系統(tǒng)調節(jié)流量至設計流量,待流量穩(wěn)定后,通過自動水位測量系統(tǒng)采集各個測點的水位,并用 ADV 流速儀測量各測點的流速,每個測點進行 3 次測量取平均值.調節(jié)流量和坡度至不同工況,重復上述試驗過程.整個試驗均保證自由出流.

    2 結果與分析

    2.1 水力特性分析

    堰槽組合有2種過流方式:當流量較小時,水流全部從量水槽內部通過,稱為槽內流;隨著流量的增加,水深超過堰頂,量水槽和堰頂同時過流,稱為堰槽流,二者的臨界流量稱為流量閾值.槽內測點水深與過槽流量具有穩(wěn)定關系,用于確定通過水槽和相鄰堰的流量.該測點處水深與槽深的比值為相對水深[12],通過相對水深可判斷過流方式.不同坡度的流量閾值分別為 19.73,20.16,19.88 L/s(對應流量閾值相對水深分別為0.880,0.930,0.873).

    2.1.1 水面線

    水面線可以直觀反映水流沿程變化情況.通過試驗測量各測點水深,2種過流方式不同工況下試驗段中垂面水深沿程變化如圖3,h為中垂面上各測點的水深,Q為流量.

    由圖3可知,水面線沿程總體變化趨勢基本一致.量水槽進口過流斷面收縮,水流在上游產生壅水,渠道坡度越大,壅水越嚴重;水流進入量水槽矩形段,水位略微下降;水流進入過渡段,過流斷面逐漸減小,水面出現(xiàn)明顯下降;水流經過加速進入梯形喉部形成控制段,水面下降最快,流態(tài)從緩流變?yōu)榧绷?,并產生臨界流.水流流出量水槽后,過流斷面突然變大,量水槽對水流的側向約束消失,水流向兩側和下游沖擊形成薄水層,水面進一步降低,渠道側壁限制水流的側向流動,在下游形成中間淺、兩側深的現(xiàn)象.同一流量在不同坡度下,水面線在出口附近基本重合,坡度對于此段水流流態(tài)影響不大.渠道兩側的水流匯流后導致下游水深增加,水面線升高,水流由急流變?yōu)榫徚鳎诒∷畬雍笮纬伤S,此時水面運動劇烈,水流出現(xiàn)明顯的摻氣現(xiàn)象.薄水層前端水深隨流量的增大而增加.槽內流時,流量越大,發(fā)生水躍的位置越靠后,坡度只影響其下游水深;發(fā)生堰槽流時,流量和坡度對水躍的位置影響不大.

    2.1.2 弗勞德數(shù)

    弗勞德數(shù)Fr用來判別明渠水流流態(tài),根據 Fr 可判斷臨界流出現(xiàn)的位置[13].當 Fr=1 時,慣性力與重力作用相等,水流為臨界流;當 Fr<1 時,慣性力小于重力作用,水流為緩流.流態(tài)分為槽內流和堰槽流,對不同的過流機制采取不同的計算方式.槽內流時,斷面平均流速可通過式 (1) 計算;堰槽流時,采用流速儀測量得到的中垂線平均流速來計算 Fr,如式(2)所示.

    Fr=vgh=QAgh,(1)

    Fr=vvcrgh,(2)

    式中:v為斷面平均流速,m/s;h 為斷面平均水深,m;g為重力加速度,取9.81 m/s2;A為過流面積,m2;vvcr為中垂線平均流速,m/s.

    槽內流時,不同工況下Fr 沿程變化如圖4a所示,F(xiàn)r沿程總體變化趨勢相似.量水槽上游Fr沿程減小且保持在0.15~0.25,滿足測流規(guī)范,F(xiàn)r隨坡度的增加而增加;水流進入量水槽后,水面下降,勢能轉化為動能,F(xiàn)r增至0.30~0.40,流量越大,F(xiàn)r越大;進入過渡段后,F(xiàn)r增速加快;梯形段內Fr增速最大,水流由緩流轉為急流.根據量水規(guī)范要求,水位測點處需滿足Fr<0.50[13].由圖4可知,發(fā)生槽內流時,測點M11上游測點均滿足量水水位測點選取標準;發(fā)生堰槽流時,測點M10上游測點滿足水位測點要求.

    2.1.3 流速分布

    不同坡度下中垂線平均流速如圖5所示.

    垂線平均流速都具有沿程增大的規(guī)律.量水槽入口過流面積減小,流速增大,過渡段內水流不斷加速,在梯形段形成臨界流.槽內流時,流量越大,各測點流速也越大,坡度對其變化趨勢影響較小.當流量超過閾值流量后,不同流量下流速分布線在測點M11附近相交,即該測點附近流速相同,而在M7—M10,流速隨流量的增加而減小.水流在進口被分水墻分成3部分,隨著流量的增大,通過堰的流量增加,水深在堰前壅高,超過分水墻高程,部分流量躍過分水墻進入量水槽,使槽內流量沿程增加,流速逐漸增大.

    2.2 流量公式

    2.2.1 測流公式

    槽內流時,臨界流位于梯形段,臨界斷面具有穩(wěn)定的水深-流量關系.試驗結果表明,臨界斷面位置會隨流量變化,無法確定臨界斷面的準確位置,因此選擇上游水流平穩(wěn)的斷面作為測量斷面,通過能量方程建立穩(wěn)定的流量水深關系式.假設忽略測點斷面至臨界斷面的能量損失,結合伯努利方程可以得到槽內流的測流公式為

    Q=vcAc=ghcbhc=bghc3/2,(3)

    Q=2323bghi+Δz+v2i2g3/2,(4)

    式中:vc為臨界流速,m/s;Ac為臨界斷面面積,m2;hc為臨界水深,m;b為梯形段底寬,m;hi為上游測點水深,m;Δz為上游測點基于臨界斷面底部高程的位置水頭,m;vi為上游測點平均流速,m/s.

    在實際測量中無法確定流速水頭,因此采用水深代替總水頭.測量斷面總水頭包括位置水頭、壓力水頭及流速水頭,將總水頭轉換成hi.考慮測量斷面和臨界斷面間的能量損失,引入綜合流量系數(shù)Cdv對流量公式進行修正,由此得到槽內流的測流公式為

    Q=0.385Cdvb2ghi3/2.(5)

    堰槽流時,克倫普堰上游水深高于槽深,水流發(fā)生橫向流動進入量水槽,槽內流量沿程增大.因此,王文娥等[11]將堰槽流下的堰槽組合設施看作2個“L”形迷宮堰,提出堰頂溢流時堰槽組合設施的流量公式為

    Q=0.385m0Lw2gh3/2w,(6)

    式中:m0為堰槽流的綜合流量系數(shù);Lw為堰槽組合的有效臂長,m;hw為堰槽組合設施上游高于堰頂高程的水深,m.

    2.2.2 流量系數(shù)

    根據試驗數(shù)據,利用槽內流流量公式計算不同流量下各測點(M7—M10)的綜合流量系數(shù),分別對M7—M10斷面相對水深與綜合流量系數(shù)進行相關性分析,發(fā)現(xiàn)綜合流量系數(shù)與測點M9處的相對水深的相關性最好,3種坡度下的相關系數(shù)分別為0.997 5,0.995 0,0.991 9,擬合得到綜合流量系數(shù)的經驗公式為

    Cdv1=0.907 3+0.475 8(h9/d)+0.039 9(h9/d)2,(7)

    Cdv4=0.857 6+0.633 8(h9/d)-0.016 4(h9/d)2,(8)

    Cdv8=1.102 9+0.177 9(h9/d)+0.198 76(h9/d)2,(9)

    式中:Cdv1,Cdv4和Cdv8分別為坡度 0.1%,0.4%,0.8%流量閾值內的綜合流量系數(shù);h9為測點M9處的水深,m;d為槽深;h9/d為測點M9處的相對水深.

    槽內流時,擬合綜合流量系數(shù)與測點M9的相對水深得到擬合曲線,如圖6a所示,相同坡度下,綜合流量系數(shù)隨測點相對水深的增加而增大,在同一相對水深下,渠道坡度越大,流量系數(shù)越大.

    堰槽流時,選擇上游水流緩慢平穩(wěn)的測點M1的水深h1作為測點水深,擬合不同坡度下的流量系數(shù)與hw/h1,對應的相關系數(shù)分別為0.999 4,0.999 9和0.999 9,擬合得到流量系數(shù)的經驗公式為

    m0-1=0.329 34(hw/h1)-1.405 65,(10)

    m0-4=0.338 84(hw/h1)-1.416 74,(11)

    m0-8=0.329 34(hw/h1)-1.461 92,(12)

    式中:m0-1,m0-4和m0-8分別為0.1%,0.4%,0.8%坡度下的流量系數(shù).

    擬合不同坡度下流量系數(shù)m0與測點M1的相對水深hw/h1,得到擬合曲線如圖6b所示.由圖可知,不同坡度下,流量系數(shù)隨hw/h1的變化趨勢相近,流量系數(shù)的變化范圍較大;相同坡度下,流量系數(shù)隨相對水深的增加而減?。煌幌鄬λ钕?,坡度越大,流量系數(shù)越大.

    2.2.3 測點M1與測點M9的水深關系

    量水槽上游流速緩慢,水流中的泥沙等沉積物會在堰槽前沉積,使測點M1處水深的測量誤差較大,因此水深監(jiān)測需要選擇沒有沉積物的區(qū)域.在試驗中發(fā)現(xiàn)量水槽內水流流速增加,不會出現(xiàn)泥沙淤積.因此,堰槽流的水位測點位置應布設在槽內,以避免淤積對測流精度的影響.

    測點M9為槽內流的測點,且與M1水深之間相關性高,因此,將測點M1的水深轉化為M9的水深,測量中對測點M9進行水位監(jiān)測,即可獲得渠道中的流量.測點M1相對水深隨測點M9相對水深的增加而增加,具有很明顯的線性關系,3種坡度對應的相關系數(shù)分別為0.999 6,0.993 5和0.996 4,擬合得到不同坡度下水深經驗轉化公式為

    h1-1/d=0.122 73+0.919 61(h9-1/d),(13)

    h1-2/d=0.112 29+0.912 53(h9-2/d),(14)

    h1-3/d=0.147 65+0.845 75(h9-3/d),(15)

    式中:h1-1,h1-2,h1-3分別為0.1%,0.4%,0.8%坡度下測點M1水深;h9-1,h9-2,h9-3分別為0.1%,0.4%,0.8%坡度下測點M9水深.

    2.2.4 流量公式精度驗證

    將綜合流量系數(shù)的經驗公式代入槽內流流量公式,水位轉化公式和流量系數(shù)經驗公式代入堰槽流流量公式,試驗中的各工況下測點M9的水深代入流量公式,得到不同工況下的計算流量,并與實測流量對比發(fā)現(xiàn),計算流量與實測流量的最大誤差為7.39%,最小為0.05%,平均誤差為1.54%. 3組坡度(0.1%,0.4%,0.8%)對應的平均誤差分別為1.20%,1.69%和1.74%.渠道越緩,流量公式精度越高.結果表明流量公式計算精度較高,滿足灌區(qū)量水要求.

    3 堰槽數(shù)值模擬

    為進一步研究堰槽組合設施內部流場分布,深入探究其水力特性,利用 Fluent 對堰槽組合設施進行模擬分析,將數(shù)值模擬結果與模型試驗進行比對,驗證模擬的適用性及準確性.

    3.1 控制方程與網格及邊界條件

    明渠水流中流體為氣液兩相流,2種流體互不相溶,采用VOF模型模擬自由水面,計算域內每個網格單元中氣體和液體的體積分數(shù)和為1.根據模型試驗設計尺寸建立數(shù)值模型,由于試驗水深較淺,為提高計算速率,構建渠道模型高度為50 cm.重力加速度為-9.81 m/s2.明渠流動是充分發(fā)展的湍流,RNG k-ε模型能更好地處理有強旋流和彎曲流線流動的湍流計算問題,更加適用于明渠三維數(shù)值模擬.使用瞬態(tài)算法,時間步長0.001 s.采用混合網格對計算域進行網格劃分,上下游渠道采用結構化網格,堰槽組合段采用非結構化網格,同時對其進行局部加密.渠道的空氣進口設置為pressure-inlet,水體進口設置為velocity-inlet;渠道及堰槽壁面默認無滑移,設置為wall;水流出口設置為pressure-outlet,選擇PISO算法進行壓力-速度的耦合求解.

    3.2 模擬結果驗證

    對比不同工況下沿程水深及流速的模擬結果與實測數(shù)據,驗證模擬結果的可靠性.以0.1%坡度為例,結果如圖7所示,對比自由出流條件下沿程水深和流速的模擬值和實測值,可以看出,模擬結果與實測結果吻合,分析數(shù)據得出實測值與模擬值誤差最大為4.16%,平均誤差值為1.70%,誤差滿足小于5.00%,驗證了數(shù)值模擬的可靠性.

    3.3 流線特征及流速分布

    槽內流以10.02 L/s流量為例,堰槽流以25.00 L/s流量為例,得到不同過流情況下的流線及流速v分布如圖8所示,圖中I為渠道坡度.由圖可知,量水槽上游水流平緩,流線基本平行;槽內流時,水流進入量水槽,受分水墻的分割,一部分流向堰前產生壅水,此處流線較為紊亂;進口處過流斷面突變,水流發(fā)生繞流,流線出現(xiàn)明顯彎曲;水流在槽內匯集,流線逐漸收縮,密度變大,水流加快;水流流出量水槽后,流線向兩側分散.堰槽流時,堰槽段流線分為3部分,一部分直接從量水槽內部通過,一部分從堰頂通過,少部分躍過分水墻進入量水槽.最大流速出現(xiàn)在出口下游區(qū)域.槽內流時,水流在出口向兩側自由發(fā)散,最大流速出現(xiàn)位置靠后;堰槽流時,通過堰頂?shù)乃鲗α克哿鞒龅乃鳟a生擠壓,水流很快由急流變成緩流,阻礙出口處水流的擴散,使最大流速出現(xiàn)位置提前.根據堰槽橫斷面流速分布圖可以看出,流速在量水槽內沿程逐漸變得集中,呈現(xiàn)對稱分布.槽內流時,流速大的地方靠近槽底中間;堰槽流時,流速大的地方靠近槽底兩側,這是由于堰前水位高于分水墻,水流橫向流動沿槽體進入槽內,使槽內靠近槽體的流量增加,流速增大.

    4 討 論

    堰槽組合設施能夠在大范圍的流量下實現(xiàn)精準測量,具有較高的靈敏度.水流在量水槽內不斷加速,泥沙等在槽內沉積的風險降低,位于量水槽內的水位測量點能很好地避免沉積對測量精度的影響.在試驗中發(fā)現(xiàn),雖然在低流量條件下無法完全避免水槽內泥沙的沉積,但是沉積層相對較薄,因此不會對精度產生嚴重的影響.渠道上游停止放水,少量泥沙會在槽內沉積,但是在下次放水的時候,水流將重新懸浮槽內的沉積物,有較好的自清潔功能.流量大于閾值流量時,測點水位轉換為上游渠道中的相應水位,不受上游渠道中沉積物沉積的影響,可以更準確地估算上游水位[14].量水槽具有穩(wěn)定的校準特性,不受相鄰堰結構變化的影響;克倫普堰相對容易建造,堰體堅固,對堰頂三角形輪廓的輕微損傷不敏感.槽內流時,建立的流量公式不受堰形的影響;堰槽流時,需根據堰形對堰槽流的流量公式進行探討.

    隨著渠道坡度的增加,水流在堰槽上游壅水更加明顯,流量公式精度逐漸降低,因此在應用中為保證其自由出流的條件,需盡量選擇坡度較緩的渠道.

    文中研究是在理想的試驗室條件下進行,實際應用中的問題往往會導致整體精度降低,比如:水位測量或設施安裝精度差、淹沒出流、堰頂出現(xiàn)沉積等.沉積物對水位測量的誤差難以量化,因此,對堰槽組合設施應進行定期檢查和校核,最大限度地減小這種影響.

    5 結 論

    1) 在渠道坡度為 0.1%,0.4% 和 0.8% 時,堰槽的閾值流量分別為 19.73,20.16,19.88 L/s,對應流量閾值相對水深 0.880,0.930,0.873.

    2) 堰槽上游水深沿程增加,坡度越大,壅水高度越大;水流進入量水槽后,水深沿程下降,渠道坡度越大,水深越小.水流進入堰槽前, Fr 沿程減小,進入量水槽后 Fr 沿程增大,在過渡段和梯形段增速最大.流量越大,F(xiàn)r 值越大;量水槽內 Fr 值隨坡度的增加而增大,但坡度對 Fr 的影響較小.

    3) 槽內流的綜合流量系數(shù)隨坡度的增加而增大,堰槽流時,不同坡度下流量系數(shù)隨相對水深的變化總體趨勢接近,變化范圍較大,流量系數(shù)隨相對水深的增加而減小.

    參考文獻(References)

    [1] PITON G, RECKING A. Effects of check dams on bed-load transport and steep-slope stream morphodynamics [J]. Geomorphology, 2017, 291: 94-105.

    [2] KALITA H M, DAS R, HAJONG A, et al. Experi-mental and numerical flow simulation over weirs [J]. Water resources, 2019, 46(6): 934-943.

    [3] 凌剛, 王文娥, 王惠, 等. 量水槽出口寬度對堰槽組合量水設施過流能力的影響 [J]. 西北農林科技大學學報(自然科學版), 2023, 51(12):132-142.

    LING Gang, WANG Wen′e, WANG Hui, et al. In-fluence of outlet width of flume on overflow capacity of weir-flume combination facility [J]. Journal of North-west Aamp;F University(natural science edition), 2023, 51(12): 132-142.(in Chinese)

    [4] 凌剛, 廖偉, 程勇, 等. 臂坡對堰槽組合設施測流影響分析 [J]. 中國農業(yè)大學學報, 2023, 28(7): 204-214.

    LING Gang, LIAO Wei, CHENG Yong, et al. Influence analysis of crest slope on the flow measurement of weir-flume combination facility [J]. Journal of China Agricultural University, 2023, 28(7): 204-214.(in Chinese)

    [5] 廖偉, 張維樂, 王文娥, 等, 臂坡對堰槽組合設施紊流結構影響試驗研究 [J]. 水科學進展, 2023, 34(3): 465-479.

    LIAO Wei, ZHANG Weile, WANG Wen′e, et al. Experimental study on the effect of arm slope on turbulent structure of weir-flume combination [J]. Advances in water science, 2023, 34(3): 465-479.

    [6] JIN L, ZHANG C, ENGEL B A, et al. Hydraulic cha-racteristics and flow measurement performance analysis of portable fish-shaped flumes in U-shaped channels [J]. Flow measurement and instrumentation, 2023, 90: 102327.

    [7] XU H, LI Z M, WU W Y, et al. Hydraulic characteris-tics of the Parshall flume entrance wingwall: numerical simulation and test [J]. Irrigation and drainage, 2022, 71 (5): 1280-1293.

    [8] AALI F, VATANKHAH A R. Experimental study of simple flumes with trapezoidal contraction [J]. Flow measurement and instrumentation, 2023, 90: 102328.

    [9] 孫斌,朱順,楊磊,等.灌區(qū)明渠無喉道量水槽水力特性的影響因素研究[J].水利水電技術(中英文),2021,52(4):105-114.

    SUN Bin,ZHU Shun,YANG Lei,et al. Research on hydraulic characteristics influencing factors of open channel cut-throat flume in irrigation district[J]. Water resources and hydropower engineering,2021,52(4) : 105-114.(in Chinese)

    [10] 王蒙,張寬地,王文娥,等. 用于梯形渠道的仿翼形便攜式量水槽水力性能 [J]. 農業(yè)工程學報, 2023, 39(3): 76-83.

    WANG Meng, ZHANG Kuandi, WANG Wen′e, et al. Hydraulic performance of portable flow measuring flume of imitating aerofoil in the trapezoidal channel [J]. Transactions of the CSAE, 2023, 39(3): 76-83.(in Chinese)

    [11] 王文娥,廖偉, 陳土成, 等. 堰槽組合設施測流機制試驗研究 [J]. 水科學進展, 2021, 32(6): 922-932.

    WANG Wen′e, LIAO Wei, CHEN Tucheng, et al. Experimental study on the flow measurement mechanism of a weir-flume combination facility [J]. Advances in water science, 2021, 32(6): 922-932.(in Chinese)

    [12] ZUIKOV A L. Hydraulics of the classical crump weir water gauge [J]. Power technology and engineering, 2017, 50(6): 611-619.

    [13] 王長德. 量水技術與設施 [M]. 北京: 中國水利水電出版社, 2006.

    [14] WESSELS P, ROOSEBOOM A. Flow-gauging struc-tures in South African rivers Part 1: an overview [J]. Water SA, 2009, 35(1): 1-9.

    (責任編輯 黃鑫鑫)

    猜你喜歡
    數(shù)值模擬
    基于AMI的雙色注射成型模擬分析
    錐齒輪精密冷擺輾成形在“材料成型數(shù)值模擬”課程教學中的應用
    科教導刊(2016年28期)2016-12-12 06:22:00
    基于氣象信息及風場信息的風機輪轂處風速預測
    鉆孔灌注樁樁底沉渣對樁體承載特性影響的模擬分析
    西南地區(qū)氣象資料測試、預處理和加工研究報告
    科技資訊(2016年18期)2016-11-15 08:01:18
    張家灣煤礦巷道無支護條件下位移的數(shù)值模擬
    科技視界(2016年18期)2016-11-03 23:14:27
    張家灣煤礦開切眼錨桿支護參數(shù)確定的數(shù)值模擬
    科技視界(2016年18期)2016-11-03 22:57:21
    跨音速飛行中機翼水汽凝結的數(shù)值模擬研究
    科技視界(2016年18期)2016-11-03 20:38:17
    姚橋煤礦采空區(qū)CO2防滅火的數(shù)值模擬分析
    雙螺桿膨脹機的流場數(shù)值模擬研究
    科技視界(2016年22期)2016-10-18 14:53:19
    九九在线视频观看精品| 最好的美女福利视频网| 精品午夜福利在线看| 亚洲aⅴ乱码一区二区在线播放| 国产精品一二三区在线看| av在线老鸭窝| 99视频精品全部免费 在线| av国产免费在线观看| 亚洲国产欧美人成| 少妇熟女欧美另类| 伦精品一区二区三区| 神马国产精品三级电影在线观看| 国产精品一区二区三区四区久久| 能在线免费观看的黄片| 久久久a久久爽久久v久久| 国产精品1区2区在线观看.| 91av网一区二区| 黄色视频,在线免费观看| 春色校园在线视频观看| 99热精品在线国产| 97碰自拍视频| 偷拍熟女少妇极品色| 国产精品一区www在线观看| 成人美女网站在线观看视频| 中文亚洲av片在线观看爽| 国产精品av视频在线免费观看| 一夜夜www| 国产亚洲欧美98| 中文字幕av在线有码专区| 成年女人看的毛片在线观看| 天堂√8在线中文| 亚洲四区av| 久久久a久久爽久久v久久| 九九在线视频观看精品| 变态另类丝袜制服| 亚洲成人久久爱视频| 老女人水多毛片| 三级经典国产精品| 国产精品免费一区二区三区在线| 欧美色视频一区免费| 国产精品一区二区免费欧美| 老熟妇仑乱视频hdxx| 精品人妻偷拍中文字幕| 日韩制服骚丝袜av| 嫩草影院入口| 国产在线精品亚洲第一网站| 1000部很黄的大片| 国产精品爽爽va在线观看网站| 欧美xxxx性猛交bbbb| 大又大粗又爽又黄少妇毛片口| 伦理电影大哥的女人| 嫩草影视91久久| 午夜福利18| 我要看日韩黄色一级片| 又爽又黄无遮挡网站| 亚洲人成网站在线观看播放| 黄色日韩在线| 亚洲成人精品中文字幕电影| 丰满的人妻完整版| 免费人成视频x8x8入口观看| 久久久久精品国产欧美久久久| 国内精品一区二区在线观看| 91麻豆精品激情在线观看国产| 亚洲婷婷狠狠爱综合网| 97人妻精品一区二区三区麻豆| 亚洲内射少妇av| 国产片特级美女逼逼视频| 级片在线观看| 我的老师免费观看完整版| 亚洲熟妇中文字幕五十中出| 一进一出抽搐动态| 特级一级黄色大片| 亚洲欧美成人综合另类久久久 | 在线免费观看的www视频| 18禁裸乳无遮挡免费网站照片| 尤物成人国产欧美一区二区三区| 国内精品美女久久久久久| 久久精品国产清高在天天线| 九九热线精品视视频播放| 中文字幕熟女人妻在线| 久久精品人妻少妇| 少妇熟女aⅴ在线视频| 日本黄色视频三级网站网址| 在线天堂最新版资源| 又粗又爽又猛毛片免费看| 一a级毛片在线观看| 国产黄a三级三级三级人| 一个人观看的视频www高清免费观看| 日韩亚洲欧美综合| 老司机福利观看| 亚洲三级黄色毛片| .国产精品久久| av在线观看视频网站免费| 国产精品,欧美在线| 我要搜黄色片| 国产片特级美女逼逼视频| 午夜福利成人在线免费观看| 看黄色毛片网站| av.在线天堂| 男人舔女人下体高潮全视频| 狂野欧美激情性xxxx在线观看| 男插女下体视频免费在线播放| 丝袜喷水一区| 1000部很黄的大片| 一夜夜www| 蜜桃久久精品国产亚洲av| 日韩精品中文字幕看吧| 国产男靠女视频免费网站| 亚洲自拍偷在线| 国产探花在线观看一区二区| 大香蕉久久网| 又粗又爽又猛毛片免费看| 亚洲18禁久久av| 97热精品久久久久久| 九九爱精品视频在线观看| 亚洲丝袜综合中文字幕| 看免费成人av毛片| 亚洲乱码一区二区免费版| 国产一区亚洲一区在线观看| 丝袜美腿在线中文| 欧美色欧美亚洲另类二区| 国内精品久久久久精免费| 99久久九九国产精品国产免费| 无遮挡黄片免费观看| 国产一区二区亚洲精品在线观看| 中国美女看黄片| 亚洲精品456在线播放app| 精品久久久久久久人妻蜜臀av| 国产精品久久久久久久电影| 国产不卡一卡二| 最近视频中文字幕2019在线8| 男人的好看免费观看在线视频| 成人二区视频| 91狼人影院| 黄片wwwwww| 99riav亚洲国产免费| 啦啦啦观看免费观看视频高清| 男人舔奶头视频| 观看美女的网站| 亚洲人成网站在线播| 在线天堂最新版资源| 男女边吃奶边做爰视频| 国产91av在线免费观看| 久久精品人妻少妇| 国产精品久久视频播放| 亚洲国产日韩欧美精品在线观看| 国产成人91sexporn| 最近在线观看免费完整版| 久久久久久久久中文| 国产精品无大码| 亚洲精品影视一区二区三区av| 亚洲人成网站高清观看| 亚洲激情五月婷婷啪啪| 精品免费久久久久久久清纯| 午夜精品一区二区三区免费看| 日韩av不卡免费在线播放| 一夜夜www| 国产 一区 欧美 日韩| av在线蜜桃| 露出奶头的视频| 精品日产1卡2卡| 小说图片视频综合网站| 偷拍熟女少妇极品色| 国产久久久一区二区三区| 精品一区二区三区人妻视频| 久久6这里有精品| 亚洲乱码一区二区免费版| 最近视频中文字幕2019在线8| 成年av动漫网址| 日本黄大片高清| 成人漫画全彩无遮挡| 久久久久久久久久久丰满| 国产女主播在线喷水免费视频网站 | АⅤ资源中文在线天堂| 午夜福利成人在线免费观看| 中文字幕av成人在线电影| 欧美一区二区亚洲| 老师上课跳d突然被开到最大视频| 91久久精品电影网| 日日摸夜夜添夜夜爱| 亚洲在线观看片| 国产av在哪里看| 久久午夜福利片| 99精品在免费线老司机午夜| 男女做爰动态图高潮gif福利片| 午夜精品在线福利| 伦理电影大哥的女人| 哪里可以看免费的av片| 亚洲人成网站在线播| 人妻夜夜爽99麻豆av| 激情 狠狠 欧美| 国产黄色视频一区二区在线观看 | or卡值多少钱| 高清毛片免费观看视频网站| 久久草成人影院| 国产精品亚洲美女久久久| 2021天堂中文幕一二区在线观| or卡值多少钱| 国产中年淑女户外野战色| 小蜜桃在线观看免费完整版高清| 国产一级毛片七仙女欲春2| 丰满人妻一区二区三区视频av| videossex国产| 亚洲第一电影网av| 国语自产精品视频在线第100页| 成人国产麻豆网| 国产精品久久久久久av不卡| 亚洲久久久久久中文字幕| 免费无遮挡裸体视频| 久久精品久久久久久噜噜老黄 | 亚洲内射少妇av| 亚洲精品456在线播放app| 国产午夜福利久久久久久| 一夜夜www| 欧美一级a爱片免费观看看| 亚洲精品日韩av片在线观看| 可以在线观看的亚洲视频| 白带黄色成豆腐渣| 国产一区二区三区在线臀色熟女| 午夜精品一区二区三区免费看| 一级a爱片免费观看的视频| 免费观看的影片在线观看| 国产午夜精品论理片| 久久中文看片网| 成年免费大片在线观看| 三级国产精品欧美在线观看| 国产精品久久视频播放| 少妇熟女aⅴ在线视频| 变态另类丝袜制服| 欧洲精品卡2卡3卡4卡5卡区| 成人三级黄色视频| av女优亚洲男人天堂| 看十八女毛片水多多多| 一级毛片久久久久久久久女| 天堂av国产一区二区熟女人妻| 日本色播在线视频| 久久人人爽人人片av| 中文亚洲av片在线观看爽| 综合色丁香网| 国产亚洲精品综合一区在线观看| 色吧在线观看| 亚洲人成网站在线播| 亚洲激情五月婷婷啪啪| 99在线视频只有这里精品首页| 成年版毛片免费区| 最近中文字幕高清免费大全6| 中文字幕久久专区| 搡女人真爽免费视频火全软件 | 国产亚洲精品久久久久久毛片| 欧美zozozo另类| 男女边吃奶边做爰视频| 日本色播在线视频| 色噜噜av男人的天堂激情| 国产一区二区三区在线臀色熟女| 日韩 亚洲 欧美在线| 国产乱人偷精品视频| 久久久久国产精品人妻aⅴ院| 五月伊人婷婷丁香| 91在线精品国自产拍蜜月| 久久精品国产99精品国产亚洲性色| 波多野结衣巨乳人妻| 麻豆久久精品国产亚洲av| 日韩精品有码人妻一区| 一夜夜www| 美女 人体艺术 gogo| 免费av观看视频| 久久人人精品亚洲av| 日韩av在线大香蕉| 又爽又黄a免费视频| 精品熟女少妇av免费看| 有码 亚洲区| 一级黄片播放器| 国产免费男女视频| 热99在线观看视频| 国产成人影院久久av| 国产欧美日韩精品一区二区| 午夜影院日韩av| 久久6这里有精品| 特大巨黑吊av在线直播| 午夜福利高清视频| 欧美成人a在线观看| 欧美人与善性xxx| 一级毛片电影观看 | 我要看日韩黄色一级片| 亚洲av中文字字幕乱码综合| 国产黄a三级三级三级人| 亚洲天堂国产精品一区在线| 色综合亚洲欧美另类图片| 免费看av在线观看网站| 亚洲激情五月婷婷啪啪| 欧美xxxx性猛交bbbb| 精品一区二区三区人妻视频| 国产精品久久视频播放| av天堂中文字幕网| 毛片一级片免费看久久久久| 国内精品一区二区在线观看| 18禁在线播放成人免费| 日日摸夜夜添夜夜添av毛片| 狂野欧美激情性xxxx在线观看| 亚洲三级黄色毛片| АⅤ资源中文在线天堂| 久久久久久久亚洲中文字幕| 亚洲欧美成人综合另类久久久 | 在线免费观看不下载黄p国产| 亚洲成人av在线免费| 欧洲精品卡2卡3卡4卡5卡区| 精品人妻视频免费看| 亚洲欧美中文字幕日韩二区| 永久网站在线| 精品欧美国产一区二区三| 变态另类成人亚洲欧美熟女| 97热精品久久久久久| 亚洲av五月六月丁香网| 国产高清有码在线观看视频| 亚洲最大成人中文| 两个人视频免费观看高清| 欧美日韩国产亚洲二区| 亚洲性夜色夜夜综合| 啦啦啦韩国在线观看视频| 丰满人妻一区二区三区视频av| 亚洲美女搞黄在线观看 | 成人无遮挡网站| 国产免费男女视频| 99九九线精品视频在线观看视频| 久久亚洲国产成人精品v| 色哟哟·www| 色视频www国产| 日韩在线高清观看一区二区三区| 久久久久久久久久黄片| 老司机午夜福利在线观看视频| 国产一级毛片七仙女欲春2| 性色avwww在线观看| h日本视频在线播放| 日日干狠狠操夜夜爽| 国产精品99久久久久久久久| 欧美日韩乱码在线| 成人高潮视频无遮挡免费网站| 亚洲av一区综合| 亚洲成av人片在线播放无| 久久久久久国产a免费观看| 人妻丰满熟妇av一区二区三区| 国产精品,欧美在线| 亚洲第一电影网av| 免费电影在线观看免费观看| 亚洲人成网站高清观看| 国产亚洲91精品色在线| av专区在线播放| 天堂√8在线中文| 天堂网av新在线| 国产一区亚洲一区在线观看| 精品久久久噜噜| 久久久久久久久久黄片| 青春草视频在线免费观看| 亚洲国产欧洲综合997久久,| 大又大粗又爽又黄少妇毛片口| 五月玫瑰六月丁香| 99热精品在线国产| 欧美性猛交黑人性爽| 免费一级毛片在线播放高清视频| 久久久精品94久久精品| 日韩在线高清观看一区二区三区| 亚洲性夜色夜夜综合| 欧美日本视频| 嫩草影院新地址| 日本熟妇午夜| av在线播放精品| 天天一区二区日本电影三级| 色播亚洲综合网| 亚洲av五月六月丁香网| 在线观看午夜福利视频| 亚洲第一区二区三区不卡| 给我免费播放毛片高清在线观看| 国产女主播在线喷水免费视频网站 | 日本熟妇午夜| 男女视频在线观看网站免费| 99视频精品全部免费 在线| 久久天躁狠狠躁夜夜2o2o| .国产精品久久| 97热精品久久久久久| 国产精品亚洲一级av第二区| 性色avwww在线观看| 亚洲中文字幕日韩| 最近在线观看免费完整版| 人妻夜夜爽99麻豆av| 麻豆成人午夜福利视频| 亚洲精品日韩在线中文字幕 | 搡老熟女国产l中国老女人| 美女xxoo啪啪120秒动态图| 久久亚洲国产成人精品v| 国产精品女同一区二区软件| 日韩精品中文字幕看吧| 成人无遮挡网站| 夜夜看夜夜爽夜夜摸| 最近手机中文字幕大全| 12—13女人毛片做爰片一| 亚洲成人久久爱视频| 国产三级中文精品| 一级av片app| www日本黄色视频网| 成年av动漫网址| 日韩欧美一区二区三区在线观看| 久久人妻av系列| 亚洲国产欧洲综合997久久,| 成人精品一区二区免费| 日韩欧美免费精品| 三级经典国产精品| 99久久精品国产国产毛片| 三级国产精品欧美在线观看| 国产中年淑女户外野战色| 免费看a级黄色片| 久久草成人影院| 午夜a级毛片| 国产乱人偷精品视频| 精品乱码久久久久久99久播| 欧美潮喷喷水| 国产精品女同一区二区软件| 免费无遮挡裸体视频| 成人无遮挡网站| 国产极品精品免费视频能看的| 国产淫片久久久久久久久| 在线观看av片永久免费下载| 麻豆精品久久久久久蜜桃| 国产真实伦视频高清在线观看| 精品日产1卡2卡| 国产蜜桃级精品一区二区三区| av在线老鸭窝| 丰满乱子伦码专区| 一个人看视频在线观看www免费| 日本爱情动作片www.在线观看 | 国内精品宾馆在线| 日本熟妇午夜| 国产 一区精品| 男女啪啪激烈高潮av片| 嫩草影院精品99| 尤物成人国产欧美一区二区三区| 国产精品福利在线免费观看| 精品人妻熟女av久视频| 亚洲经典国产精华液单| 欧美性猛交╳xxx乱大交人| 亚洲欧美日韩高清在线视频| 久久精品久久久久久噜噜老黄 | 色哟哟哟哟哟哟| 亚洲七黄色美女视频| 亚洲三级黄色毛片| 日韩制服骚丝袜av| 啦啦啦啦在线视频资源| 22中文网久久字幕| 男人的好看免费观看在线视频| 99久久九九国产精品国产免费| 欧美性猛交╳xxx乱大交人| 国产激情偷乱视频一区二区| 日韩三级伦理在线观看| 欧美国产日韩亚洲一区| 亚洲国产精品sss在线观看| 久久人人爽人人片av| 久久热精品热| 99久久精品国产国产毛片| av天堂中文字幕网| 成人午夜高清在线视频| 美女大奶头视频| 有码 亚洲区| 99热这里只有是精品50| 婷婷六月久久综合丁香| 51国产日韩欧美| 国产成人一区二区在线| 九九久久精品国产亚洲av麻豆| 九九热线精品视视频播放| 小蜜桃在线观看免费完整版高清| 岛国在线免费视频观看| 国产午夜福利久久久久久| 国产中年淑女户外野战色| 3wmmmm亚洲av在线观看| 久久久精品大字幕| 在线观看美女被高潮喷水网站| 老熟妇乱子伦视频在线观看| 亚洲欧美清纯卡通| 亚洲成人久久性| 夜夜夜夜夜久久久久| 久久精品国产亚洲av香蕉五月| 成人综合一区亚洲| 欧美日韩乱码在线| 成人特级黄色片久久久久久久| 人妻少妇偷人精品九色| 一级a爱片免费观看的视频| 国产aⅴ精品一区二区三区波| 变态另类丝袜制服| 少妇熟女aⅴ在线视频| 男插女下体视频免费在线播放| 97碰自拍视频| 亚洲精品乱码久久久v下载方式| 亚洲av成人av| 国产伦一二天堂av在线观看| 一边摸一边抽搐一进一小说| 男人舔奶头视频| 国产91av在线免费观看| 99在线人妻在线中文字幕| 色5月婷婷丁香| 久久午夜福利片| 天天躁夜夜躁狠狠久久av| 欧美成人免费av一区二区三区| av免费在线看不卡| 97超碰精品成人国产| 欧美在线一区亚洲| 成人毛片a级毛片在线播放| 最近在线观看免费完整版| 亚洲第一区二区三区不卡| 国产真实乱freesex| 国产乱人视频| 干丝袜人妻中文字幕| 色噜噜av男人的天堂激情| 国产精品永久免费网站| 免费不卡的大黄色大毛片视频在线观看 | 精品人妻一区二区三区麻豆 | 可以在线观看的亚洲视频| 国产免费一级a男人的天堂| 久久午夜亚洲精品久久| 熟女电影av网| 插阴视频在线观看视频| 网址你懂的国产日韩在线| 亚洲aⅴ乱码一区二区在线播放| 熟妇人妻久久中文字幕3abv| a级毛片a级免费在线| 久久久欧美国产精品| or卡值多少钱| 亚洲精华国产精华液的使用体验 | 99久久无色码亚洲精品果冻| 国产精品一区二区三区四区免费观看 | 99热这里只有是精品在线观看| 人人妻人人澡欧美一区二区| 国产单亲对白刺激| 日本精品一区二区三区蜜桃| 日韩国内少妇激情av| 在现免费观看毛片| 日韩中字成人| 国产精品野战在线观看| 一个人观看的视频www高清免费观看| ponron亚洲| 欧美高清性xxxxhd video| 国产精品日韩av在线免费观看| 看十八女毛片水多多多| 美女 人体艺术 gogo| 色综合亚洲欧美另类图片| 国产精品人妻久久久久久| 亚洲高清免费不卡视频| 国产精品久久久久久亚洲av鲁大| av黄色大香蕉| 偷拍熟女少妇极品色| 插逼视频在线观看| 麻豆精品久久久久久蜜桃| 免费电影在线观看免费观看| 国产精品99久久久久久久久| 精品欧美国产一区二区三| 日韩中字成人| 真人做人爱边吃奶动态| 中国美女看黄片| 2021天堂中文幕一二区在线观| 十八禁国产超污无遮挡网站| 搡老妇女老女人老熟妇| 成人毛片a级毛片在线播放| 内地一区二区视频在线| 97在线视频观看| 一级黄色大片毛片| 最新在线观看一区二区三区| 国产高清视频在线观看网站| 成人高潮视频无遮挡免费网站| 午夜爱爱视频在线播放| 菩萨蛮人人尽说江南好唐韦庄 | 91精品国产九色| 亚洲内射少妇av| 国产精品1区2区在线观看.| aaaaa片日本免费| 在线免费观看的www视频| 伊人久久精品亚洲午夜| 国产真实乱freesex| 久久婷婷人人爽人人干人人爱| 在线播放无遮挡| 亚洲四区av| avwww免费| 可以在线观看毛片的网站| 在线观看美女被高潮喷水网站| 国产精品1区2区在线观看.| 欧美zozozo另类| 欧美成人免费av一区二区三区| 在线天堂最新版资源| 精品无人区乱码1区二区| 亚洲成人中文字幕在线播放| 色尼玛亚洲综合影院| 国产精品人妻久久久影院| 久久久欧美国产精品| 天美传媒精品一区二区| 色综合色国产| 禁无遮挡网站| 国产爱豆传媒在线观看| 国产精品伦人一区二区| 色综合亚洲欧美另类图片| 午夜精品一区二区三区免费看| 久久久久久大精品| 亚洲av五月六月丁香网| 国产真实乱freesex| 五月伊人婷婷丁香| 啦啦啦观看免费观看视频高清| 俄罗斯特黄特色一大片| 国产白丝娇喘喷水9色精品| 在现免费观看毛片| 国产av在哪里看| 欧美国产日韩亚洲一区| 国产私拍福利视频在线观看| 国产色婷婷99| 国产av不卡久久| 国产精品不卡视频一区二区| 亚洲aⅴ乱码一区二区在线播放| 国产精品亚洲一级av第二区| 可以在线观看的亚洲视频| 日本色播在线视频| 丝袜喷水一区|