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

    急彎段閘孔出流的彎道水流特性研究

    2016-01-18 08:33:06郭紅民,覃閃,蔡黎明
    關(guān)鍵詞:數(shù)值模擬

    急彎段閘孔出流的彎道水流特性研究

    郭紅民1,2,覃閃1,蔡黎明1,曹光春1,張?zhí)锾?,胡文兵1

    (1. 三峽大學(xué) 水利與環(huán)境學(xué)院,湖北宜昌443002; 2. 三峽地區(qū)地質(zhì)災(zāi)害與生態(tài)環(huán)境湖北省協(xié)同創(chuàng)新中心,湖北宜昌443002)

    摘要:在河彎處修建水工建筑物時(shí),研究彎道水流特性對(duì)水工建筑物布置的影響十分重要。采用RNC κ-ε紊流模型對(duì)急彎段閘孔出流的彎道水流進(jìn)行了三維數(shù)值模擬。通過(guò)模擬計(jì)算,比較分析了急彎段不同閘位、不同流量條件下閘孔過(guò)流時(shí)的流速分布、紊流動(dòng)能、壓強(qiáng)分布和環(huán)流等相關(guān)水流特性與變化規(guī)律。分析研究結(jié)果表明,閘壩位于彎道上游時(shí),彎道對(duì)閘壩過(guò)流的影響最小,但需對(duì)閘壩下流彎道凸岸岸坡附近區(qū)域采取一定的防護(hù)措施;閘壩布置在彎道中游時(shí),過(guò)閘水流受彎道的影響最小,閘室水流的穩(wěn)定性、閘門的工作安全以及閘墩的空蝕破壞都是實(shí)際工程中需要注意的重點(diǎn)。數(shù)值模擬驗(yàn)證與物理模型試驗(yàn)結(jié)果吻合較好,計(jì)算精度較高,能夠滿足相關(guān)計(jì)算研究要求。

    關(guān)鍵詞:急彎段; 閘孔出流; 彎道水流; 水流特性; 數(shù)值模擬

    中圖分類號(hào):TV135

    文獻(xiàn)標(biāo)志碼:A

    文章編號(hào):1009-640X(2015)03-0024-07

    Abstract:Owing to the geologyical, topographical and hydrological factors, it is inevitable to build hydraulic structures in the river bend. So it is very important to carry out studies of influences given by the bendway flow on the layout of the hydraulic structures. The three dimensional numerical simulation of the bendway flow outflowing from the sluice orifices located on the rapid river bend has been made by using the RNC κ-ε turbulent model. With the aid of simulation calculation, a comparative analysis is made of the flow characteristics and variation law related to flow velocity distribution, turbulent kinetic energy and pressure distribution as well as circulation flow passing through the sluice under the conditions of different gate openings and different flow discharges. Analysis results show that as the dam sluice is located on the upper river bend, the influence given by the sluice on the flow passing through the dam sluice is the minimal, but some prevention measures should be taken for protecting the area close to the convex bank slope of the lower outlet bend against erosion; and the influence given by the river bend on the flow passing through the sluice is the maximal. Therefore, the stability of the flow flowing through the sluice chamber, the gate safety operation and the cavitation damage to the piers should be taken as the key points in consideration and to which more attention should be paid in practical engineering. The numerical simulation results, which are in good agreement with the results obtained from the physical model tests, can meet the requirements of the calculation studies.

    DOI:10.16198/j.cnki.1009-640X.2015.03.005

    收稿日期:2014-10-30

    基金項(xiàng)目:國(guó)家自然科學(xué)基金資助項(xiàng)目(51178122)

    作者簡(jiǎn)介:張麗娟(1969—), 女, 遼寧遼陽(yáng)人, 副教授, 主要從事軟土工程性質(zhì)、邊坡工程, 樁基礎(chǔ)等方面的教學(xué)與科研工作。E-mail: zhanglijuan1969@126.com

    由于地質(zhì)、地形和水文等因素的要求,水利樞紐工程有時(shí)不可避免要修建在河道彎道處,彎道水流是渠道和河道中很常見(jiàn)的一種水流運(yùn)動(dòng)現(xiàn)象,并且大多都是紊流流動(dòng)。對(duì)彎道水流特性的研究,張紅武等[1]在1984—1986年通過(guò)對(duì)4個(gè)河彎概化模型進(jìn)行了水流運(yùn)動(dòng)規(guī)律的研究,得到了形式簡(jiǎn)單且適用天然河彎的環(huán)流流速垂線分布計(jì)算式;徐桂英[2-3]通過(guò)對(duì)閘前彎道建立物理模型進(jìn)行了一系列的試驗(yàn)研究,得到了彎道水流流速分布和其他相關(guān)特性的規(guī)律;H. J. Vriend[4]對(duì)U型矩形彎道水流的各項(xiàng)特性作了詳細(xì)的試驗(yàn)研究;劉曾美等[5]導(dǎo)出了彎道緩流的縱向和橫向水面曲線計(jì)算式,提出了縱向水面線的簡(jiǎn)易計(jì)算方法;C. C. Fung 等[6]更是提出一種算法,可以自動(dòng)識(shí)別矩形渠道內(nèi)橫向環(huán)流的存在和位置;趙小娥等[7]通過(guò)試驗(yàn)研究了在U形彎道上建閘后水流運(yùn)動(dòng)特性的變化。近年來(lái),還有一些學(xué)者則通過(guò)數(shù)值模擬計(jì)算來(lái)作為理論和試驗(yàn)研究的補(bǔ)充,譚柱林等[8]有效預(yù)測(cè)了明渠180°彎道水流水面高度的變化影響范圍及水面橫比降的主要范圍并分析了不同水深的平面流場(chǎng)和典型斷面彎道環(huán)流分布;鐘杰[9]研究了彎道上低水頭攔河閘閘前的水流特性;R. D. Moser等[10-11]對(duì)小曲率低雷諾數(shù)彎曲河槽中的水流流動(dòng)進(jìn)行了直接數(shù)值模擬。對(duì)于強(qiáng)曲率的急彎河道,由于水流呈較強(qiáng)的三維特性,彎道環(huán)流發(fā)展充分,其在不同流量下對(duì)處于不同彎道位置的水工建筑物布置的影響,目前這方面研究尚不多見(jiàn)。本文利用Flow-3D軟件,采用RNG 紊流模型,對(duì)急彎段閘孔出流的彎道水流進(jìn)行了三維數(shù)值模擬,比較研究了急彎段不同閘位條件下閘孔過(guò)流時(shí)水流特性和變化規(guī)律,可為類似設(shè)計(jì)提供參考。

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

    1.1控制方程

    (1)

    (2)

    式中:ρ為體積分?jǐn)?shù)平均的密度;μ為分子黏性系數(shù);p為壓強(qiáng);μt為紊流黏性系數(shù),μt=ρCμκ2/ε,其中Cμ為經(jīng)驗(yàn)常數(shù),Cμ=0.084 5; μ+μt為廣義黏性系數(shù);ui,uj分別為流速矢量在xi,xj方向的分量。

    1.2紊流模型

    RNGκ-ε模型是基于重整化群(Renormalization Group)的理論提出來(lái)的,經(jīng)過(guò)了改進(jìn)和實(shí)用化處理[12],形式上類似于標(biāo)準(zhǔn)κ-ε模型,但在計(jì)算功能上強(qiáng)于標(biāo)準(zhǔn)κ-ε模型, RNGκ-ε模型考慮了旋轉(zhuǎn)效應(yīng),因此也提高了強(qiáng)旋轉(zhuǎn)流動(dòng)計(jì)算精度,模型中包含了計(jì)算湍流Prandtl數(shù)的解析公式,并且對(duì)近壁區(qū)進(jìn)行適當(dāng)處理后可以計(jì)算低雷諾數(shù)效應(yīng)。

    (3)

    (4)

    1.3數(shù)值方法和自由表面的確定

    采用有限差分格式對(duì)上述控制方程組進(jìn)行差分得到線性方程組,再采用廣義極小殘差算法(GMRES)對(duì)得到的線性方程組進(jìn)行迭代求解。目前,大部分處理水氣兩相流的自由表面都是采用流體體積法(Volume of Fluid VOF)[13],該方法往往會(huì)增加計(jì)算時(shí)間或使最終計(jì)算結(jié)果有較大偏差。本文對(duì)自由表面的描述采用TruVOF方法,該方法不考慮只含氣體的控制單元,只計(jì)算含有液體單元,因此很大程度上減少了模型收斂所需的時(shí)間,對(duì)自由液面描述也更為準(zhǔn)確。TruVOF法對(duì)流體界面的三維重構(gòu)控制方程[14]為:

    (5)

    2幾何模型和計(jì)算驗(yàn)證

    2.1幾何模型及網(wǎng)格劃分

    圖1 彎道幾何模型(單位: cm) Fig.1 A geometrical model for bend (unit: cm)

    彎道幾何模型構(gòu)造分為水槽和閘壩兩部分。彎道水槽底寬60 cm,邊坡系數(shù)m為0.5,水槽底面坡度為0.005,水槽彎道轉(zhuǎn)彎角度為90°,中心轉(zhuǎn)彎半徑為93 cm。閘壩軸線位置共布置3處,閘位1在彎道起點(diǎn),閘位2在彎道轉(zhuǎn)彎45°處,閘位3在彎道終點(diǎn),閘壩皆為5孔平底閘孔,閘底板與水槽底等高,閘壩邊墩上下游均以直立面與水槽邊坡相接。驗(yàn)證流量1.1×107cm3/s,由于流量較小且驗(yàn)證試驗(yàn)閘門并未全開(kāi),驗(yàn)證模型的彎道水流形態(tài)為緩流,各驗(yàn)證斷面留3個(gè)測(cè)壓孔,用測(cè)壓管測(cè)每個(gè)斷面上凸、凹岸邊和中心線水位,矩形堰測(cè)流量,光電流速儀測(cè)垂線縱向流速,采用五點(diǎn)法計(jì)算平均流速,測(cè)量誤差0.01 cm。計(jì)算流量1為3.3×107cm3/s,計(jì)算流量2為5.5×107cm3/s,計(jì)算流量3為7.3×107cm3/s,計(jì)算流量下閘孔內(nèi)的水流形態(tài)均為急流,相對(duì)緩流而言,閘孔內(nèi)的水流更加不平衡,水流對(duì)閘壩的沖擊和空蝕都會(huì)產(chǎn)生不可忽視的危害,且對(duì)下游水流流態(tài)也會(huì)有較大影響。模型布置如圖1。

    本文采用FAVOR法(Fractional Area/V ̄o ̄l ̄u ̄m ̄e ̄ Obstacle Representation)定義幾何體和劃分網(wǎng)格,采用多區(qū)塊結(jié)構(gòu)性網(wǎng)格(Multi-Block Grids)進(jìn)行處理,通過(guò)連接式(Linked block)或者巢狀式(Nested block)區(qū)塊建立網(wǎng)格,有效解決傳統(tǒng)單區(qū)塊結(jié)構(gòu)性網(wǎng)格為求模型細(xì)化而大量增加網(wǎng)格數(shù)量的問(wèn)題。本文網(wǎng)格的劃分采用兩區(qū)塊結(jié)構(gòu)化網(wǎng)格,分別包含上游河道及閘壩和下游河道,通過(guò)連接式區(qū)塊建立網(wǎng)格,網(wǎng)格為均一網(wǎng)格,網(wǎng)格尺寸1 cm。

    上游進(jìn)水口設(shè)置為流量邊界條件,下游出口設(shè)定為自由出流,河道底部及左右側(cè)壁均選擇固壁邊界;渠道頂部空氣入口設(shè)定為對(duì)稱邊界,為零流動(dòng)區(qū)域,且在邊界上無(wú)需計(jì)算剪應(yīng)力,在該邊界上不會(huì)產(chǎn)生對(duì)流或者擴(kuò)散通量。

    2.2模型驗(yàn)證

    驗(yàn)證計(jì)算選擇前述驗(yàn)證流量和多個(gè)驗(yàn)證閘位進(jìn)行,結(jié)果代表以如圖1所示的驗(yàn)證閘位為例,閘孔開(kāi)度e分別為1.5, 2.0, 2.5, 3.0, 3.5和4.0 cm。

    2.2.1流速和主流線驗(yàn)證驗(yàn)證閘位、驗(yàn)證流量下流速和主流線驗(yàn)證計(jì)算結(jié)果見(jiàn)表1,可見(jiàn),計(jì)算結(jié)果與模型試驗(yàn)基本吻合,誤差均小于5%。

    表1 斷面流速和主流線偏離系數(shù)

    2.2.2橫比降和縱比降驗(yàn)證彎頂處的橫比降與閘門相對(duì)開(kāi)度之間的關(guān)系驗(yàn)證見(jiàn)圖2。由圖2可見(jiàn),數(shù)值計(jì)算值和理論計(jì)算值符合文獻(xiàn)[2]相關(guān)線性要求。閘孔開(kāi)度e為3.5 cm時(shí)進(jìn)行縱比降的驗(yàn)證(圖3),可見(jiàn)數(shù)值計(jì)算的結(jié)果與模型試驗(yàn)實(shí)測(cè)值基本一致。驗(yàn)證計(jì)算表明,所建立的三維數(shù)值計(jì)算模型能夠滿足相關(guān)計(jì)算研究要求,計(jì)算精度較高。

    圖2 最大橫降比和閘門開(kāi)度與水深比關(guān)系 Fig.2 Relationships between the maximum transverse gradient and the ratio of gate opening to depth

    圖3 閘門開(kāi)度3.5 cm的彎道處水深 Fig.3 Depth at bend with gate opening 3.5 cm

    3計(jì)算成果分析

    3.1流速分布

    不同計(jì)算流量下不同閘位的閘孔過(guò)流時(shí)的平面流速矢量分布云圖見(jiàn)圖4。閘位1時(shí),彎道位于閘壩下游,閘孔進(jìn)流及閘室過(guò)流受彎道水流的影響不大,進(jìn)閘水流方向與閘墩基本平行,閘孔上游及閘室內(nèi)、兩側(cè)的流速分布基本均勻?qū)ΨQ,閘室內(nèi)平均流速約130 cm/s;水流流出閘孔后受彎道影響,流速重新分布,凸岸流速變大,凹岸流速變小,其差值隨流量增大而增大,最大流速比可達(dá)2.0~2.6。閘位2時(shí),閘壩處于彎道正中,閘孔進(jìn)出流均受彎道水流影響,各流量條件下,主流在進(jìn)閘前與閘墩呈40°~45°的夾角;過(guò)閘水流在閘室內(nèi)和出閘后的流速分布仍受彎道影響,閘室內(nèi)平均流速約為閘位1時(shí)室內(nèi)平均流速的90%,閘孔下游流速分布仍呈現(xiàn)凸岸大,凹岸小的情況,其流速比約為1.3~1.6。閘位3時(shí),彎道位于閘壩上游,主流在進(jìn)閘時(shí)與閘墩呈20°~25°夾角,較閘位2時(shí)主流與閘墩夾角有較大減??;水流在進(jìn)入各閘孔時(shí)與閘墩夾角的規(guī)律為越靠近凹岸夾角越小,越靠近凸岸夾角越大,變化幅度約在5°~55°,閘室內(nèi)平均流速約為閘位1時(shí)室內(nèi)平均流速的80%;閘壩下游渠段雖已為直段,但流速分布仍受前段彎道的影響,呈現(xiàn)凹岸流速大,凸岸流速小的現(xiàn)象,其流速比值約為1.1~1.2。

    圖4 不同閘位的平面流速分布(單位:cm/s) Fig.4 Distribution of plane velocities (unit:cm/s)

    3.2湍流動(dòng)能分布

    計(jì)算流量3時(shí)不同閘位的閘孔過(guò)流時(shí)的平面湍流動(dòng)能分布見(jiàn)圖5。對(duì)于同一閘位,不同流量時(shí)的湍流動(dòng)能分布規(guī)律基本相同,各部位湍流動(dòng)能隨流量增大而增大。閘位1時(shí),閘室內(nèi)湍流動(dòng)能分布均勻且對(duì)稱,最大湍流動(dòng)能約150~190 cm2/s2,閘壩出流下游的湍流動(dòng)能核心區(qū)主要位于凸岸岸坡附近區(qū)域,該區(qū)域的湍流動(dòng)能約580~690 cm2/s2。閘位2時(shí),湍流動(dòng)能核心區(qū)主要分布在閘室內(nèi)的閘墩邊墩附近和閘孔下游凸岸處,閘室內(nèi)的最大湍流動(dòng)能約為閘位1時(shí)閘室內(nèi)最大湍流動(dòng)能的1.9倍,閘壩下游凸岸區(qū)域的湍流動(dòng)能為370~460 cm2/s2。閘位3時(shí)湍流動(dòng)能核心區(qū)主要分布在閘孔下游,閘室內(nèi)的最大湍流動(dòng)能在閘室偏下游處,約為閘位1時(shí)閘室內(nèi)最大湍流動(dòng)能的1.1倍,而該閘位時(shí)閘孔下游的湍流動(dòng)能明顯小于閘位1和閘位2時(shí),該區(qū)域范圍的湍流動(dòng)能為230~310 cm2/s2,是閘位1時(shí)核心區(qū)湍流動(dòng)能的40%~50%,是閘位2時(shí)核心區(qū)湍流動(dòng)能的60%~70%。

    圖5 不同閘位下平面湍流動(dòng)能分布(單位:cm 2/s 2) Fig.5 Distribution of plane turbulent kinetic energy of flow 3 (unit: cm 2/s 2)

    3.3壓強(qiáng)分布

    圖6 不同閘位下平面壓強(qiáng)分布(單位:Pa) Fig.6 Distribution of plane pressure of flow 3 (unit: Pa)

    不同計(jì)算流量下不同閘位閘孔過(guò)流時(shí)的平面壓強(qiáng)分布見(jiàn)圖6。由圖6可見(jiàn),同一閘位不同流量時(shí)的壓強(qiáng)分布規(guī)律基本相同,最大正壓均出現(xiàn)在閘墩頭部。閘位1時(shí)閘孔上游和閘室內(nèi)壓強(qiáng)基本均勻?qū)ΨQ,未出現(xiàn)負(fù)壓,最大正壓約為1.94 kPa。閘位2時(shí)壓強(qiáng)分布規(guī)律為閘孔附近凹岸壓強(qiáng)總體上大于凸岸,凹岸平均壓強(qiáng)約為凸岸平均壓強(qiáng)的1.6倍,最大正壓約1.97 kPa,流量2和流量3在凸岸邊墩下游處出現(xiàn)負(fù)壓,其最大負(fù)壓約-0.59 kPa。閘位3時(shí)最大正壓約1.97 kPa,各計(jì)算流量下凹岸邊墩下游處均出現(xiàn)負(fù)壓,最大負(fù)壓約-0.20 kPa。

    3.4橫向環(huán)流

    流量3閘位3前徑向橫斷面的橫向環(huán)流見(jiàn)圖7。閘位1時(shí)由于閘壩位于彎道上游,其流速分布基本不受彎道影響,其閘前未出現(xiàn)橫向環(huán)流,局部橫向流速主要受到邊墩收縮河道影響而指向河道中心線。閘位2由于處于彎道正中,閘壩在過(guò)流時(shí)受彎道影響較大,閘前的橫斷面出現(xiàn)逆時(shí)針(上游向下游視角,下同)橫向環(huán)流,環(huán)流中心位置偏向凹岸。閘位3時(shí)閘壩雖已處于彎道尾端,但流速分布仍受前段彎道影響,閘前橫斷面仍出現(xiàn)逆時(shí)針橫向環(huán)流,環(huán)流中心位置也偏向凹岸。

    圖7 不同閘位橫向環(huán)流變化 Fig.7 Changes of transverse circulation

    4結(jié)語(yǔ)

    本文采用RNGκ-ε紊流模型,對(duì)急彎段閘孔出流的彎道水流進(jìn)行了三維數(shù)值模擬,數(shù)值模擬驗(yàn)證與物理模型試驗(yàn)結(jié)果吻合較好,計(jì)算精度較高,能夠滿足相關(guān)計(jì)算研究要求。通過(guò)模擬計(jì)算,比較分析了不同流量下不同閘位時(shí)閘壩過(guò)流的部分水力學(xué)特征和規(guī)律,結(jié)果表明:

    (1)閘壩位于彎道上游時(shí),進(jìn)閘水流方向與閘墩基本平行,閘室內(nèi)的流速均勻?qū)ΨQ,閘室內(nèi)的湍流動(dòng)能較小,且閘前未產(chǎn)生橫向環(huán)流,閘墩附近沒(méi)有出現(xiàn)負(fù)壓,但是閘壩出流下游彎道凸岸岸坡附近區(qū)域的湍流動(dòng)能較大,約為閘位2的1.5~1.6倍,閘位3的2.2~2.5倍,實(shí)際工程設(shè)計(jì)中需對(duì)該岸坡采取一定的防護(hù)措施。

    (2)閘壩在彎道中間時(shí),過(guò)閘水流受到彎道的影響較大,閘前有明顯橫向環(huán)流,水流在進(jìn)閘時(shí)與閘墩呈較大夾角約為40°~45°,對(duì)閘墩的橫向沖擊較大,同時(shí)該閘位相對(duì)另外兩種閘位,閘室內(nèi)湍流動(dòng)能也最大,約為閘位2的1.9倍,閘位3的1.8倍,不利于閘室水流穩(wěn)定和閘門安全工作;該閘位閘墩下游處還會(huì)出現(xiàn)較大負(fù)壓,可能對(duì)閘墩產(chǎn)生空蝕破壞;閘壩下游凸岸岸坡附近區(qū)域也出現(xiàn)了較大的湍流動(dòng)能,該處岸坡穩(wěn)定也需注意。

    (3)閘壩在彎道尾端時(shí),閘前也出現(xiàn)了較為明顯的逆時(shí)針橫向環(huán)流,水流在進(jìn)閘時(shí)仍與閘墩呈一定夾角,但較閘位2有所減??;該閘位部分閘墩下游也出現(xiàn)了負(fù)壓,同樣需要防止空蝕破壞。該位置閘壩下游出流相對(duì)平穩(wěn),湍流動(dòng)能值較小,對(duì)岸坡影響較小。

    參考文獻(xiàn):

    [1]張紅武, 呂昕. 彎道水力學(xué)[M]. 北京: 水利電力出版社, 1993. (ZHANG Hong-wu, LV Xi. Hydraulics of curved flow[M]. Beijing: Water Power Press, 1993.(in Chinese))

    [2]徐桂英. 閘前彎道水流特性的試驗(yàn)研究[J]. 武漢水利電力大學(xué)學(xué)報(bào), 1995, 28(3): 338- 341. (XU Gui-ying. Experimental study on the behavior of flow through the curved channel in front of a sluice gate[J]. Journal of Wuhan University of Hydraulic and Electric Engineering, 1995, 28(3): 338- 341. (in Chinese))

    [3]徐桂英. 閘前梯形斷面90°彎道水流垂線縱向流速分布的實(shí)驗(yàn)研究[J]. 中國(guó)農(nóng)村水利水電, 1997(1): 42- 43. (XU Gui-ying. Experimental study on the 90° curved flow of vertical longitudinal velocity at trapezoid section before the gate[J]. China Rural Water and Hydropower, 1997(1): 42- 43. (in Chinese))

    [4]VRIEND H J. Flow measurements in cured rectangular channel: Part 2. Rough bottom[R]. The Netherlands: University of Delft, 1981: 5- 81.

    [5]劉曾美, 吳俊校, 黃國(guó)如. 河渠彎道緩流水面曲線計(jì)算探討[J]. 水利水運(yùn)工程學(xué)報(bào), 2008(2): 54- 59. (LIU Zeng-mei, WU Jun-xiao, HUANG Guo-ru. Computation of water surface curve of suberitical flow in meandering river channel[J]. Hydro-Science and Engineering, 2008(2): 54- 59. (in Chinese))

    [6]FUNG C C, CHUNG K P, CHANDRATILLEKE T. Dentification of secondary flow pattern in a heated curve rectangular channel using image processing technique[C]∥IEEE Instrumentation and Measurement Technology Conference, Institute of Electrical and Electronics Engineers Inc, 2001: 2093- 2098.

    [7]趙小娥, 鐘杰, 劉興年, 等. U形彎道建閘對(duì)水流運(yùn)動(dòng)特性的影響試驗(yàn)[J]. 水利水電科技進(jìn)展, 2010, 30(3): 47- 49, 57. (ZHAO Xiao-e, ZHONG Jie, LIU Xing-nian, et al. Experimental study on influences of sluice construction in a U-shape bend on characteristics of flow motion[J]. Advances in Science and Technology of Water Resources, 2010, 30(3): 47- 49, 57. (in Chinese))

    [8]譚柱林, 彭?xiàng)? 明渠彎道水流三維數(shù)值模擬[J]. 水運(yùn)工程, 2012(3): 46- 49. (TAN Zhu-lin, PENG Yang. Three-dimensional numerical simulation of bend flow in open channel[J]. Port & Waterway Engineering, 2012(3): 46- 49. (in Chinese))

    [9]鐘杰. 彎道低水頭攔河閘閘孔分流比試驗(yàn)研究[J]. 長(zhǎng)江科學(xué)院院報(bào), 2011, 28(3): 28- 32. (ZHONG Jie. Experinental study on flow diversion ratio of low head sluice built in curved channel[J]. Journal of Yangtze River Scientific Research Institute, 2011, 28(3): 28- 32. (in Chinese))

    [10]MOSER R D, MOIN P. Direct numerical simulation of curved turbulent channel flow[R]. CA: NASA Ames Research Center, 1984.

    [11]MOSER R D, MOIN P. The effects of curvature in wall-bounded turbulent flow[J]. Journal of Fluid Mechanics, 1987, 175: 479- 510.

    [12]RIJN L C.Sediment transport, Part I: bed load transport[J]. Journal of Hydraulic Engineering, 1984, 110(10): 1431- 1456.

    [13]HIR C W, NICHOILS B D. Volume of fluid(VOF) for the dynamics of free boundaries[J]. Journal of Computational Physics, 1981, 39(1): 201- 225.

    [14]MATTHEWS B W, FLETCHER C A J, PARTRIDGE A C, et al. Computations of curved free surface water flow on spiral concentrators[J]. Canadian Metallurgical Quarterly, 1999, 125(11): 1126- 1139.

    A study of bendway flow characteristics outflowing from sluice orifice on rapid river bend

    GUO Hong-min1,2, QIN Shan1, CAI Li-ming1, CAO Guang-chun1, ZHANG Tian-tian1, HU Wen-bing1

    (1.CollegeofHydraulicandEnvironmentalEngineering,ChinaThreeGorgesUniversity,Yichang443002,China; 2.CollaborativeInnovationCenterforGeo-HazardsandEco-EnvironmentinThreeGorgesArea,Yichang443002,China)

    Key words: rapid river bend; outflow from sluice orifice; bendway flow; flow characteristics; numerical simulation

    張麗娟, 劉仁釗. 南沙港淤泥固化前后物理力學(xué)性能和微觀結(jié)構(gòu)變化[J]. 水利水運(yùn)工程學(xué)報(bào), 2015(3): 31-36. (ZHANG Li-juan, LIU Ren-zhao. Physico-mechanical properties and changes in microstructure of silt before and after solidification in Nansha port[J]. Hydro-Science and Engineering, 2015(3): 31-36.)

    猜你喜歡
    數(shù)值模擬
    基于AMI的雙色注射成型模擬分析
    錐齒輪精密冷擺輾成形在“材料成型數(shù)值模擬”課程教學(xué)中的應(yīng)用
    基于氣象信息及風(fēng)場(chǎng)信息的風(fēng)機(jī)輪轂處風(fēng)速預(yù)測(cè)
    鉆孔灌注樁樁底沉渣對(duì)樁體承載特性影響的模擬分析
    西南地區(qū)氣象資料測(cè)試、預(yù)處理和加工研究報(bào)告
    科技資訊(2016年18期)2016-11-15 08:01:18
    張家灣煤礦巷道無(wú)支護(hù)條件下位移的數(shù)值模擬
    科技視界(2016年18期)2016-11-03 23:14:27
    張家灣煤礦開(kāi)切眼錨桿支護(hù)參數(shù)確定的數(shù)值模擬
    科技視界(2016年18期)2016-11-03 22:57:21
    跨音速飛行中機(jī)翼水汽凝結(jié)的數(shù)值模擬研究
    科技視界(2016年18期)2016-11-03 20:38:17
    姚橋煤礦采空區(qū)CO2防滅火的數(shù)值模擬分析
    雙螺桿膨脹機(jī)的流場(chǎng)數(shù)值模擬研究
    科技視界(2016年22期)2016-10-18 14:53:19
    av免费在线看不卡| 久久 成人 亚洲| 亚洲国产精品成人久久小说| 乱码一卡2卡4卡精品| 国产成人午夜福利电影在线观看| 国产黄片美女视频| 国产亚洲5aaaaa淫片| 精品久久久久久电影网| 色视频在线一区二区三区| 国语对白做爰xxxⅹ性视频网站| 国产中年淑女户外野战色| 晚上一个人看的免费电影| 免费人成在线观看视频色| 国产成人精品福利久久| 国产精品偷伦视频观看了| 日韩中文字幕视频在线看片| 国模一区二区三区四区视频| 九九爱精品视频在线观看| 国产伦在线观看视频一区| 国产成人精品久久久久久| 国产精品偷伦视频观看了| 在线观看三级黄色| 国产亚洲精品久久久com| 精品99又大又爽又粗少妇毛片| 成人午夜精彩视频在线观看| 纵有疾风起免费观看全集完整版| 一本色道久久久久久精品综合| 少妇猛男粗大的猛烈进出视频| 国产黄色视频一区二区在线观看| 大又大粗又爽又黄少妇毛片口| 精品酒店卫生间| 亚洲一区二区三区欧美精品| 2018国产大陆天天弄谢| 美女福利国产在线| 亚洲图色成人| 久久鲁丝午夜福利片| 亚洲av国产av综合av卡| 日韩强制内射视频| a 毛片基地| videos熟女内射| 亚洲国产日韩一区二区| 国产亚洲5aaaaa淫片| 又粗又硬又长又爽又黄的视频| 美女xxoo啪啪120秒动态图| 国内精品宾馆在线| 精品亚洲成国产av| 波野结衣二区三区在线| 日韩不卡一区二区三区视频在线| 天美传媒精品一区二区| 在线观看www视频免费| 亚洲三级黄色毛片| 中文天堂在线官网| 男女边吃奶边做爰视频| 国产视频内射| 精品少妇黑人巨大在线播放| 欧美日韩在线观看h| 在线观看三级黄色| 成年人午夜在线观看视频| 91精品国产国语对白视频| 狂野欧美激情性bbbbbb| 日韩欧美一区视频在线观看 | 天堂8中文在线网| 插阴视频在线观看视频| 丁香六月天网| 肉色欧美久久久久久久蜜桃| 97超视频在线观看视频| 国产亚洲精品久久久com| 久久精品国产亚洲av天美| 亚洲,欧美,日韩| 成人毛片60女人毛片免费| 最新中文字幕久久久久| 国产成人免费无遮挡视频| 免费在线观看成人毛片| 亚洲欧美成人精品一区二区| 这个男人来自地球电影免费观看 | 黄色怎么调成土黄色| 男人和女人高潮做爰伦理| 又爽又黄a免费视频| 久久久久国产网址| 3wmmmm亚洲av在线观看| 一二三四中文在线观看免费高清| 99久久人妻综合| 国产一区二区三区综合在线观看 | av.在线天堂| 一区二区av电影网| 久久国内精品自在自线图片| 大又大粗又爽又黄少妇毛片口| 最近中文字幕2019免费版| 免费观看a级毛片全部| 偷拍熟女少妇极品色| 一本一本综合久久| 日韩一区二区视频免费看| 你懂的网址亚洲精品在线观看| 成人黄色视频免费在线看| 亚洲伊人久久精品综合| 丰满少妇做爰视频| 国产精品麻豆人妻色哟哟久久| 欧美最新免费一区二区三区| 亚洲精品乱久久久久久| 97精品久久久久久久久久精品| 青春草视频在线免费观看| 亚洲国产精品一区二区三区在线| 国产亚洲最大av| 亚洲国产色片| 80岁老熟妇乱子伦牲交| 国产黄频视频在线观看| 一级毛片久久久久久久久女| 狠狠精品人妻久久久久久综合| 免费看光身美女| 亚洲熟女精品中文字幕| 精品少妇内射三级| 欧美97在线视频| 人妻夜夜爽99麻豆av| 亚洲精品色激情综合| 国产精品嫩草影院av在线观看| tube8黄色片| 国产有黄有色有爽视频| 午夜久久久在线观看| 人妻少妇偷人精品九色| 欧美高清成人免费视频www| 2022亚洲国产成人精品| 免费观看无遮挡的男女| 乱码一卡2卡4卡精品| 国产白丝娇喘喷水9色精品| 一区在线观看完整版| 青春草亚洲视频在线观看| 简卡轻食公司| 国产精品人妻久久久久久| 边亲边吃奶的免费视频| 一级a做视频免费观看| 老司机亚洲免费影院| 国产一区有黄有色的免费视频| av福利片在线观看| 国内揄拍国产精品人妻在线| 一本久久精品| 在线观看一区二区三区激情| 亚洲av国产av综合av卡| 亚洲av电影在线观看一区二区三区| 亚洲欧美一区二区三区国产| 人人妻人人爽人人添夜夜欢视频 | 91在线精品国自产拍蜜月| 在线观看www视频免费| a级片在线免费高清观看视频| 纵有疾风起免费观看全集完整版| 久久久久人妻精品一区果冻| 少妇熟女欧美另类| 国产精品99久久久久久久久| 视频区图区小说| 国产成人a∨麻豆精品| 中文字幕制服av| 亚洲无线观看免费| 婷婷色综合大香蕉| 国产免费视频播放在线视频| 大又大粗又爽又黄少妇毛片口| 91午夜精品亚洲一区二区三区| 亚洲av电影在线观看一区二区三区| 一级av片app| 中文资源天堂在线| 亚洲国产精品专区欧美| 欧美精品亚洲一区二区| 国产欧美亚洲国产| 日韩中文字幕视频在线看片| 免费看av在线观看网站| 亚洲av二区三区四区| 十分钟在线观看高清视频www | 亚洲av二区三区四区| 成人国产av品久久久| 晚上一个人看的免费电影| 国产男女超爽视频在线观看| 97在线视频观看| 青春草国产在线视频| 欧美人与善性xxx| 国产一区二区在线观看av| 亚洲综合精品二区| 国产精品不卡视频一区二区| 精品国产一区二区久久| 欧美xxxx性猛交bbbb| 人人妻人人爽人人添夜夜欢视频 | 99热全是精品| 五月开心婷婷网| 精品一区二区免费观看| 最近2019中文字幕mv第一页| 亚洲欧美日韩卡通动漫| 大又大粗又爽又黄少妇毛片口| 日韩人妻高清精品专区| 激情五月婷婷亚洲| 欧美成人午夜免费资源| 久久久亚洲精品成人影院| 久久综合国产亚洲精品| 老熟女久久久| 久久99蜜桃精品久久| 亚洲婷婷狠狠爱综合网| 午夜视频国产福利| 欧美精品一区二区免费开放| 性色avwww在线观看| 男人添女人高潮全过程视频| 午夜av观看不卡| 一区二区三区乱码不卡18| 午夜免费观看性视频| 国产熟女午夜一区二区三区 | 精品亚洲成a人片在线观看| 中文天堂在线官网| 亚洲色图综合在线观看| 一级毛片aaaaaa免费看小| 一区在线观看完整版| 最近手机中文字幕大全| 亚洲国产精品999| 国产女主播在线喷水免费视频网站| 有码 亚洲区| 丝袜脚勾引网站| 久久久a久久爽久久v久久| 国产精品国产三级专区第一集| 美女福利国产在线| 久久亚洲国产成人精品v| 亚洲怡红院男人天堂| 黄色视频在线播放观看不卡| 成人亚洲欧美一区二区av| 欧美国产精品一级二级三级 | 高清av免费在线| 超碰97精品在线观看| 99精国产麻豆久久婷婷| 久久影院123| 蜜桃久久精品国产亚洲av| 欧美日韩一区二区视频在线观看视频在线| 91久久精品国产一区二区三区| 亚洲人成网站在线播| av网站免费在线观看视频| 一区二区三区四区激情视频| 国产又色又爽无遮挡免| av又黄又爽大尺度在线免费看| 一级av片app| 亚洲精品aⅴ在线观看| 亚洲,一卡二卡三卡| 日韩电影二区| 日本vs欧美在线观看视频 | 草草在线视频免费看| 97精品久久久久久久久久精品| 国产高清不卡午夜福利| 久久 成人 亚洲| 另类亚洲欧美激情| 精品午夜福利在线看| 欧美日韩亚洲高清精品| 91午夜精品亚洲一区二区三区| 久久午夜福利片| a级片在线免费高清观看视频| 十分钟在线观看高清视频www | 婷婷色麻豆天堂久久| 一本久久精品| 波野结衣二区三区在线| 日日爽夜夜爽网站| 高清视频免费观看一区二区| av一本久久久久| 中文字幕人妻熟人妻熟丝袜美| 亚洲精品,欧美精品| 2021少妇久久久久久久久久久| 欧美激情国产日韩精品一区| 我要看日韩黄色一级片| a级一级毛片免费在线观看| 久久久久久久久大av| 99热这里只有是精品在线观看| 日本黄大片高清| 国产免费视频播放在线视频| 国产精品99久久久久久久久| √禁漫天堂资源中文www| 韩国av在线不卡| 国产在线免费精品| 日日啪夜夜撸| 欧美精品亚洲一区二区| 少妇高潮的动态图| 亚洲av电影在线观看一区二区三区| 精品99又大又爽又粗少妇毛片| 黑丝袜美女国产一区| 人人妻人人添人人爽欧美一区卜| 久久精品国产亚洲av天美| 91aial.com中文字幕在线观看| 久久午夜福利片| 久久这里有精品视频免费| 亚洲av.av天堂| 熟女人妻精品中文字幕| 看十八女毛片水多多多| 丰满迷人的少妇在线观看| 国产一级毛片在线| 一级毛片黄色毛片免费观看视频| 亚洲成人手机| av视频免费观看在线观看| 亚洲精品色激情综合| 99热这里只有是精品在线观看| 精品酒店卫生间| 一级毛片久久久久久久久女| 新久久久久国产一级毛片| 天堂中文最新版在线下载| 美女福利国产在线| 在线亚洲精品国产二区图片欧美 | 亚洲av日韩在线播放| 久久久久人妻精品一区果冻| 久久婷婷青草| 亚洲av成人精品一二三区| 国产精品久久久久久精品古装| 最近中文字幕2019免费版| 高清不卡的av网站| 赤兔流量卡办理| 国产亚洲91精品色在线| 狠狠精品人妻久久久久久综合| av播播在线观看一区| 18禁在线无遮挡免费观看视频| 欧美老熟妇乱子伦牲交| 日韩三级伦理在线观看| 精品亚洲成国产av| 免费观看a级毛片全部| a级毛色黄片| 国产精品一区二区在线观看99| 9色porny在线观看| 精品一区二区三卡| 中文精品一卡2卡3卡4更新| 三级国产精品片| 国产成人freesex在线| 成人漫画全彩无遮挡| 精品亚洲成a人片在线观看| 亚洲av成人精品一二三区| 熟女人妻精品中文字幕| 国产成人精品久久久久久| 天美传媒精品一区二区| 亚洲av不卡在线观看| 免费人妻精品一区二区三区视频| 午夜91福利影院| 黄色配什么色好看| 国产成人精品婷婷| 国产高清三级在线| 最新中文字幕久久久久| 亚洲精品成人av观看孕妇| 在线看a的网站| 久久午夜综合久久蜜桃| 亚洲无线观看免费| 亚洲一级一片aⅴ在线观看| 亚洲图色成人| 国产一区二区三区av在线| 国产精品一区www在线观看| 人妻制服诱惑在线中文字幕| 丰满迷人的少妇在线观看| av在线老鸭窝| 日韩三级伦理在线观看| 日本-黄色视频高清免费观看| 一区二区三区免费毛片| 国产毛片在线视频| 一边亲一边摸免费视频| a级一级毛片免费在线观看| 久久精品久久久久久噜噜老黄| 如何舔出高潮| 国内少妇人妻偷人精品xxx网站| 美女xxoo啪啪120秒动态图| 老女人水多毛片| 国产成人91sexporn| 99热6这里只有精品| 国产白丝娇喘喷水9色精品| 国产高清国产精品国产三级| 18禁动态无遮挡网站| 国产熟女午夜一区二区三区 | 美女视频免费永久观看网站| 欧美精品人与动牲交sv欧美| av女优亚洲男人天堂| 看十八女毛片水多多多| 午夜老司机福利剧场| 国产白丝娇喘喷水9色精品| 18禁在线播放成人免费| 亚洲人成网站在线观看播放| 777米奇影视久久| 亚洲av成人精品一区久久| 亚洲熟女精品中文字幕| 99热这里只有精品一区| 国产精品成人在线| 晚上一个人看的免费电影| 深夜a级毛片| 国产精品一区www在线观看| 久久久久久久久久久久大奶| 亚洲综合精品二区| 久久99一区二区三区| 一二三四中文在线观看免费高清| 久久国产精品大桥未久av | 日韩欧美一区视频在线观看 | 国产一区有黄有色的免费视频| 爱豆传媒免费全集在线观看| 亚洲国产成人一精品久久久| 成人国产麻豆网| 色视频www国产| 一区二区三区精品91| 免费观看的影片在线观看| 中文字幕人妻丝袜制服| 国产精品久久久久久av不卡| 在线天堂最新版资源| 高清不卡的av网站| 欧美变态另类bdsm刘玥| 一级片'在线观看视频| 精品一品国产午夜福利视频| av又黄又爽大尺度在线免费看| 成人综合一区亚洲| 看非洲黑人一级黄片| 最近2019中文字幕mv第一页| 欧美精品一区二区免费开放| 欧美 日韩 精品 国产| 天天躁夜夜躁狠狠久久av| av卡一久久| 国产探花极品一区二区| 国产日韩一区二区三区精品不卡 | 久久精品熟女亚洲av麻豆精品| 亚洲国产av新网站| 纯流量卡能插随身wifi吗| 国产欧美日韩综合在线一区二区 | 国产在线免费精品| 久久午夜综合久久蜜桃| 国产在线视频一区二区| a级毛片免费高清观看在线播放| 国产精品人妻久久久久久| 亚洲欧美精品专区久久| 亚洲综合精品二区| av网站免费在线观看视频| 99久久综合免费| 国产毛片在线视频| 一区二区三区四区激情视频| 午夜视频国产福利| 精品久久久久久久久av| 另类亚洲欧美激情| 国产男人的电影天堂91| 中文字幕久久专区| 各种免费的搞黄视频| 美女视频免费永久观看网站| 欧美日韩精品成人综合77777| 在现免费观看毛片| 国产在线免费精品| 国产精品三级大全| 少妇精品久久久久久久| 九色成人免费人妻av| 亚洲av不卡在线观看| 久久99精品国语久久久| 国产欧美日韩一区二区三区在线 | 日韩大片免费观看网站| 久久久久久人妻| 精品一区二区免费观看| 国产91av在线免费观看| 69精品国产乱码久久久| 亚洲美女搞黄在线观看| 91aial.com中文字幕在线观看| 国产欧美日韩一区二区三区在线 | 国产亚洲午夜精品一区二区久久| 人人妻人人看人人澡| av又黄又爽大尺度在线免费看| 三上悠亚av全集在线观看 | 在线观看人妻少妇| 日韩制服骚丝袜av| 中文字幕人妻熟人妻熟丝袜美| 久久精品国产自在天天线| 国产老妇伦熟女老妇高清| 中文天堂在线官网| 欧美日韩av久久| 久久久久国产网址| 亚洲精品国产av成人精品| 欧美高清成人免费视频www| 国产成人免费观看mmmm| 亚洲国产欧美日韩在线播放 | 51国产日韩欧美| 国产av精品麻豆| 老司机亚洲免费影院| 插逼视频在线观看| 少妇 在线观看| 成人毛片a级毛片在线播放| 男人和女人高潮做爰伦理| 最新的欧美精品一区二区| 麻豆成人午夜福利视频| 最近中文字幕2019免费版| 简卡轻食公司| 在线观看av片永久免费下载| 99热全是精品| 22中文网久久字幕| 亚洲va在线va天堂va国产| 在线观看av片永久免费下载| 水蜜桃什么品种好| 91aial.com中文字幕在线观看| 一级黄片播放器| 大话2 男鬼变身卡| 国产无遮挡羞羞视频在线观看| 免费高清在线观看视频在线观看| 久久久久久人妻| 婷婷色综合大香蕉| 中文乱码字字幕精品一区二区三区| 一级,二级,三级黄色视频| 欧美日韩av久久| 欧美日韩视频精品一区| 99热6这里只有精品| 最新中文字幕久久久久| 国产91av在线免费观看| 久久久国产一区二区| 男女无遮挡免费网站观看| 在线观看三级黄色| 午夜福利,免费看| 久久久精品94久久精品| 三级经典国产精品| 在线观看免费高清a一片| 黄色一级大片看看| av不卡在线播放| 国产精品久久久久久久电影| 亚洲精品日本国产第一区| 伊人亚洲综合成人网| 欧美精品亚洲一区二区| 少妇被粗大猛烈的视频| 国产91av在线免费观看| 免费久久久久久久精品成人欧美视频 | 一个人免费看片子| 国产精品免费大片| 少妇被粗大的猛进出69影院 | 久久av网站| av卡一久久| 欧美日本中文国产一区发布| 日韩成人av中文字幕在线观看| av天堂中文字幕网| 亚洲高清免费不卡视频| 久久久亚洲精品成人影院| 国产成人aa在线观看| 精品人妻偷拍中文字幕| 精品99又大又爽又粗少妇毛片| 26uuu在线亚洲综合色| 亚洲欧美中文字幕日韩二区| 国产午夜精品久久久久久一区二区三区| 建设人人有责人人尽责人人享有的| 国产精品秋霞免费鲁丝片| 狠狠精品人妻久久久久久综合| 国内揄拍国产精品人妻在线| 99热这里只有是精品50| 久久久国产一区二区| 久久久久久久久久成人| 又粗又硬又长又爽又黄的视频| 夜夜爽夜夜爽视频| 激情五月婷婷亚洲| 黑丝袜美女国产一区| 嫩草影院新地址| 大话2 男鬼变身卡| 国产淫片久久久久久久久| 久久国产乱子免费精品| 亚洲国产成人一精品久久久| 欧美精品国产亚洲| 男女无遮挡免费网站观看| 高清不卡的av网站| 国产精品一区二区在线观看99| 最近的中文字幕免费完整| 日日啪夜夜撸| a级毛色黄片| 18禁裸乳无遮挡动漫免费视频| 亚洲av成人精品一区久久| 亚洲精品久久午夜乱码| 色94色欧美一区二区| 亚洲欧美日韩另类电影网站| av在线观看视频网站免费| 久久久a久久爽久久v久久| 日韩大片免费观看网站| 在线精品无人区一区二区三| av国产精品久久久久影院| 女人久久www免费人成看片| 最近2019中文字幕mv第一页| 国产精品国产三级专区第一集| 涩涩av久久男人的天堂| 国产亚洲欧美精品永久| 少妇高潮的动态图| 成人无遮挡网站| 国产精品国产三级专区第一集| 久久精品久久久久久久性| 国产黄频视频在线观看| 黄色欧美视频在线观看| 丰满乱子伦码专区| 黑丝袜美女国产一区| 老司机影院成人| 五月伊人婷婷丁香| 涩涩av久久男人的天堂| 日韩成人av中文字幕在线观看| www.av在线官网国产| 欧美日韩av久久| 2022亚洲国产成人精品| 国产精品麻豆人妻色哟哟久久| 一边亲一边摸免费视频| 亚洲丝袜综合中文字幕| av国产久精品久网站免费入址| 免费观看av网站的网址| 女的被弄到高潮叫床怎么办| 一区二区av电影网| 中文资源天堂在线| 黄色毛片三级朝国网站 | 六月丁香七月| 在线观看免费日韩欧美大片 | 老女人水多毛片| 777米奇影视久久| 久久久久久久久久久免费av| 欧美 亚洲 国产 日韩一| 精品久久久久久久久亚洲| 精品亚洲乱码少妇综合久久| 国产高清国产精品国产三级| 人妻制服诱惑在线中文字幕| 51国产日韩欧美| 日产精品乱码卡一卡2卡三| 如何舔出高潮| av一本久久久久| 色5月婷婷丁香| 亚洲va在线va天堂va国产| 纵有疾风起免费观看全集完整版| 乱人伦中国视频| 国产成人精品婷婷| av福利片在线观看| 亚洲国产色片| 日日撸夜夜添| 午夜免费观看性视频| 午夜福利影视在线免费观看| a级毛片在线看网站| 久久久a久久爽久久v久久| 国产亚洲最大av| 日本黄大片高清| 99热6这里只有精品| 国产精品国产三级国产av玫瑰| 精品视频人人做人人爽| 午夜老司机福利剧场|