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

    基于CFD 的泵站進(jìn)水流道優(yōu)化

    2024-01-01 00:00:00田堯天江全星彭昆萬川羅莎陳慧符向前
    農(nóng)業(yè)工程 2024年6期
    關(guān)鍵詞:計(jì)算流體力學(xué)數(shù)值模擬泵站

    摘 要: 受限于用地規(guī)模、地形等因素,我國排水泵站流道的設(shè)計(jì)趨于復(fù)雜,普遍具有流道深、狹長的特點(diǎn)。這導(dǎo)致泵站進(jìn)水流道存在回流、流速分布不均勻等不良流態(tài),影響泵站安全穩(wěn)定運(yùn)行。以某泵站為實(shí)例,采用計(jì)算流體力學(xué)(computational fluid dynamics,CFD)技術(shù),通過改變肘形流道的彎曲段長度、底邊線傾角和喉管高度,擬定優(yōu)化設(shè)計(jì)方案,采用基于CFD 的數(shù)值模擬方法分析流道內(nèi)部流場分布,并根據(jù)流道出口的流速分布均勻度、速度加權(quán)平均角度和流道水力損失等水力指標(biāo)比選,優(yōu)選出最優(yōu)方案。該研究為進(jìn)水流道設(shè)計(jì)提供了更精細(xì)的優(yōu)化手段,彌補(bǔ)了傳統(tǒng)經(jīng)驗(yàn)設(shè)計(jì)方法只能定性分析的不足,得到的最優(yōu)設(shè)計(jì)方案可供類似工程參考,保障泵站機(jī)組的穩(wěn)定運(yùn)行。

    關(guān)鍵詞:泵站;流道設(shè)計(jì);流態(tài)優(yōu)化;數(shù)值模擬;計(jì)算流體力學(xué)

    中圖分類號:S277 文獻(xiàn)標(biāo)識(shí)碼:A 文章編號:2095-1795(2024)06-0092-08

    DOI:10.19998/j.cnki.2095-1795.2024.06.016

    0 引言

    泵站進(jìn)水流道由進(jìn)水池、前池和吸水室組成,前池是連接進(jìn)水池和吸水室的泵站建筑物,前池設(shè)計(jì)的合理與否關(guān)系到水流平穩(wěn)、均勻地進(jìn)入吸水室,否則影響泵站機(jī)組的穩(wěn)定與安全[1-3]。泵站進(jìn)水流道不合理會(huì)導(dǎo)致回流、橫向流動(dòng)等不良流態(tài),引發(fā)氣蝕,影響水泵運(yùn)行安全和工作效率[4]。我國泵站的設(shè)計(jì)與建造以《泵站設(shè)計(jì)規(guī)范》為依據(jù)[5]。在設(shè)計(jì)前池進(jìn)水流道時(shí),通常要考慮將流道渠線設(shè)計(jì)得盡可能順直。目前泵站的布局受到用地指標(biāo)、地形等因素的限制,進(jìn)水流道的設(shè)計(jì)偏狹小,多彎道,水流在進(jìn)水流道中不能充分?jǐn)U散,容易產(chǎn)生渦流、水面超高、流速分布不均等不良流態(tài)[6]。因此,在泵站設(shè)計(jì)的過程中分析其進(jìn)水流道流態(tài)和水力特性,并提出合理的整流方案對其進(jìn)行優(yōu)化是必要的。對于改善泵站前池進(jìn)水流道的水流條件,提高泵站機(jī)組運(yùn)行效率、延長水泵使用壽命,以及節(jié)約工程造價(jià)具有重大意義[7-12]。

    隨著計(jì)算流體力學(xué)(computational fluid dynamics,CFD)技術(shù)的發(fā)展,越來越多的研究人員將 CFD 與實(shí)際工程設(shè)計(jì)相聯(lián)系,用計(jì)算機(jī)模擬技術(shù)對泵站進(jìn)水流道內(nèi)流態(tài)進(jìn)行研究。資丹等[13] 對廣東省永湖泵站進(jìn)行了數(shù)值計(jì)算和現(xiàn)場測試,研究了組合式導(dǎo)流墩在改善大型泵站前池、進(jìn)水池流態(tài)方面的效果。張睿等[14] 基于計(jì)算流體動(dòng)力學(xué),研究了斜向管涵進(jìn)流對泵站水力流態(tài)特性的影響,并發(fā)現(xiàn)采用分流墩、組合梁及相背布置短導(dǎo)流墩的組合式整流措施可以顯著改善不良流態(tài)。李志祥等[15] 對泵站同時(shí)做了模型試驗(yàn)和數(shù)值計(jì)算,分析了引水構(gòu)筑物內(nèi)水流流態(tài)及其配水特性,認(rèn)為采用分流墩、橫梁和張角漸擴(kuò)導(dǎo)流墩的組合式整流措施對分水箱涵、擴(kuò)散段、攔污柵處的流態(tài)有明顯優(yōu)化效果。

    本研究以湖北省某泵站進(jìn)水流道的優(yōu)化設(shè)計(jì)為目標(biāo),通過改變肘形流道的彎曲段長度、底邊線傾角和喉管高度,以及肘形流道的邊線傾角和喇叭口高度,擬定優(yōu)化設(shè)計(jì)方案,采用標(biāo)準(zhǔn)k-ε 模型和多面體網(wǎng)格建立數(shù)學(xué)模型,運(yùn)用Fluent 軟件求解,分析流道內(nèi)部流場分布,并通過比較流道出口截面的流速分布均勻度、速度加權(quán)平均角度和流道水力損失等水力指標(biāo),優(yōu)選出最佳方案,以期為合理設(shè)計(jì)泵站進(jìn)水流道提供參考依據(jù)。

    1 數(shù)值計(jì)算模型

    1.1 控制方程與湍流模型

    流體流動(dòng)的規(guī)律必須滿足質(zhì)量守恒、動(dòng)量守恒和能量守恒3 大定理,將流體的特點(diǎn)結(jié)合到3 大定理中,可以得到流體流動(dòng)的基本控制方程,本研究的CFD 流場計(jì)算針對不可壓縮黏性流動(dòng),所以只涉及質(zhì)量守恒、動(dòng)量守恒。

    (1)質(zhì)量守恒方程。連續(xù)方程是質(zhì)量守恒定理在流體流動(dòng)中應(yīng)用。針對流場的任何一個(gè)控制體,流體在流動(dòng)的過程中,源源不斷地出入控制面,控制體中流體的質(zhì)量可能隨時(shí)間的變化而變化,但是,控制體自身不會(huì)產(chǎn)生質(zhì)量,也不可能消除質(zhì)量。也就是說,在單位時(shí)間內(nèi),出入控制面流體的總質(zhì)量等于控制體中流體的質(zhì)量增量。

    式中 ρ——流體密度

    ui——流體i 方向運(yùn)動(dòng)速度

    t——時(shí)間

    xi——三維笛卡爾坐標(biāo)對象i 方向距離

    (2)動(dòng)量守恒方程。運(yùn)動(dòng)方程是動(dòng)量守恒定理在流體流動(dòng)中的應(yīng)用。針對流場的任何一個(gè)控制體,其控制面受到外界環(huán)境施加的表面力及質(zhì)量力場施加的質(zhì)量力。單位時(shí)間控制體受到的質(zhì)量力與表面力的合力等于控制體的動(dòng)量變化。運(yùn)動(dòng)方程的微分形式為

    式中 fi——單位質(zhì)量流體所受質(zhì)量力

    σij——流體應(yīng)力張量

    uj——流體j 方向運(yùn)動(dòng)速度

    xj——三維笛卡爾坐標(biāo)對象j 方向距離

    (3)湍流模型。RNG k-ε 湍流模型通過引入兩個(gè)附加方程(k 方程和ε 方程)來計(jì)算出湍流黏性系數(shù),進(jìn)一步計(jì)算出雷諾應(yīng)力。其中,k 方程是湍動(dòng)能方程,ε 方程是湍流耗散率方程。

    式中 k、ε——湍動(dòng)能和湍動(dòng)耗散率

    μeff——湍流渦黏系數(shù)

    ak、aε——k 和ε 有效普朗特?cái)?shù)倒數(shù)

    Gb——浮力引起的湍動(dòng)能k 產(chǎn)生項(xiàng)

    YM——可壓湍流脈動(dòng)擴(kuò)張項(xiàng)

    Gk——速度梯度引起應(yīng)力源項(xiàng)

    C1ε、C2ε、C3ε——經(jīng)驗(yàn)常數(shù)項(xiàng)

    1.2 邊界條件與網(wǎng)格劃分

    本研究對進(jìn)水流道、葉輪、導(dǎo)葉和出水流道采用計(jì)算精度更高的多面體網(wǎng)格進(jìn)行離散化處理。在保證計(jì)算精度和計(jì)算效率的前提下,進(jìn)行網(wǎng)格尺寸控制。其中,對于結(jié)構(gòu)較為規(guī)則的進(jìn)水流道和出水流道部分,最大網(wǎng)格尺寸控制100 mm,對于結(jié)構(gòu)較為復(fù)雜的葉輪和導(dǎo)葉部分,最大網(wǎng)格尺寸控制60 mm。最終確定下來的不同方案的網(wǎng)格數(shù)量控制250 萬,網(wǎng)格質(zhì)量在element quality 評判標(biāo)準(zhǔn)下達(dá)到0.4 以上,滿足計(jì)算要求,網(wǎng)格示意如圖1 所示。

    進(jìn)口邊界設(shè)置為質(zhì)量流量進(jìn)口(mass-flow-inlet),Q=34 938 kg/s(35 m3/s);出口邊界設(shè)置為自由出流(outflow);進(jìn)?出水流道及導(dǎo)葉區(qū)域的壁面采用靜止無滑移邊界, 軸流泵葉輪所在的區(qū)域通過MovingMesh 模型進(jìn)行建模,葉輪轉(zhuǎn)速125 r/min。操作條件:本研究所有CFD 數(shù)值模擬均認(rèn)為泵站在設(shè)計(jì)揚(yáng)程下運(yùn)行,泵站前池水位高程22.15 m,流道進(jìn)口操作壓力等效給定118 325 Pa(絕對壓力)。另外,考慮重力對水流的做功,浮力項(xiàng)操作密度設(shè)定為0。

    速度壓力耦合方法采用Coupled 方法,壓力差分格式采用PRESTO 格式,動(dòng)量、湍動(dòng)能及湍流耗散率的差分格式均采用二階迎風(fēng)格式(second order upwind)。其中二階差分是高階求解格式,其特點(diǎn)是求解準(zhǔn)確,結(jié)果可信性高,另外迎風(fēng)格式收斂性好。

    2 水力性能評價(jià)指標(biāo)

    進(jìn)水流道的出口斷面實(shí)際上是水泵葉輪室的進(jìn)口斷面,水泵葉輪進(jìn)口斷面的水力特性直接影響到水泵的性能,為葉輪室提供均勻的速度場和壓力場是進(jìn)水流道設(shè)計(jì)最為重要的目標(biāo)。因此,對進(jìn)水流道水力優(yōu)化設(shè)計(jì),應(yīng)結(jié)合對進(jìn)水流道的水力設(shè)計(jì)要求,確定科學(xué)有效合理的優(yōu)選準(zhǔn)則,否則對進(jìn)水流道的優(yōu)化水力設(shè)計(jì)沒有實(shí)際意義和應(yīng)用價(jià)值。一般而言,進(jìn)水流道的設(shè)計(jì)質(zhì)量通過以下3 個(gè)方面反映。

    2.1 葉輪進(jìn)口軸向流速分布均勻度Vu

    速度均勻度是反應(yīng)截面速度分布是否均勻的一種量化指標(biāo),其數(shù)學(xué)意義與速度的標(biāo)準(zhǔn)差和速度加權(quán)平均值的比值相關(guān)。

    式中 uˉa——某特定截面平均軸向速度

    uai——斷面第i 個(gè)節(jié)點(diǎn)軸向速度

    m——截面網(wǎng)格節(jié)點(diǎn)個(gè)數(shù)

    2.2 速度加權(quán)平均角度θ

    流道過流斷面速度加權(quán)平均角度是評價(jià)流道的另一個(gè)指標(biāo)。對于進(jìn)口流道而言,出口斷面的加權(quán)平均角度越接近90°,則進(jìn)水流道內(nèi)部的水流轉(zhuǎn)向越好,引入葉輪的水流也滿足設(shè)計(jì)要求。計(jì)算公式為

    式中 uti——出口截面第i 個(gè)網(wǎng)格節(jié)點(diǎn)切向速度

    2.3 流道水力損失Δh

    水流的黏滯性會(huì)使水流內(nèi)部產(chǎn)生摩擦力,以此來抵抗水流之間的相互作用,從而產(chǎn)生水力損失,通常采用伯努利方程來計(jì)算流道水力損失。計(jì)算公式為

    式中 P1、P2——流道進(jìn)出口斷面靜壓

    Z1、Z2——流道進(jìn)出口斷面高程

    u1、u2——流道進(jìn)出口斷面速度

    3 肘形進(jìn)水流道優(yōu)化計(jì)算

    由于肘形流道涉及的幾何參數(shù)眾多,在肘形流道的優(yōu)化設(shè)計(jì)工作中,部分結(jié)構(gòu)尺寸將按照提供的圖紙《泵址泵房結(jié)構(gòu)布置圖》進(jìn)行約束,包括進(jìn)口高度、寬度、流道長度和出口直徑。流道進(jìn)口底部高程(13.3 m)和出口高程(16.1 m、與喇叭口銜接處)保持不變,參與優(yōu)化的結(jié)構(gòu)參數(shù)包括彎曲段長度(W=L?Z)、底邊線傾角和喉管高度。肘形流道結(jié)構(gòu)參數(shù)如表1 所示,肘形流道結(jié)構(gòu)如圖2 所示。

    采用CFD 數(shù)值建模方法(考慮軸流泵的旋轉(zhuǎn)流動(dòng)),對以下6 種設(shè)計(jì)方案進(jìn)行流體動(dòng)力學(xué)分析,以驗(yàn)證優(yōu)化方案的可行性,進(jìn)而對優(yōu)化方案提出更合理的改進(jìn)依據(jù)。

    方案1:彎曲段長度W=3 040 mm、喉管高度H=2 680 mm、底邊線傾角=7°;根據(jù)CAD 圖紙《泵址泵房結(jié)構(gòu)布置圖》確定;導(dǎo)流隔墩參數(shù)與方案3 相同。

    方案2:彎曲段長度W=3 040 mm、喉管高度H=2 857 mm、底邊線傾角4°;通過曲面響應(yīng)法計(jì)算獲?。桓舳臻L度7 000 mm、首端直徑1 200 mm、末端直徑600 mm。

    方案3:導(dǎo)流隔墩長度延長至9 000 mm,其余參數(shù)均與方案2 相同。

    方案4:導(dǎo)流隔墩長度延長至11 000 mm,其余參數(shù)均與方案2 相同。

    方案5: 流道寬度逐漸縮小, 彎曲段初始寬度7 000 mm,其余參數(shù)均與方案3 相同。

    方案6:底邊線傾角設(shè)置0 °,其余參數(shù)均與方案3 相同。

    3.1 流道出口截面流態(tài)

    肘形進(jìn)水流道出口截面軸向速度分布對比如圖3所示,方案1~方案6 的最大軸向速度依次為5.4、5.8、5.8、5.9、5.9 和5.8 m/s,均位于靠近流道進(jìn)口的一側(cè);從軸向速度分布形態(tài)來看,方案1 的軸向速度表現(xiàn)出了明顯的不對稱性特征,而方案6 的軸向速度對稱效果最好;從方案3~方案5 的結(jié)果來看,不同的導(dǎo)流隔墩長度對出口截面的軸向速度對稱效果的影響并不明顯??傮w上來說,肘形流道出口截面的軸向速度分布并不是均勻分布的,其速度變化3.0~6.0 m/s。

    不同方案出口截面切向速度如圖4 所示,從分布形態(tài)來看,切向速度關(guān)于軸流泵旋轉(zhuǎn)中心呈現(xiàn)正負(fù)值對稱的分布規(guī)律,這是水流在外側(cè)壁面的約束下向中心匯聚的結(jié)果,亦可以表明,軸流泵的旋轉(zhuǎn)并不會(huì)引起流道出口的水流發(fā)生旋轉(zhuǎn)。出口的切向速度范圍在±1.5 m/s,從空間速度矢量的角度來考慮,切向速度峰值的分布區(qū)域越靠近流道進(jìn)口的一側(cè),則意味著水流在彎曲段的引導(dǎo)下更快的轉(zhuǎn)向,因此方案5 的效果最好,方案1 的效果最差。

    流道出口截面湍動(dòng)能對比結(jié)果如圖5 所示,從分布形態(tài)來看,出口截面的湍動(dòng)能呈現(xiàn)對稱分布規(guī)律,靠近流道進(jìn)口一側(cè)的湍動(dòng)能較高,而靠近出口一側(cè)的湍動(dòng)能則處于較低水平??傮w上來說,肘形進(jìn)水流道出口截面的湍動(dòng)能變化0.1~0.4 m2/s2。

    3.2 水力性能評價(jià)指標(biāo)

    在肘形進(jìn)水流道的水力指標(biāo)計(jì)算中,所有CFD 算例均包括軸流泵的旋轉(zhuǎn)流場,以及泵后導(dǎo)葉的流場。由于缺乏軸流泵的葉型圖紙,數(shù)值模擬中選用的軸流泵模型通過CF-turbo 軟件進(jìn)行設(shè)計(jì),通過初步驗(yàn)證,當(dāng)轉(zhuǎn)速125 r/min、流量35 m3/s 時(shí),模擬中采用的軸流泵揚(yáng)程5.2~5.8 m,基本滿足實(shí)際情況。

    由表2 可知,方案1 的出口速度分布均勻度最高,Vu=88.81%,而曲面響應(yīng)優(yōu)化獲取的方案3 則比方案1低了0.11%;從方案2~方案4 的對比來看,過長或過短的導(dǎo)流隔墩都將導(dǎo)致速度分布均勻度降低,合理的隔墩長度應(yīng)在9 000 mm 左右。另外,方案5 和方案6的結(jié)果顯示,繼續(xù)減小底邊線傾角或采用逐漸縮小的直管段對速度均勻度的優(yōu)化效果都不理想。

    出口速度加權(quán)平均角θ 計(jì)算結(jié)果如表3 所示,總體上來看,所有方案的θ 均保持在88.5 °以上的水平,不同方案之間的差異并不明顯。對比結(jié)果表明,方案3中的水流轉(zhuǎn)向效果最好(θ 最接近90 °),而方案6 的水流轉(zhuǎn)向效果最差。

    水力損失(Δh)特性對比如表4 所示(在Δh 的計(jì)算中,流道進(jìn)口的高程取整個(gè)截面高程的加權(quán)平均值,進(jìn)?出口高程差0.3 m)。對比結(jié)果顯示,方案2 的水力損失最大,Δh=0.265 m,k=2.163e-04;方案5 的水力損失最小,Δh=0.226 m,k=1.845e-04;總體上來說,6種肘形進(jìn)水流道的水力損失差異在0.04 m 以內(nèi),故在流道設(shè)計(jì)中可將該指標(biāo)的優(yōu)先級降低。

    由表5 可知,方案1、方案3 和方案4 出口截面的湍動(dòng)能差異并不顯著,而方案5 和方案6 的湍動(dòng)能則明顯提高,該結(jié)果表明,減小底邊線的傾角,或采用截面積逐漸縮小的直管段,會(huì)加劇流道出口的不穩(wěn)定流動(dòng)。

    3.3 流道全局流動(dòng)特性

    由圖6 可知,不同方案的流速分布規(guī)律基本一致,隨著水流發(fā)展的方向,流速逐漸增大,在水泵吸入口達(dá)到最高值。需要注意的是,由于導(dǎo)流隔墩的存在,彎曲段內(nèi)的速度變化梯度較大,水流在剪切作用下可能會(huì)產(chǎn)生大尺度旋渦,因此隔墩長度不宜過長。

    流場湍動(dòng)能分布如圖7 所示,為了體現(xiàn)水流由前池到達(dá)流道進(jìn)口時(shí)存在大尺度的渦運(yùn)動(dòng),在CFD 數(shù)值模型中,對肘形流道進(jìn)口給定了較高的湍動(dòng)能初值。計(jì)算結(jié)果顯示,流場的湍動(dòng)能沿流動(dòng)發(fā)展方向迅速減小,在喉管附近達(dá)到最低值;到達(dá)彎曲段后,湍動(dòng)能在軸流泵的作用下又逐漸開始增大。總體上來看,5 種方案的肘形流道均對水流具有良好的整流效果,相比較而言,方案1、方案3 和方案4 在喇叭口區(qū)域內(nèi)的湍流強(qiáng)度要略低于其他3 種方案,這一結(jié)論與表5 中,流道出口湍動(dòng)能的平均計(jì)算結(jié)果相互對應(yīng)。

    3.4 沿流動(dòng)方向速度變化

    沿流動(dòng)方向各過流斷面速度變化曲線如圖8 所示,各方案沿流道方向過流斷面速度單調(diào)增加,曲線基本保持光滑,說明了在流量不變的情況下,過流斷面面積保持均勻減小,流道是逐漸光滑過渡的收縮流道,不會(huì)出現(xiàn)突變的流道,從而保證了水力損失較小,還保證了水泵進(jìn)口流體具有一定的速度。

    4 結(jié)束語

    本研究以湖北省某泵站進(jìn)水流道的優(yōu)化設(shè)計(jì)為目標(biāo),通過改變肘形流道的彎曲段長度、底邊線傾角和喉管高度,擬定優(yōu)化設(shè)計(jì)方案,采用標(biāo)準(zhǔn)k-ε 紊流模型和多面體網(wǎng)格模型,運(yùn)用Fluent 軟件計(jì)算分析流道內(nèi)部流場分布、水力性能和沿流動(dòng)方向速度變化。

    對于肘形進(jìn)水流道,方案1 和方案3 均具有良好的水力性能, 對比發(fā)現(xiàn), 方案1( Vu=88.81%,Ek=0.116 m2/s2)出口的流速分布均勻度和湍動(dòng)能優(yōu)于方案3( Vu=88.70%, Ek=0.119 m2/s2) , 而方案3(θ=88.83°,Δh=0.249 mm)的出口速度加權(quán)平均角和水力損失優(yōu)于方案1(θ=88.77°,Δh=0.255 mm)。總體上來看,兩種方案的水力性能指標(biāo)差異不大,還需根據(jù)進(jìn)?出水流道全局流場的CFD 計(jì)算結(jié)果來確定優(yōu)選方案。

    參考文獻(xiàn)

    [1]郁片紅,李君菡,付小莉.大型雨水泵站進(jìn)水流道水力特性研究[J].中國農(nóng)村水利水電,2020(10):126-130.

    YU Pianhong,LI Junhan,F(xiàn)U Xiaoli.Research on the hydraulic characteristicsof the inlet flow-passage of large rainwater pumping stations[J].China Rural Water and Hydropower,2020(10):126-130.

    [2]盧春暉,張永麗,唐劍韜.泵站前池水力性能的數(shù)值模擬[J].流體機(jī)械,2012,40(7):20-23,28.

    LU Chunhui, ZHANG Yongli, TANG Jiantao. Numerical simulationof hydraulic performance in the"forebay of pumping station[J]. FluidMachinery,2012,40(7):20-23,28.

    [3]江文,于建忠,傅宗甫,等.彎道河段下游泵站進(jìn)水前池流態(tài)及整流措施[J].人民黃河,2019,41(4):84-87.

    JIANG Wen, YU Jianzhong, FU Zongfu, et al. Flow patterns andrectification measures in a forebay of pumping station designed at bendriver[J].Yellow River,2019,41(4):84-87.

    [4]蘇正洋, 陸嘉偉, 張志韜, 等. 泵站前池整流技術(shù)研究綜述[J].人民珠江,2020,41(1):69-75.

    SU Zhengyang,LU Jiawei,ZHANG Zhitao,et al.Study on rectifiertechnology in forebay of pumping station[J]. Pearl River, 2020,41(1):69-75.

    [5] 水利部.泵站設(shè)計(jì)規(guī)范[M].北京:中國計(jì)劃出版社,1997.

    [6] 黃繼宏,王曉升,張曉毅.城市排水泵站前池導(dǎo)流墩整流措施的數(shù)值模擬[J].江蘇農(nóng)業(yè)科學(xué),2013,41(5):360-362.

    [7]羅海軍,張睿,徐輝.改善城市排水泵站進(jìn)水流態(tài)的試驗(yàn)研究[J].中國農(nóng)村水利水電,2019(1):176-179.

    LUO Haijun, ZHANG Rui, XU Hui. Experimental research on theimprovement in the inflow pattern of an urban drainage pumpingstation[J].China Rural Water and Hydropower,2019(1):176-179.

    [8]王翌.基于數(shù)值計(jì)算和模型試驗(yàn)的泵站側(cè)向進(jìn)水結(jié)構(gòu)優(yōu)化[J].水利科學(xué)與寒區(qū)工程,2023,6(11):73-76.

    WANG Yi. Optimization of lateral inlet structure of pumping stationsbased on numerical calculation and model experiments[J]. Hydro Scienceand Cold Zone Engineering,2023,6(11):73-76.

    [9]奚斌,陳志剛,陳葉欣,等.復(fù)雜進(jìn)流條件下泵站進(jìn)水流態(tài)優(yōu)化試驗(yàn)[J].揚(yáng)州大學(xué)學(xué)報(bào)(自然科學(xué)版),2023,26(6):73-78.

    XI Bin,CHEN Zhigang,CHEN Yexin,et al.Experimental study onoptimization of inlet flow pattern in pumping station under complex inletconditions[J]. Journal of Yangzhou University ( Natural Science Edition),2023,26(6):73-78.

    [10]王為術(shù),甄娟,郭嘉偉,等.并聯(lián)循環(huán)泵站前池及吸水室整流優(yōu)化[J].排灌機(jī)械工程學(xué)報(bào),2023,41(11):1119-1126.

    WANG Weishu, ZHEN Juan, GUO Jiawei, et al. Optimization ofrectification for forebay and suction chamber of parallel circulation pumpstation[J]. Journal of Drainage and Irrigation Machinery Engineering,2023,41(11):1119-1126.

    [11]胥維纖,張曉斐,孫斌.基于模型試驗(yàn)的泵站進(jìn)水前池優(yōu)化和特性研究[J].西北水電,2023(3):82-87.

    XU Weixian,ZHANG Xiaofei,SUN Bin.Study on hydraulic characteristicsof pumping station intake forebay undershape optimization based on hydraulic model test[J]. Northwest Hydropower, 2023( 3) : 82-87.

    [12]施偉,倪春,于賢磊,等.閘站并列布置工程泵站前池流態(tài)改善[J].南水北調(diào)與水利科技,2023,21(4):810-820.

    SHI Wei,NI Chun,YU Xianlei,et al.Flow pattern improvement inthe forebay of pumping station with parallel sluice[J]. South-to-NorthWater Transfers and Water Science amp; Technology,2023,21(4):810-820.

    [13]資丹,王福軍,姚志峰,等.大型泵站進(jìn)水流場組合式導(dǎo)流墩整流效果分析[J].農(nóng)業(yè)工程學(xué)報(bào),2015,31(16):71-77.

    ZI Dan,WANG Fujun,YAO Zhiwei,et al.Effects analysis on rectifyingintake flow field for large scale pumping station with combined diversionpiers[J]. Transactions of the Chinese Society of AgriculturalEngineering,2015,31(16):71-77.

    [14]張睿,徐輝,陳毓陵,等.斜向管涵進(jìn)流城市雨水泵站箱涵流態(tài)分析及整流措施研究[J].水利學(xué)報(bào),2018,49(5):598-607.

    ZHANG Rui,XU Hui,CHEN Yuling,et al.Study on the flow patternsand rectification measures of box culvert of urban storm waterpumping station with oblique pipe culvert[J].Journal of Hydraulic Engineering,2018,49(5):598-607.

    [15]李志祥,馮建剛,錢尚拓,等.斜向進(jìn)流城市雨水泵站引水建筑物流態(tài)改善數(shù)值模擬[J].中國農(nóng)村水利水電,2020(6):139-143,147.

    LI Zhixiang, FENG Jiangang, QIAN Shangtuo, et al. Numericalsimulation of improvement of flow pattern in the diversion structure ofthe urban storm water pumping station with slanted intake[J]. ChinaRural Water and Hydropower,2020(6):139-143,147.

    基金項(xiàng)目: 湖北省重點(diǎn)研發(fā)計(jì)劃項(xiàng)目(2023BBB085)

    猜你喜歡
    計(jì)算流體力學(xué)數(shù)值模擬泵站
    張家邊涌泵站建設(shè)難點(diǎn)及技術(shù)創(chuàng)新實(shí)踐
    2016年河南省己建成泵站數(shù)量
    《工程流體力學(xué)》教學(xué)方法探討
    基于預(yù)條件技術(shù)的風(fēng)力機(jī)葉片計(jì)算方法研究
    跨音速飛行中機(jī)翼水汽凝結(jié)的數(shù)值模擬研究
    科技視界(2016年18期)2016-11-03 20:38:17
    姚橋煤礦采空區(qū)CO2防滅火的數(shù)值模擬分析
    雙螺桿膨脹機(jī)的流場數(shù)值模擬研究
    科技視界(2016年22期)2016-10-18 14:53:19
    一種基于液壓緩沖的減震管卡設(shè)計(jì)與性能分析
    科技視界(2016年20期)2016-09-29 11:08:27
    全省已建成泵站數(shù)量
    民用飛機(jī)靜壓源位置誤差修正設(shè)計(jì)研究①
    科技資訊(2015年17期)2015-10-09 21:02:59
    av在线播放免费不卡| 亚洲精品一二三| 欧美激情 高清一区二区三区| 色综合站精品国产| 精品国产乱码久久久久久男人| 国产精华一区二区三区| 成人影院久久| 黑人巨大精品欧美一区二区mp4| 18禁观看日本| 亚洲熟妇中文字幕五十中出 | 免费日韩欧美在线观看| 亚洲国产精品999在线| 久久香蕉国产精品| 久9热在线精品视频| 91大片在线观看| 亚洲午夜精品一区,二区,三区| 女人高潮潮喷娇喘18禁视频| 久久久久久久午夜电影 | 黑人巨大精品欧美一区二区蜜桃| av有码第一页| 欧美乱妇无乱码| 免费久久久久久久精品成人欧美视频| 一本综合久久免费| 亚洲av成人不卡在线观看播放网| 亚洲狠狠婷婷综合久久图片| aaaaa片日本免费| 日日夜夜操网爽| 亚洲精品在线美女| 午夜福利,免费看| 免费看十八禁软件| 成人18禁高潮啪啪吃奶动态图| 亚洲午夜精品一区,二区,三区| 美女大奶头视频| 黄色视频,在线免费观看| 久久久国产成人免费| 国产在线精品亚洲第一网站| 亚洲精品中文字幕在线视频| 国产午夜精品久久久久久| 中文欧美无线码| 黄色丝袜av网址大全| 又黄又爽又免费观看的视频| 欧美黑人精品巨大| 又大又爽又粗| 最新美女视频免费是黄的| 亚洲五月色婷婷综合| 日本 av在线| 18禁美女被吸乳视频| 亚洲成a人片在线一区二区| 亚洲精品av麻豆狂野| 亚洲精品中文字幕在线视频| 久久久久久大精品| 日韩一卡2卡3卡4卡2021年| 麻豆一二三区av精品| 麻豆久久精品国产亚洲av | 99国产精品99久久久久| 97碰自拍视频| 日韩视频一区二区在线观看| 国产区一区二久久| 好看av亚洲va欧美ⅴa在| 男人舔女人的私密视频| 欧美日韩中文字幕国产精品一区二区三区 | 午夜影院日韩av| 淫秽高清视频在线观看| 国产欧美日韩综合在线一区二区| 在线国产一区二区在线| 波多野结衣av一区二区av| 91成年电影在线观看| 中文欧美无线码| 最新美女视频免费是黄的| 久久国产精品影院| 88av欧美| 99精品久久久久人妻精品| 亚洲情色 制服丝袜| 亚洲精品成人av观看孕妇| av在线播放免费不卡| 亚洲国产精品一区二区三区在线| 国产激情久久老熟女| a级毛片黄视频| 99re在线观看精品视频| 亚洲国产欧美一区二区综合| 国产高清国产精品国产三级| 不卡一级毛片| 黄色视频不卡| 日韩有码中文字幕| 午夜日韩欧美国产| 亚洲国产精品sss在线观看 | 丝袜美足系列| 搡老岳熟女国产| 午夜亚洲福利在线播放| svipshipincom国产片| av免费在线观看网站| 国产av一区二区精品久久| 一级片免费观看大全| 校园春色视频在线观看| 欧美成人午夜精品| 一二三四社区在线视频社区8| 亚洲色图 男人天堂 中文字幕| 一级毛片女人18水好多| 纯流量卡能插随身wifi吗| 日韩三级视频一区二区三区| 国产精品美女特级片免费视频播放器 | 91老司机精品| 好男人电影高清在线观看| 午夜福利免费观看在线| 国产区一区二久久| 九色亚洲精品在线播放| 日本黄色日本黄色录像| e午夜精品久久久久久久| 亚洲专区字幕在线| 亚洲激情在线av| 日韩大码丰满熟妇| 看免费av毛片| 女人高潮潮喷娇喘18禁视频| 一级片免费观看大全| 最新美女视频免费是黄的| 久久久久久免费高清国产稀缺| 亚洲熟妇中文字幕五十中出 | 久久久国产欧美日韩av| 国产精品 国内视频| 久久午夜综合久久蜜桃| 久久久久亚洲av毛片大全| 久久中文字幕一级| 成人影院久久| 国产精品影院久久| 老鸭窝网址在线观看| 日本三级黄在线观看| 亚洲欧美日韩另类电影网站| 久久久国产精品麻豆| 夜夜夜夜夜久久久久| 欧美精品啪啪一区二区三区| 亚洲久久久国产精品| 日韩欧美一区二区三区在线观看| 99在线视频只有这里精品首页| 精品卡一卡二卡四卡免费| 在线视频色国产色| 18禁裸乳无遮挡免费网站照片 | 在线观看免费日韩欧美大片| 97超级碰碰碰精品色视频在线观看| 色尼玛亚洲综合影院| 在线观看免费日韩欧美大片| 在线观看免费午夜福利视频| 成年人黄色毛片网站| 十分钟在线观看高清视频www| 人妻丰满熟妇av一区二区三区| 久久国产精品影院| 人人妻,人人澡人人爽秒播| 日韩精品青青久久久久久| 国产av一区二区精品久久| 欧美黄色片欧美黄色片| 亚洲一区二区三区欧美精品| 如日韩欧美国产精品一区二区三区| 最新美女视频免费是黄的| 中文字幕精品免费在线观看视频| 欧美一区二区精品小视频在线| 国产蜜桃级精品一区二区三区| 亚洲性夜色夜夜综合| 操出白浆在线播放| 国产aⅴ精品一区二区三区波| 777久久人妻少妇嫩草av网站| 在线观看一区二区三区| 中亚洲国语对白在线视频| www.熟女人妻精品国产| av免费在线观看网站| 亚洲国产精品999在线| 精品免费久久久久久久清纯| 精品一区二区三区四区五区乱码| 国产精品一区二区免费欧美| 免费在线观看影片大全网站| 天堂动漫精品| 18禁黄网站禁片午夜丰满| 亚洲成国产人片在线观看| 老熟妇乱子伦视频在线观看| 人人妻,人人澡人人爽秒播| 国产三级黄色录像| av福利片在线| 亚洲情色 制服丝袜| 一级黄色大片毛片| 欧美激情 高清一区二区三区| 大香蕉久久成人网| 中出人妻视频一区二区| 国产99久久九九免费精品| 国产精品99久久99久久久不卡| 色哟哟哟哟哟哟| 亚洲精品国产精品久久久不卡| 欧美日韩视频精品一区| 国产精品电影一区二区三区| 国产无遮挡羞羞视频在线观看| 男女做爰动态图高潮gif福利片 | 俄罗斯特黄特色一大片| 亚洲一区高清亚洲精品| 国产成人啪精品午夜网站| 亚洲欧美日韩另类电影网站| avwww免费| 亚洲五月天丁香| 国产97色在线日韩免费| 免费观看精品视频网站| 久久婷婷成人综合色麻豆| 99热只有精品国产| 视频在线观看一区二区三区| 色综合欧美亚洲国产小说| 99国产精品一区二区三区| 亚洲国产精品sss在线观看 | 大码成人一级视频| 无限看片的www在线观看| av中文乱码字幕在线| 露出奶头的视频| 久久这里只有精品19| 精品卡一卡二卡四卡免费| 人人妻,人人澡人人爽秒播| 成人免费观看视频高清| 久久久国产成人精品二区 | 在线免费观看的www视频| 国产av精品麻豆| 欧美乱妇无乱码| 黑人猛操日本美女一级片| 高清欧美精品videossex| 精品久久久久久,| 亚洲九九香蕉| 十分钟在线观看高清视频www| 91在线观看av| 免费看a级黄色片| 亚洲人成77777在线视频| 黄色片一级片一级黄色片| 99国产精品99久久久久| 嫩草影视91久久| 琪琪午夜伦伦电影理论片6080| 亚洲人成电影免费在线| 老司机亚洲免费影院| 日本黄色日本黄色录像| 国产片内射在线| 亚洲伊人色综图| 在线观看午夜福利视频| 日韩国内少妇激情av| 88av欧美| 精品午夜福利视频在线观看一区| 欧美+亚洲+日韩+国产| a级毛片在线看网站| 久久久久久久久中文| 两个人免费观看高清视频| 无人区码免费观看不卡| 身体一侧抽搐| 国产成人影院久久av| 女人爽到高潮嗷嗷叫在线视频| 黄片小视频在线播放| 欧美中文日本在线观看视频| 自线自在国产av| 18美女黄网站色大片免费观看| 精品高清国产在线一区| 男女床上黄色一级片免费看| 国内久久婷婷六月综合欲色啪| 国产av在哪里看| 成人手机av| 不卡av一区二区三区| 长腿黑丝高跟| 级片在线观看| 日本五十路高清| 国产色视频综合| 99国产精品一区二区三区| 日本黄色日本黄色录像| 成人18禁高潮啪啪吃奶动态图| 在线视频色国产色| 99久久国产精品久久久| 99久久久亚洲精品蜜臀av| 在线观看一区二区三区激情| 久久香蕉激情| 999久久久精品免费观看国产| 亚洲一码二码三码区别大吗| 欧美 亚洲 国产 日韩一| 国产黄a三级三级三级人| 黑人猛操日本美女一级片| 午夜影院日韩av| 午夜精品久久久久久毛片777| 成年版毛片免费区| 亚洲欧美一区二区三区黑人| 在线观看舔阴道视频| 人人妻人人爽人人添夜夜欢视频| 国产片内射在线| 亚洲国产精品999在线| 欧美日韩福利视频一区二区| 一级a爱视频在线免费观看| 日韩欧美免费精品| 黑人操中国人逼视频| 日韩欧美三级三区| 天堂影院成人在线观看| 日本 av在线| 欧美丝袜亚洲另类 | a级片在线免费高清观看视频| 国产单亲对白刺激| 91成人精品电影| 一区在线观看完整版| 亚洲专区国产一区二区| 精品福利永久在线观看| 黑人操中国人逼视频| 亚洲成av片中文字幕在线观看| 欧美日韩亚洲国产一区二区在线观看| 国产男靠女视频免费网站| 久久午夜亚洲精品久久| 亚洲精品一卡2卡三卡4卡5卡| 免费观看人在逋| 最新美女视频免费是黄的| 国产精品久久久av美女十八| 免费不卡黄色视频| www.自偷自拍.com| av福利片在线| 午夜视频精品福利| 一区二区三区精品91| 老司机亚洲免费影院| 视频区欧美日本亚洲| 精品第一国产精品| 久久精品aⅴ一区二区三区四区| 免费不卡黄色视频| 波多野结衣高清无吗| 如日韩欧美国产精品一区二区三区| 亚洲av五月六月丁香网| 精品免费久久久久久久清纯| 国产亚洲精品一区二区www| 欧美日本中文国产一区发布| 人人澡人人妻人| 婷婷六月久久综合丁香| 婷婷丁香在线五月| 日本欧美视频一区| 亚洲成国产人片在线观看| 日韩精品青青久久久久久| 欧美成狂野欧美在线观看| 老熟妇乱子伦视频在线观看| 日日干狠狠操夜夜爽| 亚洲国产欧美网| 国产成年人精品一区二区 | 99在线人妻在线中文字幕| 国产人伦9x9x在线观看| 亚洲精品一二三| 精品久久久久久成人av| 久久狼人影院| 成人亚洲精品av一区二区 | 制服人妻中文乱码| 另类亚洲欧美激情| 最近最新中文字幕大全电影3 | 在线观看66精品国产| 久久精品国产亚洲av高清一级| 女人爽到高潮嗷嗷叫在线视频| 淫秽高清视频在线观看| 99久久99久久久精品蜜桃| 成人免费观看视频高清| 看黄色毛片网站| 午夜福利,免费看| 正在播放国产对白刺激| 久久中文看片网| 国产高清videossex| 一区福利在线观看| 免费久久久久久久精品成人欧美视频| 美女午夜性视频免费| 精品久久蜜臀av无| 亚洲国产看品久久| 亚洲国产中文字幕在线视频| 国产免费av片在线观看野外av| 成熟少妇高潮喷水视频| 欧美黑人精品巨大| 久久久久精品国产欧美久久久| 精品一品国产午夜福利视频| 色尼玛亚洲综合影院| 亚洲片人在线观看| 高清在线国产一区| 亚洲成人国产一区在线观看| 国产成人av激情在线播放| 久久中文看片网| 正在播放国产对白刺激| 88av欧美| 窝窝影院91人妻| 亚洲精品国产精品久久久不卡| 亚洲 欧美一区二区三区| 长腿黑丝高跟| 在线观看免费高清a一片| 亚洲av日韩精品久久久久久密| 国产激情欧美一区二区| 国产免费现黄频在线看| www国产在线视频色| 成人特级黄色片久久久久久久| 国产一卡二卡三卡精品| 久久午夜亚洲精品久久| 水蜜桃什么品种好| 亚洲成人国产一区在线观看| 性欧美人与动物交配| 99久久人妻综合| 日本一区二区免费在线视频| 亚洲 欧美 日韩 在线 免费| 亚洲国产欧美一区二区综合| 亚洲国产精品合色在线| 国产精品av久久久久免费| 真人做人爱边吃奶动态| 亚洲视频免费观看视频| 人妻久久中文字幕网| 欧美一区二区精品小视频在线| a级毛片在线看网站| 十八禁人妻一区二区| 757午夜福利合集在线观看| 啦啦啦 在线观看视频| 久久人妻av系列| 精品国产美女av久久久久小说| 亚洲精品国产一区二区精华液| 在线视频色国产色| 精品无人区乱码1区二区| 国产精品99久久99久久久不卡| 在线观看一区二区三区| 亚洲中文av在线| 国产国语露脸激情在线看| 国产精品一区二区在线不卡| 色哟哟哟哟哟哟| 久久中文字幕人妻熟女| 搡老乐熟女国产| 亚洲国产精品sss在线观看 | 久久久国产一区二区| 身体一侧抽搐| 极品教师在线免费播放| 91老司机精品| 精品人妻1区二区| 亚洲色图 男人天堂 中文字幕| 亚洲自偷自拍图片 自拍| 长腿黑丝高跟| 精品日产1卡2卡| 国产一区二区三区视频了| 一个人观看的视频www高清免费观看 | 国产免费男女视频| 国产男靠女视频免费网站| 天天添夜夜摸| 久久久久国产一级毛片高清牌| 亚洲aⅴ乱码一区二区在线播放 | 99久久综合精品五月天人人| 精品电影一区二区在线| 日本a在线网址| 丝袜美足系列| 黄色女人牲交| 亚洲人成电影免费在线| 日韩中文字幕欧美一区二区| 国产av精品麻豆| 丰满人妻熟妇乱又伦精品不卡| 村上凉子中文字幕在线| 1024视频免费在线观看| 欧美日韩国产mv在线观看视频| 免费高清视频大片| 欧美乱码精品一区二区三区| 国内久久婷婷六月综合欲色啪| 天堂俺去俺来也www色官网| 欧美成狂野欧美在线观看| 精品一区二区三区视频在线观看免费 | 很黄的视频免费| 999久久久精品免费观看国产| 女警被强在线播放| 国产深夜福利视频在线观看| 欧美一级毛片孕妇| 美女午夜性视频免费| 中文字幕最新亚洲高清| 亚洲男人天堂网一区| 国产片内射在线| 国产欧美日韩精品亚洲av| 免费观看精品视频网站| 久久人妻福利社区极品人妻图片| 色老头精品视频在线观看| 男女下面插进去视频免费观看| 国产在线精品亚洲第一网站| 亚洲av电影在线进入| 精品卡一卡二卡四卡免费| 满18在线观看网站| 国产精品99久久99久久久不卡| 成人手机av| 999久久久精品免费观看国产| 热99re8久久精品国产| 99久久人妻综合| 国产欧美日韩一区二区精品| 久久久国产欧美日韩av| 精品日产1卡2卡| 亚洲av美国av| 国产99久久九九免费精品| 丝袜美腿诱惑在线| 又黄又粗又硬又大视频| 丝袜美腿诱惑在线| 免费在线观看视频国产中文字幕亚洲| 亚洲狠狠婷婷综合久久图片| 国产乱人伦免费视频| 一级片免费观看大全| 老汉色∧v一级毛片| 在线观看免费午夜福利视频| 老熟妇仑乱视频hdxx| 视频区图区小说| 国产深夜福利视频在线观看| 国产麻豆69| 欧美+亚洲+日韩+国产| 一进一出抽搐动态| 国产av在哪里看| 国产精品 欧美亚洲| 日本精品一区二区三区蜜桃| 久久 成人 亚洲| 国产三级黄色录像| 老熟妇乱子伦视频在线观看| 丝袜人妻中文字幕| 99国产极品粉嫩在线观看| 性少妇av在线| 丰满的人妻完整版| 黄片播放在线免费| 国产精品av久久久久免费| av电影中文网址| 99riav亚洲国产免费| 天堂俺去俺来也www色官网| 在线观看免费午夜福利视频| 女人高潮潮喷娇喘18禁视频| 欧美成人免费av一区二区三区| 他把我摸到了高潮在线观看| 50天的宝宝边吃奶边哭怎么回事| 成人免费观看视频高清| 午夜日韩欧美国产| 电影成人av| 一级a爱视频在线免费观看| 久99久视频精品免费| 国产精品影院久久| 丰满人妻熟妇乱又伦精品不卡| 国产91精品成人一区二区三区| 国产精品永久免费网站| 视频区欧美日本亚洲| 日韩免费高清中文字幕av| 国产精品久久视频播放| 757午夜福利合集在线观看| 久久精品国产亚洲av高清一级| 岛国在线观看网站| 视频在线观看一区二区三区| 日韩精品青青久久久久久| 夜夜夜夜夜久久久久| 久久精品国产清高在天天线| 人人澡人人妻人| 91精品三级在线观看| 亚洲欧美日韩高清在线视频| 天堂√8在线中文| 热99国产精品久久久久久7| 国产精品久久久久久人妻精品电影| 悠悠久久av| 久久这里只有精品19| 欧美黑人精品巨大| 久久人妻av系列| 中文字幕另类日韩欧美亚洲嫩草| 亚洲午夜精品一区,二区,三区| 亚洲色图 男人天堂 中文字幕| ponron亚洲| 欧美中文综合在线视频| 亚洲久久久国产精品| 免费日韩欧美在线观看| 1024视频免费在线观看| 亚洲少妇的诱惑av| 免费久久久久久久精品成人欧美视频| 久久香蕉国产精品| 亚洲自偷自拍图片 自拍| 后天国语完整版免费观看| 999久久久国产精品视频| 欧洲精品卡2卡3卡4卡5卡区| 午夜福利,免费看| 欧美性长视频在线观看| 一个人观看的视频www高清免费观看 | √禁漫天堂资源中文www| 亚洲激情在线av| 一本综合久久免费| 亚洲精华国产精华精| 午夜免费成人在线视频| 午夜久久久在线观看| 国产精品爽爽va在线观看网站 | av福利片在线| 99国产极品粉嫩在线观看| 久久国产精品人妻蜜桃| 久久香蕉国产精品| 搡老熟女国产l中国老女人| 欧美精品一区二区免费开放| 欧美+亚洲+日韩+国产| 中文字幕最新亚洲高清| 欧美激情极品国产一区二区三区| 黄色成人免费大全| 国产成人av教育| av在线播放免费不卡| 最近最新中文字幕大全电影3 | 日本黄色视频三级网站网址| 性欧美人与动物交配| 在线观看免费视频网站a站| 亚洲av成人一区二区三| 一个人观看的视频www高清免费观看 | 成人亚洲精品av一区二区 | 一边摸一边抽搐一进一出视频| 免费在线观看黄色视频的| 一进一出好大好爽视频| 国产成+人综合+亚洲专区| а√天堂www在线а√下载| 午夜精品国产一区二区电影| 亚洲精品国产色婷婷电影| 日韩三级视频一区二区三区| 在线观看66精品国产| 国产欧美日韩一区二区精品| 日韩欧美一区视频在线观看| 一级毛片女人18水好多| cao死你这个sao货| 悠悠久久av| 国产片内射在线| 五月开心婷婷网| 亚洲三区欧美一区| 操出白浆在线播放| 在线观看一区二区三区激情| 免费一级毛片在线播放高清视频 | 91在线观看av| 夜夜爽天天搞| 村上凉子中文字幕在线| 精品国产乱子伦一区二区三区| videosex国产| 天堂俺去俺来也www色官网| 性少妇av在线| 免费高清视频大片| 国产亚洲精品一区二区www| 99在线视频只有这里精品首页| 欧美av亚洲av综合av国产av| 久久 成人 亚洲| 国产成人精品久久二区二区91| 黑人巨大精品欧美一区二区蜜桃| 亚洲成国产人片在线观看|