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

    仿舵型量水槽體形模擬研究

    2022-06-07 10:36:42張新燕劉朝陽(yáng)
    灌溉排水學(xué)報(bào) 2022年5期
    關(guān)鍵詞:測(cè)流線形水頭

    李 偉,張新燕,2*,劉朝陽(yáng)

    (1.西北農(nóng)林科技大學(xué) 水利與建筑工程學(xué)院,陜西 楊凌712100;2.西北農(nóng)林科技大學(xué) 旱區(qū)農(nóng)業(yè)水土工程教育部重點(diǎn)試驗(yàn)室,陜西 楊凌 712100)

    仿舵型量水槽體形模擬研究

    李 偉1,張新燕1,2*,劉朝陽(yáng)1

    (1.西北農(nóng)林科技大學(xué) 水利與建筑工程學(xué)院,陜西 楊凌712100;2.西北農(nóng)林科技大學(xué) 旱區(qū)農(nóng)業(yè)水土工程教育部重點(diǎn)試驗(yàn)室,陜西 楊凌 712100)

    【目的】探索具有水流平順、流態(tài)穩(wěn)定、測(cè)流準(zhǔn)確、水頭損失小等優(yōu)良特性的量水槽線形?!痉椒ā酷槍?duì)基于水下裝備流線型體構(gòu)造的量水槽,選取4 種量水槽體形,運(yùn)用數(shù)值模擬方法研究不同體形的量水槽在不同流量條件下的流速分布、水面線、佛汝德數(shù)、水頭損失、壅水高度等水力性能?!窘Y(jié)果】流線型體擴(kuò)散段較陡時(shí),仿舵型量水槽的壅水高度和水頭損失相對(duì)較小。不同流量條件下,4 種體形的量水槽上游佛汝德數(shù)均低于0.5;壅水高度均低于3 cm;計(jì)算流量與模擬流量的平均誤差為2.42%,測(cè)流精度較高,且均具有較小的水頭損失。【結(jié)論】卡克斯仿舵型量水槽在各流量條件下的水面線均較為平穩(wěn)且流速分布更加均勻,在不同流量條件下其水頭損失均小于0.35 cm,測(cè)流精度大于98.27%,水力性能最優(yōu),可作為仿舵型量水槽的基本體形。

    矩形渠道;仿舵型量水槽;數(shù)值模擬;流線型體

    0 引 言

    【研究意義】灌區(qū)量水是灌區(qū)管理部門進(jìn)行正確引水、輸水和水量調(diào)配的主要手段,是發(fā)展節(jié)水型農(nóng)業(yè)和實(shí)現(xiàn)農(nóng)業(yè)高效用水的關(guān)鍵技術(shù)。因此,探索結(jié)構(gòu)簡(jiǎn)單、水力性能優(yōu)良、量水精度高、滿足灌區(qū)需求的渠道量水設(shè)施具有重要意義[1-2]。堰槽類量水設(shè)施具有經(jīng)濟(jì)實(shí)用且穩(wěn)定可靠的特點(diǎn),適合在灌區(qū)大規(guī)模推廣應(yīng)用?!狙芯窟M(jìn)展】在過(guò)去,量水設(shè)施的體型由初期的巴歇爾量水槽、無(wú)喉道量水槽、直壁式量水槽等易于產(chǎn)生邊界層分離、水力特性復(fù)雜、量測(cè)難度較大的折線形體型發(fā)展到了拋物線形量水槽[3]、半圓柱形量水槽[4]、機(jī)翼形量水槽[5-6]等水流平順、流態(tài)穩(wěn)定、測(cè)量精度高的曲線形體型。近年來(lái),隨著數(shù)值模擬技術(shù)的發(fā)展,CFD 技術(shù)被廣泛應(yīng)用于量水槽體型設(shè)計(jì)。許虎[7]以巴歇爾量水槽為研究對(duì)象,模擬分析了量水槽進(jìn)口連接段形式對(duì)水面線、水頭損失、測(cè)流精度、流速及壓強(qiáng)等水力性能的影響,確定出水力性能最優(yōu)的進(jìn)口段形式。宋金妍[8]通過(guò)對(duì)經(jīng)典翼形量水槽與優(yōu)化翼形量水槽2 種體型的量水槽的水力特性進(jìn)行對(duì)比分析,得出優(yōu)化翼形量水槽的流速分布、水面線、佛汝德數(shù)、測(cè)流精度等水力性能較優(yōu)。楊珮珮[9]通過(guò)對(duì)適用于弧底梯形渠道的新型無(wú)喉道量水槽在不同收縮比和喉口邊坡系數(shù)條件下的水流特性研究,分析了其收縮比和喉口邊坡系數(shù)的最佳取值范圍,確定了量水槽的最優(yōu)體型參數(shù)。以往研究均在現(xiàn)有量水設(shè)施的基礎(chǔ)上,通過(guò)改變量水槽局部形狀(連接段)或參數(shù)(邊坡系數(shù))進(jìn)行量水槽體型的設(shè)計(jì)優(yōu)化?!厩腥朦c(diǎn)】而結(jié)合灌區(qū)實(shí)際需求,根據(jù)渠道水力特性,研究量測(cè)精度高、適應(yīng)性好、水力性能優(yōu)良的新型量水設(shè)備——仿舵型量水槽的研究較少?!緮M解決的關(guān)鍵問(wèn)題】為此,本研究基于潛水設(shè)備外形構(gòu)造,通過(guò)幾何組合法、源匯法、線形擬合等方法設(shè)計(jì)并擬合4 條流線型構(gòu)造的量水槽體型,基于FLOW-3D 軟件對(duì)量水槽水力性能進(jìn)行數(shù)值模擬研究,尋求仿舵型量水槽的最優(yōu)體型,以期為推動(dòng)智慧型現(xiàn)代化灌區(qū)建設(shè)和精細(xì)化管理提供技術(shù)支撐。

    1 模型建立

    1.1 量水槽結(jié)構(gòu)與流線型

    仿舵型量水槽的基本結(jié)構(gòu)由上游收縮段、下游擴(kuò)散段和喉口段3 部分組成,如圖1 所示。基于潛水設(shè)備的外形構(gòu)造選取的4 種流線型如圖2 所示,分別為:半橢圓+半橢圓線型(EE 型)、半橢圓+拋物線線型(EP 型)、卡克斯線型(KKS 型)、CB/Z-237 引申的離散線型(CBD 型)[10]。4 種線型均由數(shù)學(xué)方程控制,其中,EE 型與EP 型的線形(圖2(a))收縮段相同,但擴(kuò)散段較陡,曲率半徑大;KKS 型與CBD型的線形(圖2(b))收縮段及擴(kuò)散段均較陡,曲率半徑均較大。

    圖1 仿舵型量水槽結(jié)構(gòu)Fig.1 Structure diagram of rudder-like flume

    圖2 流線型線形Fig.2 Streamline shapes

    1.2 數(shù)學(xué)模型

    1.2.1 控制方程和湍流模型

    將過(guò)槽水流視為不可壓縮的三維黏性液體,根據(jù)水流運(yùn)動(dòng)特性,采用RNG k-ε 紊流模型進(jìn)行模擬研究,控制方程如下:

    連續(xù)性方程:

    動(dòng)量方程:

    湍動(dòng)能k方程:

    湍動(dòng)能耗散率ε方程:

    式中:t為時(shí)間(s);ρ為水密度(kg/m3),取值為1 000;ui、uj分別為流速矢量在xi、xj方向上的分量(i,j=1,2,3)(m/s);μ為流體動(dòng)力黏滯系數(shù)(kg/(m·s)),取值為1.005×10-3((N·s)/m2);k為湍動(dòng)能(m2/s2);ε為湍動(dòng)耗散率((kg·m2)/s3);μeff為流體的有效動(dòng)力黏滯系數(shù);Gk為平均流速梯度引起的湍動(dòng)能k的產(chǎn)生項(xiàng)。

    1.2.2 邊界條件及網(wǎng)格劃分

    本研究建立的量水槽三維立體模型如圖3 所示。仿舵型量水槽模型的槽長(zhǎng)為0.6 m,槽厚為0.12 m,收縮比為0.6。模型放置在長(zhǎng)12 m、寬0.6 m、高0.3 m 的矩形平底渠道上,量水槽的槽首距離渠道進(jìn)口4 m。渠道進(jìn)口邊界設(shè)定為流速進(jìn)口,設(shè)定流速及水位以改變流量;出口設(shè)為自由出流;渠道邊壁選擇無(wú)滑移固體邊壁;自由水面與空氣交界面設(shè)為對(duì)稱邊界。

    圖3 仿舵型量水槽三維模型及邊界設(shè)置Fig.3 Three-dimensional model and boundary setting of rudder-like flume

    通過(guò)多次網(wǎng)格劃分并以相同初始條件計(jì)算出口流量進(jìn)行網(wǎng)格密度分析,結(jié)果如圖4 所示??紤]在保證精度的前提下節(jié)約計(jì)算資源,最終取網(wǎng)格數(shù)為357 000,即計(jì)算區(qū)域單元網(wǎng)格長(zhǎng)度整體設(shè)置為0.02 m。由于量水槽曲線較為復(fù)雜且附近流場(chǎng)變化劇烈,因此將量水槽處的網(wǎng)格進(jìn)行加密,量水槽局部加密網(wǎng)格劃分如圖5 所示,設(shè)置為0.01 m,最終整個(gè)計(jì)算區(qū)域網(wǎng)格總數(shù)為504 000。

    圖4 網(wǎng)格無(wú)關(guān)性檢驗(yàn)Fig.4 Grid independence check

    圖5 局部加密網(wǎng)格劃分Fig.5 Local grid division diagram

    2 模擬結(jié)果分析

    2.1 水流流態(tài)

    圖6 Fr 沿程變化Fig.6 The longitudinal variation of Froude number

    圖6 為不同體形量水槽Fr沿程變化(過(guò)流量Q=44.92 L/s)。不同線形的量水槽過(guò)槽水流變化趨勢(shì)基本一致,槽前水流由于受到量水槽的阻擋,水位壅高,流速減小,佛汝德數(shù)Frlt;0.5。4 種線型的量水槽進(jìn)水前端線形圓滑,進(jìn)槽水流均無(wú)邊界層分離現(xiàn)象,水流順暢,緊貼槽壁,流態(tài)平穩(wěn)。在收縮段由于束窄作用,水流流速增大,水深減小,發(fā)生水跌過(guò)渡到下游急流,并在喉口斷面下游附近發(fā)生臨界流,在擴(kuò)散段末端附近形成水深最小、流速最大的收縮斷面,隨后發(fā)生水躍與下游渠道緩流相接。

    槽前Fr變化較小且未超過(guò)0.4,從槽前(x=3.6 m)到槽進(jìn)口(x=4 m)處Fr開(kāi)始逐漸增大,在槽進(jìn)口處達(dá)到0.43,仍小于0.5,即各線型量水槽均滿足測(cè)流規(guī)范中的上游佛汝德數(shù)不大于0.5 的基本要求。在槽收縮段斷面收縮,水面跌落,水深減小,F(xiàn)r增大,在喉口下游附近形成臨界流,在出口處形成收縮斷面,F(xiàn)r達(dá)到最大。不同體形量水槽收縮斷面位置略有不同,其Fr最大值也有所差異。Fr最大值由大到小依次為EE 型、KKS 型、EP 型和CBD 型,其中EE 型Fr最大,為1.22;CBD 型Fr最小,為1.12。4 種線型量水槽躍前斷面均處于槽擴(kuò)散段,其水躍發(fā)生位置及躍前斷面Fr與線形之間存在聯(lián)系。EE 型和KKS型擴(kuò)散段更陡,其水躍發(fā)生位置離喉口更遠(yuǎn),且其躍前斷面Fr更大。由此可見(jiàn),流線型體量水槽擴(kuò)散段越陡時(shí),其水躍越接近擴(kuò)散段下游發(fā)生,躍前斷面Fr越大。

    2.2 水面線

    由于形體差異,不同流量條件下量水槽水面線變化規(guī)律基本相似,但相同流量條件下各量水槽之間仍存 在差異。以流量Q=44.92 L/s 時(shí)量水槽水面線變化為例(圖7)??梢钥闯?,水面線整體變化趨勢(shì)一致,但收縮水深大小有所差異。量水槽上游水流平緩,在槽前水深稍有降低,受量水槽收縮影響,水面在槽進(jìn)口(x=4 m)開(kāi)始迅速降落,水深快速減小。喉口斷面面積最小,理論上在該斷面會(huì)形成臨界流,但由于受邊界影響,臨界流斷面發(fā)生在喉口下游某處,此時(shí)水流由緩流過(guò)渡到臨界流,之后水流受擴(kuò)散段作用繼續(xù)跌落,在距槽進(jìn)口4.5 m 附近形成水深最小、流速最大的收縮斷面,之后發(fā)生水躍與下游緩流相銜接。

    圖7 水面線變化Fig.7 Variation of water surface line

    仿舵型量水槽特征水深及斷面位置如表1 所示。4 種體型量水槽臨界流斷面發(fā)生在喉口下游3~5 cm處,其中CBD 型臨界水深值最大,為12.11 cm,而KKS 型臨界流斷面發(fā)生位置距喉口斷面最近,可提高臨界流斷面位置的確定及臨界水深的測(cè)量精度。各種體型量水槽收縮斷面均發(fā)生在量水槽出口附近,KKS 型收縮水深最小,為7.74 cm,最大為EP 型,8.35 cm。由此可以得出,流線型體量水槽擴(kuò)散段越陡,收縮斷面位置距槽喉口越遠(yuǎn),收縮水深越小。

    表1 仿舵型量水槽特征水深及斷面位置Table 1 Characteristic depths of rudder-like flume

    2.3 流速分布

    由圖8 可知,4 種線型的量水槽上游流速基本相同。流速在槽前(x=3.6 m)處開(kāi)始快速增大,在收縮斷面達(dá)到最大值,流速最大值由大到小依次為KKS型、EE 型、CBD 型和EP 型。KKS 型流速最大,1.24 m/s,CBD 型最小,為1.18 m/s。經(jīng)過(guò)收縮斷面之后流速迅速減小,在槽后逐漸穩(wěn)定。從線形來(lái)看,擴(kuò)散段較陡的EE 型和KKS 型流速上限較高。即流線型體量水槽擴(kuò)散段越陡時(shí),流速最大值越大。

    水頭損失是衡量量水設(shè)施優(yōu)劣的重要技術(shù)指標(biāo)之一。明渠流水頭損失由沿程水頭損失和局部水頭損失組成。由于量水槽內(nèi)水頭損失主要為水躍產(chǎn)生的能量損失,故只計(jì)算水躍產(chǎn)生的局部水頭損失。表2 為量水槽在流量Q=44.92 L/s 時(shí)的水躍特征參數(shù)及能量損失。4 種仿舵型量水槽的躍前斷面均位于量水槽擴(kuò)散段,躍前水深略有差異,最大為EP 型,其次為EE型;KKS 型躍長(zhǎng)最小,CBD 型躍長(zhǎng)最大;4 種線形量水槽的最大躍高相差不大。能量損失最少的是KKS型量水槽,最大的是EP 型量水槽。

    圖8 中軸線縱向流速沿程變化Fig.8 Variation of longitudinal velocity along axis

    2.4 水頭損失

    表2 仿舵型量水槽水躍特征參數(shù)及能量損失Table 2 Hydraulic jump energy loss of rudder-like flume

    2.5 壅水高度

    水流流經(jīng)渠道砌筑量水槽時(shí)必然會(huì)形成收縮,從而抬高上游水深,造成一定程度壅水。壅水高度為某一流量水流流經(jīng)量水槽時(shí)相對(duì)應(yīng)的上游水深與未修筑量水槽渠道通過(guò)同一流量水流的正常水深差值,壅水高度過(guò)大,水流溢出,造成流量損失及量測(cè)誤差。4 種體形量水槽在Q=44.92 L/s 時(shí)的上、下游水深及壅水高度如表3 所示。可以看出,量水槽的壅水高度均不大,其中KKS 型的壅水高度最小,為2.78 cm,之后由小到大依次為EE 型、CBD 型、EP 型。不同流量條件下壅水高度有所不同,但均小于3 cm,滿足規(guī)范標(biāo)準(zhǔn)要求。從體形上看,EE 型擴(kuò)散段較EP 型陡,KKS 型較CBD 型擴(kuò)散段也更陡,而KKS 型和EE 型的壅水高度都更小,即流線型體量水槽擴(kuò)散段越陡時(shí),壅水高度越小。

    表3 仿舵型量水槽壅水高度Table 3 Backwater height of rudder-like flume

    2.6 測(cè)流精度

    量水槽的測(cè)流精度是衡量其性能優(yōu)劣的重要指標(biāo)之一。表4 為計(jì)算流量與模擬流量的誤差??梢钥闯?,計(jì)算流量與模擬流量的最大相對(duì)誤差為3.76%,最小相對(duì)誤差僅為0.69%,4 種體形量水槽的平均誤差為2.42%,其中KKS 型量水槽在各流量條件下的平均誤差最小,為1.73%,測(cè)流精度最高。

    表4 仿舵型量水槽計(jì)算流量與模擬流量誤差Table 4 Comparison between calculated flow and simulated flow of rudder-like flume

    3 討 論

    近年來(lái),隨著灌區(qū)量水技術(shù)的不斷發(fā)展,量水槽形式由折線形式逐漸轉(zhuǎn)變?yōu)榍€形式,主要原因是折線形容易產(chǎn)生流線與邊壁分離現(xiàn)象,且水力特性復(fù)雜,量測(cè)難度較大,而曲線形式的外部形狀符合水流流線分布、邊界層分離現(xiàn)象弱,水流阻力小。本文采用水下設(shè)備的設(shè)計(jì)思路與方法,利用灌區(qū)渠道水流水力特性去探求最佳流線型體以應(yīng)用于新型的測(cè)量精準(zhǔn)、適用度高、應(yīng)用方便、性能優(yōu)良的渠道量水設(shè)備——仿舵型量水槽。

    本研究中,卡克斯型量水槽在各流量條件下的水面線均較為平穩(wěn)且流速分布更加均勻,水力性能較好,因此建議將其作為仿舵型量水槽基本形體。這與肖華攀等[11]認(rèn)為半橢圓+圓弧流線型的阻力特性最好的研究結(jié)果存在差異。原因在于后者為了探求流線型體結(jié)構(gòu)的設(shè)計(jì)方法,在幾種常用流線型體輪廓方程基礎(chǔ)上,對(duì)比分析不同結(jié)構(gòu)參數(shù)下4 種線形二維流線型體的水下阻力特性研究,而本研究是將4 種線形與矩形渠道結(jié)合構(gòu)成仿舵型量水槽,進(jìn)一步研究其水力性能得出結(jié)論;且其針對(duì)各線形的阻力特性進(jìn)行研究,而本研究綜合考慮流速分布、水面線、佛汝德數(shù)、水頭損失、壅水高度等水力性能;另外由于線形收縮段長(zhǎng)度的變化也會(huì)對(duì)阻力性能產(chǎn)生影響,因此得到的結(jié)論有一定的差異。賀夢(mèng)楊等[12]對(duì)矩形渠道機(jī)翼柱形量水槽進(jìn)行了模型試驗(yàn)和數(shù)值模擬研究,其流量公式計(jì)算誤差為6.34%,高于本研究的平均測(cè)流誤差(2.42%),這是由于其考慮的是圓柱體放置在中間的情況,與本文放置方式的不同導(dǎo)致了測(cè)流精度的差異。劉鴻濤[13]基于Fluent 軟件,對(duì)U 形渠道半圓柱量水槽進(jìn)行了數(shù)值模擬。紊流模型采用標(biāo)準(zhǔn)k-ε模型,VOF 法追蹤自由液面的模擬方法,試驗(yàn)流量與模擬流量平均誤差達(dá)到3.77%,同樣滿足數(shù)值模擬對(duì)精度的要求,但誤差略高于使用本研究誤差,也側(cè)面驗(yàn)證了本方法對(duì)量水槽研究的適用性。

    另外,研究通過(guò)對(duì)上游佛汝德數(shù)、流速分布、水頭損失等水力性能分析后發(fā)現(xiàn)在相同收縮比條件下,量水槽擴(kuò)散段越陡,壅水高度和水頭損失越小,與宋金妍[8]通過(guò)對(duì)量水槽的結(jié)構(gòu)形式、測(cè)流機(jī)理進(jìn)行分析研究將經(jīng)典翼形量水槽與優(yōu)化翼形量水槽的水力特性對(duì)比模擬分析之后的結(jié)果具有一致性,說(shuō)明量水槽線形曲線對(duì)量水槽水力性能影響較大,可以通過(guò)調(diào)整線形來(lái)達(dá)到提高水力性能的目的,也可為量水槽體形的進(jìn)一步優(yōu)化和簡(jiǎn)化提供便利[14-17]。

    4 結(jié) 論

    1)4 種形體量水槽模型在各流量條件下公式計(jì)算流量與模擬流量平均誤差為2.42%,測(cè)量精度較高。上游佛汝德數(shù)均小于0.5,上游易形成平穩(wěn)緩流。壅水高度和水頭損失均較小,滿足既成渠道的應(yīng)用要求。

    2)在相同收縮比條件下,量水槽擴(kuò)散段越陡,壅水高度和水頭損失越小。

    3)KKS 型在試驗(yàn)條件下,流速分布更為均勻,流態(tài)更為穩(wěn)定,測(cè)流精度最高;其模型計(jì)算流量與模擬流量平均誤差僅為1.73%,水頭損失小于0.35 cm,壅水高度小于3 cm;水力性能最優(yōu),可作為仿舵型量水槽的基本體型。

    [1]陳志愷. 中國(guó)水資源的可持續(xù)利用問(wèn)題[J]. 水文, 2003 (1): 1-5.

    CHEN Zhikai. Sustainable utilization of water resources in China[J].Hydrology, 2003(1): 1-5.

    [2]康紹忠. 加快推進(jìn)灌區(qū)現(xiàn)代化改造,補(bǔ)齊國(guó)家糧食安全短板[J]. 中國(guó)水利, 2020(9): 1-5.

    KANG Shaozhong. Accelerate the modernization of irrigation districts to complement national food security short board[J]. China Water Conservancy, 2020(9): 1-5.

    [3]呂宏興, 朱曉群, 張春娟, 等. U 形渠道拋物線形喉口式量水槽選形與設(shè)計(jì)[J]. 灌溉排水, 2001, 20(2): 55-57.

    LYU Hongxing, ZHU Xiaoqun, ZHANG Chunjuan, et al. Selection and design of U-shaped channel parabolic throat flume[J]. Irrigation and Drainage, 2001, 20 (2): 55-57.

    [4]蔣沛, 陳新明. U 形渠道半圓柱形量水槽的數(shù)值模擬[J]. 灌溉排水學(xué)報(bào), 2016, 35(3): 68-72.

    JIANG Pei, CHEN Xinming. Numerical simulation of semi-cylindrical U-shaped channel flume[J]. Journal of Irrigation and Drainage, 2016,35(3): 68-72.

    [5]呂宏興, 劉煥芳, 朱曉群, 等. 機(jī)翼形量水槽的試驗(yàn)研究[J]. 農(nóng)業(yè)工程學(xué)報(bào), 2006, 22(9): 119-123.

    LYU Hongxing, LIU Huanfang, ZHU Xiaoqun, et al. Experimental study on airfoil-shaped flume[J]. Agricultural Engineering Journal, 2006,22(9): 119-123.

    [6]孫斌. 矩形渠道機(jī)翼形量水槽水力特性數(shù)值模擬[D]. 楊凌: 西北農(nóng)林科技大學(xué), 2010.

    SUN Bin. Numerical simulation of hydraulic characteristics of rectangular channel airfoil-shaped flume[D]. Yangling: Northwest Agricultural and Forestry University, 2010.

    [7]許虎. 量測(cè)一體化巴歇爾槽數(shù)值模擬與結(jié)構(gòu)優(yōu)化研究[D]. 石河子:石河子大學(xué), 2021.

    XU Hu. Numerical simulation and structural optimization of integrated Bachelles[D]. Shihezi: Shihezi University, 2021.

    [8] 宋金妍. 翼型量水槽水力特性數(shù)值模擬及體型類比研究[D]. 鄭州:鄭州大學(xué), 2017.

    SONG Jinyan. Numerical simulation of hydraulic characteristics and shape analogy of airfoil flume[D]. Zhengzhou: Zhengzhou University,2017.

    [9] 楊珮珮. 弧底梯形無(wú)喉道量水槽的數(shù)值模擬與體型優(yōu)化研究[D]. 楊凌: 西北農(nóng)林科技大學(xué), 2016.

    YANG Peipei. Numerical simulation and shape optimization of an arc bottom trapezoidal throatless flume[D]. Yangling: Northwest A amp; F University, 2016.

    [10] 722 研究所. CB/Z 237—87[S]. 北京: 中國(guó)船舶工業(yè)總公司, 1987.722 Research Institute. CB/Z 237—87[S]. Beijing: China State Shipbuilding Industry Corporation, 1987.

    [11] 肖華攀, 方子帆, 朱陳何, 等. 水下流線型體的結(jié)構(gòu)阻力特性及其應(yīng)用研究[J]. 三峽大學(xué)學(xué)報(bào)(自然科學(xué)版), 2013, 35(6): 92-96.

    XIAO Huapan, FANG Zifan, ZHU Chenhe, et al. Structural resistance characteristics of underwater streamlines and their applications[J].Journal of Three Gorges University (Natural Science Edition), 2013, 35(6): 92-96.

    [12] 賀夢(mèng)楊, 尚海鑫, 張寬地, 等. 仿機(jī)翼形便攜式量水槽水力特性試驗(yàn)與數(shù)值模擬[J]. 農(nóng)業(yè)工程學(xué)報(bào), 2021, 37(12): 117-124.

    HE Mengyang, SHANG Haixin, ZHANG Kuandi, et al. Hydraulic characteristics test and numerical simulation of airfoil-like portable flume[J]. Agricultural Engineering Journal, 2021, 37(12): 117-124.

    [13] 劉鴻濤, 于明舟, 牛炎, 等. 梯形渠道翼柱型量水槽試驗(yàn)研究與數(shù)值模擬[J]. 灌溉排水學(xué)報(bào), 2019, 38(9): 101-107.

    LIU Hongtao, YU Mingzhou, NIU Yan, et al. Experimental study and numerical simulation of trapezoidal channel airfoil-column flume[J].Journal of Irrigation and Drainage, 2019, 38(9): 101-107.

    [14] 胡云進(jìn), 萬(wàn)五一, 蔡甫款, 等. 窄深矩形斷面明渠流速分布的研究[J].浙江大學(xué)學(xué)報(bào)(工學(xué)版), 2008(1): 183-187.

    HU Yunjin, WAN Wuyi, CAI Fukuan, et al. Study on velocity distribution of narrow and deep rectangular open channel[J]. Journal of Zhejiang University (Engineering Edition), 2008(1): 183-187.

    [15] 槐文信, 趙明登, 童漢毅. 河道及近海水流的數(shù)值模擬[M]. 北京:科學(xué)出版社, 2005.

    HUAI Wenxin, ZHAO Mingden, TONG Hanyi. Numerical simulation of river and offshore currents[M]. Beijing: Science Press, 2005.

    [16] 冉聃頡, 王文娥, 胡笑濤. 梯形喉口無(wú)喉道量水槽設(shè)計(jì)及其水力性能模擬與試驗(yàn)[J]. 農(nóng)業(yè)工程學(xué)報(bào), 2017, 33(15): 177-183.

    RAN Danjie, WANG Wen’e, HU Xiaotao. Design of trapezoidal throatless flow tank and its hydraulic performance simulation and test[J]. Journal of Agricultural Engineering, 2017, 33(15): 177-183.

    [17] 杜月中, 閔健, 郭字洲. 流線形回轉(zhuǎn)體外形設(shè)計(jì)綜述與線形擬合[J].聲學(xué)技術(shù), 2004, 23(2): 93-102.

    DU Yuezhong, MIN Jian, GUO Zizhou. Overview of shape design and linear fitting of streamlined revolution[J]. Acoustic technology, 2004,23(2): 93-102.

    The Impact of Geometry of Rudder-like Flume on Its Performance and Accuracy

    LI Wei1, ZHANG Xinyan1,2*, LIU Zhaoyang1
    (1. School of Water Conservancy and Construction Engineering, Northwest A amp; F University, Yangling 712100, China;2.Key Laboratory of Agricultural Soil and Water Engineering Ministry of Education Northwest A amp; F University, Yangling 712100, China)

    【Objective】Flume is a hydraulic device widely used in hydraulic engineering projects. Its performance depends on many factors. In this paper, we explore the alignment of measuring flume with smooth flow, stable flow pattern, accurate flow measurement and small head loss, based on the streamline structure of underwater equipment.【Method】 Four flumes with different streamline shapes were studied. Numerical simulation was used to study their hydraulic performance under different flow rates, including velocity distribution, water surface profile, Froude number, water head loss and backwater height.【Result】It is found that the backwater height of the flume and the head loss both decreased when the diffusion became steeper. The upper Froude numbers of the four flumes were less than 0.5, regardless of the flow rates. The backwater heights were less than 3 cm. The average error between the measured and simulated flow rate was less than 2.42%, indicating that the flow measurement is accurate, and the head loss is small.【Conclusion】The water surface profile of the KKS type flume is more stable, and the velocity distribution is more uniform under different flow rates, compared to other flumes. Its head loss is less than 0.35 cm,the flow measurement accuracy is more than 98.27%. Its hydraulic performance is hence the best among the flumes we compared.

    rectangular channel; rudder-like measuring flume; numerical simulation; streamline shape

    S565

    A

    10.13522/j.cnki.ggps.2021551

    李偉, 張新燕, 劉朝陽(yáng). 仿舵型量水槽體形模擬研究[J]. 灌溉排水學(xué)報(bào), 2022, 41(5): 118-123.

    LI Wei, ZHANG Xinyan, LIU Zhaoyang. The Impact of Geometry of Rudder-like Flume on Its Performance and Accuracy[J]. Journal of Irrigation and Drainage, 2022, 41(5): 118-123.

    1672-3317(2022)05-0118-06

    2021-11-10

    國(guó)家重點(diǎn)研發(fā)計(jì)劃項(xiàng)目(2021YFD1900701)

    李偉(1998-),男,陜西延安人。碩士研究生,主要從事水力學(xué)與水工建筑物研究。E-mail: lw18821702791@163.com

    張新燕(1972-),女,河南開(kāi)封人。副教授,博士,主要從事節(jié)水灌溉技術(shù)與設(shè)備研究工作。E-mail: xnzxy@ 163.com

    責(zé)任編輯:韓 洋

    猜你喜歡
    測(cè)流線形水頭
    渠道斷面自動(dòng)測(cè)流系統(tǒng)在位山灌區(qū)測(cè)水量水中的應(yīng)用
    水文測(cè)流技術(shù)方法與進(jìn)展分析
    石河子科技(2022年4期)2022-03-24 05:45:28
    玉龍水電站機(jī)組額定水頭選擇設(shè)計(jì)
    小水電(2021年6期)2021-12-15 02:00:06
    短線法預(yù)制節(jié)段梁線形綜合控制技術(shù)研究
    大跨度連續(xù)剛構(gòu)橋線形控制分析
    泵房排水工程中剩余水頭的分析探討
    曹店灌區(qū)渠首測(cè)流存在的問(wèn)題及對(duì)策
    山東水利(2018年6期)2018-03-24 13:00:35
    M9在建設(shè)在線雷達(dá)測(cè)流設(shè)備選址中的應(yīng)用
    加成固化型烯丙基線形酚醛樹(shù)脂研究
    溪洛渡水電站機(jī)組運(yùn)行水頭處理
    免费观看的影片在线观看| 交换朋友夫妻互换小说| 欧美日韩国产mv在线观看视频| 黄片无遮挡物在线观看| 夫妻性生交免费视频一级片| 免费看光身美女| 成人综合一区亚洲| 精品99又大又爽又粗少妇毛片| av卡一久久| 大香蕉久久成人网| 国产一区有黄有色的免费视频| 青青草视频在线视频观看| 国产视频首页在线观看| 久久人人爽av亚洲精品天堂| 久久久久精品久久久久真实原创| 日韩熟女老妇一区二区性免费视频| 999精品在线视频| 九色亚洲精品在线播放| 少妇猛男粗大的猛烈进出视频| 免费av中文字幕在线| 日韩不卡一区二区三区视频在线| 一区二区三区精品91| 男的添女的下面高潮视频| 十八禁网站网址无遮挡| 国产精品99久久99久久久不卡 | 午夜激情久久久久久久| 精品少妇黑人巨大在线播放| 妹子高潮喷水视频| 日韩成人av中文字幕在线观看| 人成视频在线观看免费观看| .国产精品久久| 免费黄频网站在线观看国产| 成人影院久久| av又黄又爽大尺度在线免费看| 国产av码专区亚洲av| 国产精品人妻久久久久久| 国产一区有黄有色的免费视频| 国产深夜福利视频在线观看| 成人二区视频| 亚洲精品国产av成人精品| 精品少妇久久久久久888优播| 2018国产大陆天天弄谢| 成人国产av品久久久| 亚洲国产精品999| 欧美 亚洲 国产 日韩一| 欧美日韩成人在线一区二区| 精品人妻熟女毛片av久久网站| 免费黄网站久久成人精品| 国产在视频线精品| 搡老乐熟女国产| 少妇人妻久久综合中文| 这个男人来自地球电影免费观看 | 久久99热6这里只有精品| 欧美日韩视频高清一区二区三区二| 亚洲av福利一区| 久久精品人人爽人人爽视色| 一区二区日韩欧美中文字幕 | 久久久久久伊人网av| 亚洲精品,欧美精品| 我的老师免费观看完整版| 欧美精品一区二区大全| 一区二区日韩欧美中文字幕 | 欧美人与性动交α欧美精品济南到 | 男女高潮啪啪啪动态图| 制服诱惑二区| 国产精品国产三级专区第一集| 人成视频在线观看免费观看| 高清黄色对白视频在线免费看| 99视频精品全部免费 在线| 伊人久久精品亚洲午夜| 亚洲av二区三区四区| 国产在线视频一区二区| www.av在线官网国产| 高清视频免费观看一区二区| 91精品三级在线观看| 18+在线观看网站| 婷婷色综合www| 国产成人91sexporn| 你懂的网址亚洲精品在线观看| 欧美精品高潮呻吟av久久| 亚洲高清免费不卡视频| 久久精品国产自在天天线| 国产成人精品福利久久| 国产午夜精品一二区理论片| 亚洲av不卡在线观看| 国内精品宾馆在线| 男人爽女人下面视频在线观看| 一级毛片我不卡| 性色avwww在线观看| 久久韩国三级中文字幕| 亚洲欧美日韩卡通动漫| 亚洲色图综合在线观看| 狂野欧美激情性xxxx在线观看| 超色免费av| 精品熟女少妇av免费看| 老熟女久久久| 999精品在线视频| 久久精品久久精品一区二区三区| 国产成人免费观看mmmm| 91aial.com中文字幕在线观看| 国产熟女欧美一区二区| 日日摸夜夜添夜夜添av毛片| 大又大粗又爽又黄少妇毛片口| av黄色大香蕉| 在线亚洲精品国产二区图片欧美 | 免费大片黄手机在线观看| 久久久欧美国产精品| 免费大片黄手机在线观看| 永久免费av网站大全| 国产熟女欧美一区二区| 日韩熟女老妇一区二区性免费视频| 人妻一区二区av| 精品一品国产午夜福利视频| av黄色大香蕉| 久久精品国产鲁丝片午夜精品| 久久精品国产鲁丝片午夜精品| 免费黄网站久久成人精品| 亚洲图色成人| 久久久久久久亚洲中文字幕| 精品一区二区三卡| 人妻 亚洲 视频| 久久久久久久亚洲中文字幕| 嘟嘟电影网在线观看| 2021少妇久久久久久久久久久| 毛片一级片免费看久久久久| 黄色毛片三级朝国网站| 能在线免费看毛片的网站| 亚洲av不卡在线观看| 日韩强制内射视频| 久久午夜综合久久蜜桃| av天堂久久9| 亚洲精品久久午夜乱码| 18禁在线播放成人免费| 亚洲色图综合在线观看| 午夜福利视频精品| 精品国产国语对白av| videossex国产| 国产精品久久久久久精品电影小说| 777米奇影视久久| 五月伊人婷婷丁香| 熟女电影av网| 日韩视频在线欧美| 日韩一区二区三区影片| 夜夜骑夜夜射夜夜干| 色哟哟·www| 在线观看免费日韩欧美大片 | av一本久久久久| 精品久久国产蜜桃| av一本久久久久| 最近中文字幕高清免费大全6| 国产精品久久久久久av不卡| 性色avwww在线观看| 国产精品嫩草影院av在线观看| 欧美日韩av久久| 国产精品一二三区在线看| 国产爽快片一区二区三区| 日韩一区二区三区影片| 在线观看www视频免费| 免费观看的影片在线观看| 国产成人av激情在线播放 | 亚州av有码| 又大又黄又爽视频免费| 免费观看a级毛片全部| 欧美精品一区二区免费开放| 国产女主播在线喷水免费视频网站| 天堂中文最新版在线下载| 国产成人精品久久久久久| 大话2 男鬼变身卡| 在线看a的网站| 插逼视频在线观看| 美女脱内裤让男人舔精品视频| 久久人妻熟女aⅴ| 成人亚洲精品一区在线观看| 五月开心婷婷网| 欧美激情极品国产一区二区三区 | 国产免费视频播放在线视频| 国产日韩一区二区三区精品不卡 | 久久综合国产亚洲精品| 美女福利国产在线| 亚洲欧美中文字幕日韩二区| 日韩欧美精品免费久久| 国产高清有码在线观看视频| 精品少妇久久久久久888优播| 18禁在线无遮挡免费观看视频| 亚洲精品国产av成人精品| 秋霞在线观看毛片| 一边摸一边做爽爽视频免费| 韩国av在线不卡| 国产欧美日韩一区二区三区在线 | 在线看a的网站| 久久精品国产自在天天线| av免费在线看不卡| 一级毛片电影观看| 午夜免费鲁丝| 夜夜骑夜夜射夜夜干| 麻豆乱淫一区二区| 三上悠亚av全集在线观看| 丁香六月天网| 亚洲成人一二三区av| 亚洲欧美日韩另类电影网站| 日韩 亚洲 欧美在线| 精品亚洲成a人片在线观看| 精品视频人人做人人爽| 午夜精品国产一区二区电影| 亚洲图色成人| 一级黄片播放器| 免费看av在线观看网站| 一本一本综合久久| 丝袜美足系列| 性色avwww在线观看| 亚洲国产精品成人久久小说| 最近的中文字幕免费完整| 欧美 日韩 精品 国产| 在线观看一区二区三区激情| 国产精品国产av在线观看| 成人免费观看视频高清| 永久免费av网站大全| 国产 一区精品| av免费在线看不卡| 免费av中文字幕在线| 国产日韩一区二区三区精品不卡 | 日韩一本色道免费dvd| 精品一品国产午夜福利视频| 九色亚洲精品在线播放| 在线观看免费高清a一片| 99久国产av精品国产电影| 国产黄色视频一区二区在线观看| 韩国av在线不卡| 欧美变态另类bdsm刘玥| 亚洲成人av在线免费| 你懂的网址亚洲精品在线观看| 精品人妻偷拍中文字幕| 中文欧美无线码| 婷婷色麻豆天堂久久| 欧美人与性动交α欧美精品济南到 | 国产精品99久久久久久久久| av有码第一页| 肉色欧美久久久久久久蜜桃| 美女大奶头黄色视频| 男人操女人黄网站| 日韩三级伦理在线观看| 啦啦啦视频在线资源免费观看| 久久久久网色| 亚洲av在线观看美女高潮| 热99久久久久精品小说推荐| 国产有黄有色有爽视频| av福利片在线| 亚洲国产精品一区三区| 日本爱情动作片www.在线观看| av在线播放精品| 美女国产高潮福利片在线看| 18禁动态无遮挡网站| 久久这里有精品视频免费| 久久精品人人爽人人爽视色| av.在线天堂| 日本与韩国留学比较| 七月丁香在线播放| 亚洲情色 制服丝袜| 婷婷色综合www| 香蕉精品网在线| 日韩视频在线欧美| 黄色配什么色好看| 伊人亚洲综合成人网| 日韩精品有码人妻一区| 欧美少妇被猛烈插入视频| 日韩亚洲欧美综合| 亚洲av免费高清在线观看| 午夜激情av网站| 亚洲国产色片| 久久精品夜色国产| 人妻少妇偷人精品九色| a级毛片黄视频| 亚洲av在线观看美女高潮| 久久久久久人妻| 欧美人与善性xxx| 高清视频免费观看一区二区| 大话2 男鬼变身卡| 大香蕉久久成人网| 女人精品久久久久毛片| 亚洲精品久久午夜乱码| 在线免费观看不下载黄p国产| 国产亚洲av片在线观看秒播厂| 免费大片黄手机在线观看| 日产精品乱码卡一卡2卡三| 国产男人的电影天堂91| 国产黄片视频在线免费观看| 亚洲人与动物交配视频| 在线观看一区二区三区激情| 国产白丝娇喘喷水9色精品| av在线观看视频网站免费| 多毛熟女@视频| 制服人妻中文乱码| 少妇猛男粗大的猛烈进出视频| 在线免费观看不下载黄p国产| 久久精品久久久久久久性| av在线老鸭窝| kizo精华| 三级国产精品欧美在线观看| 五月伊人婷婷丁香| 国产成人精品婷婷| 国产免费视频播放在线视频| 99热全是精品| 亚洲国产毛片av蜜桃av| 大香蕉久久网| 日产精品乱码卡一卡2卡三| 一级毛片 在线播放| 亚洲欧美日韩卡通动漫| 久久久久久久大尺度免费视频| 国产日韩欧美在线精品| 99re6热这里在线精品视频| 22中文网久久字幕| 精品久久久久久久久亚洲| √禁漫天堂资源中文www| 欧美日韩国产mv在线观看视频| 久热久热在线精品观看| 在线精品无人区一区二区三| 女人久久www免费人成看片| 国产成人精品在线电影| 五月开心婷婷网| 国产高清不卡午夜福利| .国产精品久久| 久久久久久久久久久免费av| 肉色欧美久久久久久久蜜桃| 三级国产精品欧美在线观看| 大陆偷拍与自拍| videosex国产| 欧美老熟妇乱子伦牲交| 亚洲欧美成人精品一区二区| 久久av网站| 久久女婷五月综合色啪小说| 欧美xxⅹ黑人| 国产一区二区在线观看日韩| 国产亚洲欧美精品永久| 99九九在线精品视频| 精品视频人人做人人爽| 亚洲美女黄色视频免费看| 一区二区三区免费毛片| 男女无遮挡免费网站观看| 亚洲美女黄色视频免费看| 春色校园在线视频观看| 国产毛片在线视频| 91久久精品国产一区二区三区| 夜夜爽夜夜爽视频| 七月丁香在线播放| 欧美日韩亚洲高清精品| 久久99一区二区三区| 欧美少妇被猛烈插入视频| 国产欧美日韩一区二区三区在线 | 一级毛片电影观看| 香蕉精品网在线| 免费观看在线日韩| 人人澡人人妻人| 精品亚洲成a人片在线观看| 国产精品 国内视频| 亚洲精品av麻豆狂野| 久久精品熟女亚洲av麻豆精品| 久久精品国产自在天天线| 亚洲成人一二三区av| 只有这里有精品99| 免费看av在线观看网站| 制服诱惑二区| 麻豆精品久久久久久蜜桃| 日韩人妻高清精品专区| 97精品久久久久久久久久精品| 边亲边吃奶的免费视频| 九色亚洲精品在线播放| 婷婷色麻豆天堂久久| 最近的中文字幕免费完整| 母亲3免费完整高清在线观看 | 成人国语在线视频| 亚洲欧美色中文字幕在线| 一级毛片 在线播放| 97在线人人人人妻| 桃花免费在线播放| 国产日韩欧美在线精品| 日韩av免费高清视频| 男人添女人高潮全过程视频| 国产精品欧美亚洲77777| 精品一品国产午夜福利视频| 国产男人的电影天堂91| 九九爱精品视频在线观看| 超色免费av| 免费日韩欧美在线观看| 精品少妇黑人巨大在线播放| 黑人高潮一二区| 亚洲精品456在线播放app| 成人毛片a级毛片在线播放| 一区二区三区免费毛片| 国产精品99久久久久久久久| 久久久久久久久大av| a级毛色黄片| 在线天堂最新版资源| 人人澡人人妻人| 国产色婷婷99| 26uuu在线亚洲综合色| 国产黄色视频一区二区在线观看| 亚洲精品亚洲一区二区| 丰满饥渴人妻一区二区三| 国产精品国产三级专区第一集| 少妇熟女欧美另类| 制服诱惑二区| 亚洲国产成人一精品久久久| 亚洲国产av影院在线观看| 18禁观看日本| 九九久久精品国产亚洲av麻豆| 99久久人妻综合| 99热网站在线观看| 精品久久久噜噜| 777米奇影视久久| 日韩电影二区| 亚洲欧洲精品一区二区精品久久久 | 日韩,欧美,国产一区二区三区| 视频中文字幕在线观看| 伦理电影大哥的女人| 美女福利国产在线| 夜夜看夜夜爽夜夜摸| 在现免费观看毛片| 日韩av在线免费看完整版不卡| 熟女人妻精品中文字幕| av在线观看视频网站免费| 啦啦啦啦在线视频资源| 美女主播在线视频| 亚洲国产精品国产精品| 日韩大片免费观看网站| 在线观看免费视频网站a站| av一本久久久久| 尾随美女入室| 91精品国产国语对白视频| 免费高清在线观看日韩| 老司机亚洲免费影院| 国产成人精品无人区| 国产 精品1| 亚洲av中文av极速乱| 99热这里只有是精品在线观看| 午夜av观看不卡| 中文字幕免费在线视频6| 亚洲怡红院男人天堂| 午夜福利视频在线观看免费| 亚洲欧洲日产国产| 人人妻人人爽人人添夜夜欢视频| 国产免费视频播放在线视频| 国产精品久久久久久久电影| 欧美日韩av久久| 亚洲精品久久久久久婷婷小说| 国产69精品久久久久777片| 哪个播放器可以免费观看大片| 成年av动漫网址| 日日撸夜夜添| 日韩一区二区视频免费看| 国产视频首页在线观看| 亚洲国产日韩一区二区| 人妻人人澡人人爽人人| 中文天堂在线官网| 欧美3d第一页| 菩萨蛮人人尽说江南好唐韦庄| 亚洲欧美色中文字幕在线| 成人影院久久| 精品少妇黑人巨大在线播放| 精品国产国语对白av| 男女免费视频国产| 国产日韩欧美视频二区| 久久精品久久久久久噜噜老黄| 亚洲少妇的诱惑av| 国产成人精品婷婷| 亚洲国产精品999| 亚洲欧美日韩卡通动漫| 中文字幕人妻丝袜制服| 日韩 亚洲 欧美在线| 91aial.com中文字幕在线观看| 亚洲国产成人一精品久久久| 各种免费的搞黄视频| 久久这里有精品视频免费| 久久综合国产亚洲精品| 国产男女超爽视频在线观看| 亚洲国产精品成人久久小说| 国产一区二区三区av在线| 99热这里只有精品一区| 在线观看三级黄色| 中国美白少妇内射xxxbb| 免费观看性生交大片5| 久久国产精品大桥未久av| 中文乱码字字幕精品一区二区三区| 汤姆久久久久久久影院中文字幕| 99热这里只有精品一区| 午夜激情久久久久久久| 亚洲,欧美,日韩| 亚洲婷婷狠狠爱综合网| 男的添女的下面高潮视频| 欧美3d第一页| 女人精品久久久久毛片| 热99久久久久精品小说推荐| 国产免费视频播放在线视频| 亚洲av日韩在线播放| 麻豆精品久久久久久蜜桃| 日韩精品有码人妻一区| 青春草视频在线免费观看| 国产日韩一区二区三区精品不卡 | 成人影院久久| 国产免费现黄频在线看| a级毛片黄视频| 亚洲美女黄色视频免费看| 亚洲精品中文字幕在线视频| 国产成人91sexporn| 在线观看免费视频网站a站| 午夜免费鲁丝| av.在线天堂| 一级黄片播放器| 日韩精品有码人妻一区| 日韩精品免费视频一区二区三区 | 久久久久久久精品精品| 全区人妻精品视频| 人妻系列 视频| 草草在线视频免费看| 国产日韩欧美在线精品| 国产成人av激情在线播放 | 人人妻人人添人人爽欧美一区卜| 国产成人精品福利久久| 久久久久网色| 五月伊人婷婷丁香| 天天影视国产精品| 人人妻人人添人人爽欧美一区卜| 日本免费在线观看一区| 日本黄色片子视频| 我的女老师完整版在线观看| 国产不卡av网站在线观看| 国产精品久久久久久av不卡| 母亲3免费完整高清在线观看 | 午夜免费鲁丝| 久久人妻熟女aⅴ| 日本黄色片子视频| 成人国语在线视频| 久久av网站| 能在线免费看毛片的网站| 日韩av不卡免费在线播放| 中文字幕免费在线视频6| 国产国拍精品亚洲av在线观看| 久久精品国产亚洲av涩爱| 免费观看的影片在线观看| 一区二区三区精品91| 十八禁高潮呻吟视频| 一级毛片电影观看| 亚洲av二区三区四区| 国产成人aa在线观看| 免费播放大片免费观看视频在线观看| 国产淫语在线视频| av国产精品久久久久影院| 在线看a的网站| 91精品国产九色| 亚洲av福利一区| 精品人妻熟女毛片av久久网站| 大香蕉97超碰在线| 大香蕉久久成人网| 又黄又爽又刺激的免费视频.| 9色porny在线观看| 插逼视频在线观看| 欧美丝袜亚洲另类| 五月玫瑰六月丁香| 国产在线视频一区二区| av女优亚洲男人天堂| 精品久久久久久久久亚洲| 国产淫语在线视频| 午夜福利视频精品| 视频在线观看一区二区三区| 亚洲精品一区蜜桃| 国产高清国产精品国产三级| 插逼视频在线观看| 超碰97精品在线观看| 各种免费的搞黄视频| 99re6热这里在线精品视频| 国产精品熟女久久久久浪| 欧美人与性动交α欧美精品济南到 | 一级片'在线观看视频| 国产欧美亚洲国产| 国产精品99久久99久久久不卡 | 日本91视频免费播放| 黄片播放在线免费| 满18在线观看网站| 亚洲人成网站在线播| 黑人猛操日本美女一级片| 亚洲人成77777在线视频| 国产乱来视频区| 18禁裸乳无遮挡动漫免费视频| 亚洲在久久综合| 亚洲国产色片| 国产精品秋霞免费鲁丝片| 日韩在线高清观看一区二区三区| 亚洲欧美清纯卡通| 哪个播放器可以免费观看大片| 免费人成在线观看视频色| 国产成人a∨麻豆精品| 最新的欧美精品一区二区| 丝袜喷水一区| 综合色丁香网| 成人毛片a级毛片在线播放| 99久久人妻综合| 黑人欧美特级aaaaaa片| 在线看a的网站| 自线自在国产av| 成人18禁高潮啪啪吃奶动态图 | 卡戴珊不雅视频在线播放| 国产黄色免费在线视频| 亚洲国产毛片av蜜桃av| 亚洲精品av麻豆狂野| 亚洲精品视频女| 老女人水多毛片| 99九九线精品视频在线观看视频| 寂寞人妻少妇视频99o| 中文字幕免费在线视频6| 久久ye,这里只有精品| 亚洲第一区二区三区不卡| 狠狠精品人妻久久久久久综合| 国产亚洲最大av| 午夜福利在线观看免费完整高清在| 国产欧美日韩一区二区三区在线 | 伦理电影大哥的女人|