• <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)行水頭處理
    国产精品香港三级国产av潘金莲| 男女之事视频高清在线观看| 亚洲成人免费av在线播放| 一边摸一边做爽爽视频免费| 欧美乱码精品一区二区三区| 亚洲精品av麻豆狂野| 高清在线国产一区| 美女扒开内裤让男人捅视频| 极品人妻少妇av视频| 久久久久久亚洲精品国产蜜桃av| 精品无人区乱码1区二区| 中文字幕人妻熟女乱码| 国产成人影院久久av| 国产精品 国内视频| 中文字幕人妻丝袜制服| 国产一区二区激情短视频| 久久中文字幕人妻熟女| 国产麻豆69| 精品国产乱子伦一区二区三区| 国产精品一区二区精品视频观看| 久热这里只有精品99| 99国产综合亚洲精品| 久久久国产成人免费| 搡老乐熟女国产| 咕卡用的链子| 色在线成人网| 欧美激情高清一区二区三区| 首页视频小说图片口味搜索| 欧美精品啪啪一区二区三区| 国产成人系列免费观看| 精品免费久久久久久久清纯| 亚洲欧美精品综合久久99| netflix在线观看网站| 老司机靠b影院| 欧美老熟妇乱子伦牲交| 大香蕉久久成人网| xxxhd国产人妻xxx| 日韩免费av在线播放| 在线免费观看的www视频| 电影成人av| 女人被狂操c到高潮| 欧美日韩精品网址| 韩国精品一区二区三区| 淫秽高清视频在线观看| 桃色一区二区三区在线观看| 欧美性长视频在线观看| 久久亚洲精品不卡| 亚洲一区二区三区欧美精品| 欧美大码av| 国产精品99久久99久久久不卡| av在线播放免费不卡| 女性生殖器流出的白浆| x7x7x7水蜜桃| 亚洲专区字幕在线| 国产一区在线观看成人免费| 精品福利永久在线观看| 免费av中文字幕在线| 老熟妇乱子伦视频在线观看| 男女做爰动态图高潮gif福利片 | 在线免费观看的www视频| 午夜亚洲福利在线播放| 亚洲成国产人片在线观看| 麻豆一二三区av精品| 身体一侧抽搐| 黄色视频不卡| 人人妻人人澡人人看| 脱女人内裤的视频| 亚洲少妇的诱惑av| 中文字幕另类日韩欧美亚洲嫩草| 人妻久久中文字幕网| 一本大道久久a久久精品| 国产亚洲精品第一综合不卡| 国产精品爽爽va在线观看网站 | 十八禁网站免费在线| 黄片大片在线免费观看| 可以在线观看毛片的网站| 丝袜在线中文字幕| av在线天堂中文字幕 | 午夜激情av网站| avwww免费| 国产精品久久久久久人妻精品电影| 久久久国产一区二区| 老司机亚洲免费影院| 国产在线观看jvid| 十八禁人妻一区二区| 久久人妻av系列| 午夜两性在线视频| 19禁男女啪啪无遮挡网站| 国产精品美女特级片免费视频播放器 | 在线观看午夜福利视频| 色尼玛亚洲综合影院| 久久久久国产精品人妻aⅴ院| 久久国产精品男人的天堂亚洲| 最近最新中文字幕大全免费视频| 最近最新中文字幕大全电影3 | 免费av中文字幕在线| av免费在线观看网站| 80岁老熟妇乱子伦牲交| 精品熟女少妇八av免费久了| 男男h啪啪无遮挡| 99国产极品粉嫩在线观看| 亚洲性夜色夜夜综合| 国产亚洲欧美98| 不卡av一区二区三区| 国产一区在线观看成人免费| 91国产中文字幕| 91老司机精品| 伊人久久大香线蕉亚洲五| 久久久久久大精品| 国产精品久久电影中文字幕| 麻豆成人av在线观看| 日韩欧美三级三区| 99热国产这里只有精品6| 美女福利国产在线| 国产xxxxx性猛交| 久久亚洲真实| 亚洲精品久久午夜乱码| 村上凉子中文字幕在线| 日日干狠狠操夜夜爽| 免费久久久久久久精品成人欧美视频| 亚洲成人免费电影在线观看| 国产成人精品在线电影| 成人18禁在线播放| 成人av一区二区三区在线看| 国产伦人伦偷精品视频| 国产亚洲精品久久久久5区| 性欧美人与动物交配| 岛国在线观看网站| 国产精品一区二区三区四区久久 | 五月开心婷婷网| 久久午夜亚洲精品久久| 国产成+人综合+亚洲专区| 免费少妇av软件| a级毛片在线看网站| 91大片在线观看| 麻豆久久精品国产亚洲av | 欧美日韩瑟瑟在线播放| 久久狼人影院| 午夜福利在线观看吧| 精品第一国产精品| 婷婷精品国产亚洲av在线| 亚洲七黄色美女视频| 多毛熟女@视频| 亚洲人成网站在线播放欧美日韩| 男女床上黄色一级片免费看| 中出人妻视频一区二区| 身体一侧抽搐| 嫩草影院精品99| 欧美成人免费av一区二区三区| 国产精品98久久久久久宅男小说| www.999成人在线观看| 国产欧美日韩综合在线一区二区| 欧美日韩瑟瑟在线播放| 欧美黑人欧美精品刺激| aaaaa片日本免费| 老司机深夜福利视频在线观看| 午夜福利欧美成人| 欧美日韩福利视频一区二区| 色老头精品视频在线观看| 亚洲精品一二三| 18禁黄网站禁片午夜丰满| 欧美乱色亚洲激情| 一区二区日韩欧美中文字幕| 操出白浆在线播放| 欧美日韩精品网址| 一个人免费在线观看的高清视频| 久久伊人香网站| 亚洲一区二区三区色噜噜 | 级片在线观看| 久久欧美精品欧美久久欧美| 亚洲国产欧美一区二区综合| 热re99久久国产66热| 波多野结衣av一区二区av| 国产精品久久久久久人妻精品电影| 最近最新中文字幕大全电影3 | 国产成+人综合+亚洲专区| 日本 av在线| 麻豆国产av国片精品| 91在线观看av| 丰满的人妻完整版| 久99久视频精品免费| 亚洲专区字幕在线| 国产精品二区激情视频| 麻豆国产av国片精品| 国产成人av激情在线播放| 男女下面插进去视频免费观看| 香蕉丝袜av| 日本a在线网址| 高清毛片免费观看视频网站 | √禁漫天堂资源中文www| 日韩欧美免费精品| 日日干狠狠操夜夜爽| 日本a在线网址| 后天国语完整版免费观看| 久久久久久久精品吃奶| 久久欧美精品欧美久久欧美| 99re在线观看精品视频| 黑丝袜美女国产一区| 99香蕉大伊视频| 男女高潮啪啪啪动态图| 久久精品亚洲熟妇少妇任你| 神马国产精品三级电影在线观看 | 亚洲av电影在线进入| 亚洲色图 男人天堂 中文字幕| 国产高清国产精品国产三级| 久久精品亚洲精品国产色婷小说| 女人高潮潮喷娇喘18禁视频| 久久久久久久午夜电影 | 国产高清视频在线播放一区| 中文字幕人妻丝袜制服| 国产激情欧美一区二区| 另类亚洲欧美激情| 精品久久蜜臀av无| 99国产极品粉嫩在线观看| 色精品久久人妻99蜜桃| 91精品国产国语对白视频| 天天影视国产精品| 老司机午夜十八禁免费视频| 黄网站色视频无遮挡免费观看| 国产精品爽爽va在线观看网站 | 黄片大片在线免费观看| 美女福利国产在线| a级片在线免费高清观看视频| 操美女的视频在线观看| 黄网站色视频无遮挡免费观看| 天天添夜夜摸| 久久人人爽av亚洲精品天堂| 亚洲国产毛片av蜜桃av| 99久久国产精品久久久| 妹子高潮喷水视频| 久久亚洲真实| 好看av亚洲va欧美ⅴa在| 欧美午夜高清在线| av片东京热男人的天堂| 老熟妇仑乱视频hdxx| 亚洲av熟女| 免费在线观看日本一区| 午夜福利欧美成人| 国产单亲对白刺激| 波多野结衣高清无吗| 久久精品成人免费网站| 国产精品乱码一区二三区的特点 | 国产成人av教育| 亚洲精品一卡2卡三卡4卡5卡| 久久久国产精品麻豆| 日韩欧美在线二视频| 99久久人妻综合| 亚洲一区二区三区欧美精品| 国产精品乱码一区二三区的特点 | 真人一进一出gif抽搐免费| 欧美国产精品va在线观看不卡| 天堂√8在线中文| 大码成人一级视频| 波多野结衣高清无吗| 免费看十八禁软件| 欧美不卡视频在线免费观看 | av在线播放免费不卡| 午夜精品久久久久久毛片777| 一边摸一边抽搐一进一出视频| 国产精品九九99| √禁漫天堂资源中文www| 国产男靠女视频免费网站| 男人的好看免费观看在线视频 | 中文欧美无线码| 欧美人与性动交α欧美软件| 露出奶头的视频| 中文字幕最新亚洲高清| 波多野结衣av一区二区av| 免费在线观看黄色视频的| 久久国产精品影院| 精品国产乱码久久久久久男人| 国产在线精品亚洲第一网站| 国产精品国产av在线观看| 日本撒尿小便嘘嘘汇集6| 亚洲欧美激情在线| 免费高清视频大片| 黑人猛操日本美女一级片| 午夜免费成人在线视频| 一进一出抽搐动态| 国产精品99久久99久久久不卡| 叶爱在线成人免费视频播放| 久久热在线av| 国产精品久久久人人做人人爽| 桃色一区二区三区在线观看| 97碰自拍视频| 欧美日韩一级在线毛片| 每晚都被弄得嗷嗷叫到高潮| 丝袜人妻中文字幕| 亚洲av五月六月丁香网| 免费女性裸体啪啪无遮挡网站| 久久 成人 亚洲| av欧美777| 两人在一起打扑克的视频| 免费在线观看影片大全网站| 中国美女看黄片| 免费在线观看视频国产中文字幕亚洲| 亚洲性夜色夜夜综合| 国产激情欧美一区二区| 十八禁网站免费在线| 黑人猛操日本美女一级片| 亚洲第一av免费看| 桃色一区二区三区在线观看| 亚洲av成人不卡在线观看播放网| 欧美日韩福利视频一区二区| 纯流量卡能插随身wifi吗| 这个男人来自地球电影免费观看| 欧美日韩瑟瑟在线播放| 满18在线观看网站| 亚洲自偷自拍图片 自拍| 久久精品人人爽人人爽视色| 丁香六月欧美| av有码第一页| 99精品在免费线老司机午夜| 日本精品一区二区三区蜜桃| 夜夜躁狠狠躁天天躁| 两个人免费观看高清视频| 国产亚洲欧美98| 18禁观看日本| 久久久精品国产亚洲av高清涩受| 成人av一区二区三区在线看| 一区在线观看完整版| 亚洲人成网站在线播放欧美日韩| 中国美女看黄片| 欧美激情久久久久久爽电影 | 女人被狂操c到高潮| 999久久久国产精品视频| 久久性视频一级片| 在线观看www视频免费| 国产精品 国内视频| 亚洲va日本ⅴa欧美va伊人久久| 久久久久久久午夜电影 | 男女之事视频高清在线观看| 一本大道久久a久久精品| 亚洲av日韩精品久久久久久密| 亚洲在线自拍视频| 亚洲熟妇中文字幕五十中出 | 琪琪午夜伦伦电影理论片6080| 欧美久久黑人一区二区| 1024香蕉在线观看| 久久午夜综合久久蜜桃| 欧美黑人欧美精品刺激| 黄色a级毛片大全视频| 亚洲激情在线av| 大码成人一级视频| 午夜福利欧美成人| 欧美另类亚洲清纯唯美| 国产一区在线观看成人免费| 成熟少妇高潮喷水视频| 婷婷六月久久综合丁香| 欧美中文综合在线视频| 男女做爰动态图高潮gif福利片 | 大码成人一级视频| 国产一区在线观看成人免费| 国产亚洲欧美精品永久| 一a级毛片在线观看| 女人被狂操c到高潮| 久久欧美精品欧美久久欧美| 久久久久久亚洲精品国产蜜桃av| 好男人电影高清在线观看| 一级作爱视频免费观看| 18禁观看日本| av片东京热男人的天堂| 精品久久久久久,| 99久久精品国产亚洲精品| 波多野结衣高清无吗| bbb黄色大片| ponron亚洲| 18禁黄网站禁片午夜丰满| 国产高清视频在线播放一区| bbb黄色大片| 97碰自拍视频| 两人在一起打扑克的视频| 人人妻,人人澡人人爽秒播| 亚洲欧美激情综合另类| 亚洲成人免费av在线播放| av片东京热男人的天堂| 国产真人三级小视频在线观看| 中文字幕色久视频| 免费少妇av软件| 熟女少妇亚洲综合色aaa.| 丁香欧美五月| 成人18禁在线播放| 国产精品免费一区二区三区在线| 久久人人爽av亚洲精品天堂| 亚洲熟妇熟女久久| 看黄色毛片网站| 一个人免费在线观看的高清视频| 国产熟女xx| 国内毛片毛片毛片毛片毛片| 亚洲va日本ⅴa欧美va伊人久久| 满18在线观看网站| 777久久人妻少妇嫩草av网站| 夜夜爽天天搞| 亚洲中文字幕日韩| 国产熟女午夜一区二区三区| 成人国产一区最新在线观看| 日本wwww免费看| 大型黄色视频在线免费观看| 国产亚洲欧美98| 久久久国产欧美日韩av| 少妇的丰满在线观看| 日本a在线网址| 欧美日韩瑟瑟在线播放| 久久久久久久久久久久大奶| 国产精品亚洲一级av第二区| 真人一进一出gif抽搐免费| 久久婷婷成人综合色麻豆| 99精品久久久久人妻精品| 午夜精品国产一区二区电影| 日本撒尿小便嘘嘘汇集6| 精品久久久久久,| 欧美一级毛片孕妇| 欧美黑人精品巨大| 9色porny在线观看| 老汉色av国产亚洲站长工具| 美女福利国产在线| 天堂中文最新版在线下载| 欧美日韩瑟瑟在线播放| 1024视频免费在线观看| 999精品在线视频| 亚洲一卡2卡3卡4卡5卡精品中文| 成人国产一区最新在线观看| 亚洲第一欧美日韩一区二区三区| 一区二区日韩欧美中文字幕| 丰满迷人的少妇在线观看| 日本黄色日本黄色录像| 色播在线永久视频| 日本a在线网址| 国产精品1区2区在线观看.| av超薄肉色丝袜交足视频| 少妇 在线观看| 如日韩欧美国产精品一区二区三区| 精品国产亚洲在线| 精品熟女少妇八av免费久了| 日韩成人在线观看一区二区三区| 久久久国产一区二区| 午夜福利一区二区在线看| 热99re8久久精品国产| 97碰自拍视频| 一级片'在线观看视频| 亚洲视频免费观看视频| 怎么达到女性高潮| 极品教师在线免费播放| 欧美乱色亚洲激情| 久久精品成人免费网站| 99久久久亚洲精品蜜臀av| 日本vs欧美在线观看视频| 国产成人系列免费观看| 伊人久久大香线蕉亚洲五| 国产黄色免费在线视频| 悠悠久久av| 人人妻人人澡人人看| 欧美av亚洲av综合av国产av| 欧美成人免费av一区二区三区| 久久精品aⅴ一区二区三区四区| 成人影院久久| 美女高潮到喷水免费观看| 欧美日韩国产mv在线观看视频| 黄频高清免费视频| 中文亚洲av片在线观看爽| 久9热在线精品视频| 波多野结衣高清无吗| 精品午夜福利视频在线观看一区| 久久精品国产亚洲av香蕉五月| 欧洲精品卡2卡3卡4卡5卡区| 国产精品一区二区三区四区久久 | 操美女的视频在线观看| avwww免费| 久久人妻av系列| 亚洲av日韩精品久久久久久密| 国产伦一二天堂av在线观看| 久久久久久久久久久久大奶| 不卡一级毛片| 免费在线观看亚洲国产| 男女下面进入的视频免费午夜 | 一区在线观看完整版| 神马国产精品三级电影在线观看 | 久久精品国产清高在天天线| 精品国内亚洲2022精品成人| 国内毛片毛片毛片毛片毛片| 一级a爱视频在线免费观看| 男人操女人黄网站| 欧美日本亚洲视频在线播放| 丰满人妻熟妇乱又伦精品不卡| 日本精品一区二区三区蜜桃| 欧美日韩国产mv在线观看视频| 亚洲av熟女| 精品一区二区三区视频在线观看免费 | xxx96com| 免费高清在线观看日韩| 欧美日韩精品网址| 操出白浆在线播放| 亚洲第一av免费看| 黄色怎么调成土黄色| 国产精品爽爽va在线观看网站 | 一级片'在线观看视频| 欧美老熟妇乱子伦牲交| 亚洲av片天天在线观看| 人成视频在线观看免费观看| 日本黄色视频三级网站网址| 日本撒尿小便嘘嘘汇集6| 日韩精品免费视频一区二区三区| 人人澡人人妻人| 大陆偷拍与自拍| 久久久久国产精品人妻aⅴ院| 男人的好看免费观看在线视频 | 人妻久久中文字幕网| 日韩欧美三级三区| 精品卡一卡二卡四卡免费| 日韩av在线大香蕉| 老司机靠b影院| 首页视频小说图片口味搜索| 日韩欧美一区二区三区在线观看| 久久久久久久久中文| 69av精品久久久久久| 亚洲 欧美 日韩 在线 免费| 最新在线观看一区二区三区| 久久天堂一区二区三区四区| 国产成人精品久久二区二区免费| 色精品久久人妻99蜜桃| 日本免费一区二区三区高清不卡 | 99国产精品一区二区三区| 十八禁人妻一区二区| 午夜免费成人在线视频| 高清在线国产一区| 欧美日韩中文字幕国产精品一区二区三区 | 免费一级毛片在线播放高清视频 | 美女扒开内裤让男人捅视频| 国产成年人精品一区二区 | 黄片大片在线免费观看| 人人妻,人人澡人人爽秒播| 国产片内射在线| av有码第一页| 亚洲色图av天堂| 操出白浆在线播放| 啦啦啦免费观看视频1| 91老司机精品| 满18在线观看网站| 99久久久亚洲精品蜜臀av| 久久久国产一区二区| 欧美日本中文国产一区发布| 国产免费男女视频| 久9热在线精品视频| 在线观看一区二区三区激情| 我的亚洲天堂| 亚洲av熟女| 12—13女人毛片做爰片一| 国产成人欧美在线观看| 99国产精品免费福利视频| 高清欧美精品videossex| a在线观看视频网站| 午夜影院日韩av| 丰满迷人的少妇在线观看| 欧美日韩乱码在线| 欧美日韩黄片免| 首页视频小说图片口味搜索| 亚洲熟妇中文字幕五十中出 | 欧美精品一区二区免费开放| 亚洲国产精品sss在线观看 | 国产高清激情床上av| 看片在线看免费视频| 1024视频免费在线观看| 久久欧美精品欧美久久欧美| 国产aⅴ精品一区二区三区波| 露出奶头的视频| 亚洲五月色婷婷综合| 国产97色在线日韩免费| 88av欧美| 老司机亚洲免费影院| 免费在线观看黄色视频的| 久久精品aⅴ一区二区三区四区| 精品卡一卡二卡四卡免费| 日韩高清综合在线| 涩涩av久久男人的天堂| 精品久久久久久久久久免费视频 | 亚洲av片天天在线观看| 这个男人来自地球电影免费观看| 久久久国产欧美日韩av| 人人妻人人添人人爽欧美一区卜| 久久久久久久精品吃奶| 亚洲专区中文字幕在线| 法律面前人人平等表现在哪些方面| 一个人免费在线观看的高清视频| 国产精品亚洲一级av第二区| 国产97色在线日韩免费| 久久久水蜜桃国产精品网| 欧美日韩av久久| 日本vs欧美在线观看视频| 在线免费观看的www视频| 国产精品久久久久成人av| 国产欧美日韩一区二区精品| avwww免费| 咕卡用的链子| 香蕉国产在线看| 久久狼人影院| 9191精品国产免费久久| 亚洲人成电影观看| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲精品久久午夜乱码| 啪啪无遮挡十八禁网站| 午夜免费成人在线视频| 91老司机精品| 人成视频在线观看免费观看| av中文乱码字幕在线| 淫妇啪啪啪对白视频| 精品福利永久在线观看| 国产精品影院久久| 精品国产乱子伦一区二区三区| 亚洲一码二码三码区别大吗| 成人国产一区最新在线观看| 欧美久久黑人一区二区| 国产成人精品久久二区二区91| 嫩草影视91久久| 我的亚洲天堂| 母亲3免费完整高清在线观看|