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

    葉片包角對(duì)離心泵流場(chǎng)及脈動(dòng)特性的影響

    2017-11-13 08:27:58,,2,,,
    關(guān)鍵詞:包角尾跡蝸殼

    , ,2, , ,

    (1.上海理工大學(xué) 能源與動(dòng)力工程學(xué)院,上海 200093; 2.上海理工大學(xué) 上海市動(dòng)力工程多相流動(dòng)與傳熱重點(diǎn)實(shí)驗(yàn)室,上海 200093; 3.上海船舶設(shè)備研究所,上海 200031)

    葉片包角對(duì)離心泵流場(chǎng)及脈動(dòng)特性的影響

    謝志賓1,陳二云1,2,丁京濤1,楊愛(ài)玲1,李國(guó)平3

    (1.上海理工大學(xué) 能源與動(dòng)力工程學(xué)院,上海200093;2.上海理工大學(xué) 上海市動(dòng)力工程多相流動(dòng)與傳熱重點(diǎn)實(shí)驗(yàn)室,上海200093;3.上海船舶設(shè)備研究所,上海200031)

    針對(duì)離心泵非定常流動(dòng)壓力脈動(dòng)特性,采用滑移網(wǎng)格的大渦模擬技術(shù)對(duì)葉片包角分別為95°,100°,105°,108°的4副葉輪進(jìn)行數(shù)值模擬.分析了葉片包角對(duì)離心泵水力性能、葉輪出口“射流-尾跡”、測(cè)點(diǎn)壓力脈動(dòng)頻譜特性和葉輪徑向力的影響關(guān)系.結(jié)果表明:隨著包角的增大,離心泵的水力性能下降;包角適當(dāng)增大,會(huì)使葉輪射流-尾跡流動(dòng)結(jié)構(gòu)變?nèi)?在設(shè)計(jì)工況下,蝸舌附近測(cè)點(diǎn)壓力脈動(dòng)最大;在蝸殼螺旋段壓力脈動(dòng)強(qiáng)度沿流動(dòng)方向逐漸變?nèi)?而在葉輪流道內(nèi)壓力脈動(dòng)沿流動(dòng)方向逐漸增強(qiáng),在葉輪出口處達(dá)到最大;而離心泵葉輪所受徑向力隨著包角的增大而減小,適當(dāng)?shù)卦龃蟀强梢蕴岣唠x心泵運(yùn)行的可靠性.

    葉片包角; 射流-尾跡; 壓力脈動(dòng); 徑向力

    離心泵是一種將機(jī)械能轉(zhuǎn)換成流體動(dòng)能和壓能的通用機(jī)械,在國(guó)民生產(chǎn)中占有重要地位[1].離心泵在運(yùn)行過(guò)程中,由于非對(duì)稱的蝸殼結(jié)構(gòu)、葉輪的高速旋轉(zhuǎn)及其內(nèi)部流動(dòng)的非定常激勵(lì)力影響,使得離心泵內(nèi)部產(chǎn)生周期性的壓力脈動(dòng),這將引起振動(dòng)和噪聲,降低泵運(yùn)行的可靠性.因此,離心泵內(nèi)部的壓力脈動(dòng)特性引起了許多學(xué)者的高度關(guān)注[2-4].

    葉片包角作為離心泵的重要設(shè)計(jì)參數(shù)之一,對(duì)離心泵的性能和壓力脈動(dòng)特性有十分重要的影響[5].過(guò)小的包角降低了葉片對(duì)流體的控制能力,不利于提高泵的效率;過(guò)大的包角增大了葉片摩擦損失,也不利于性能的提高.目前已有很多學(xué)者針對(duì)葉片包角對(duì)離心泵內(nèi)部流動(dòng)的影響進(jìn)行了廣泛的研究,楊華等[6]采用試驗(yàn)和數(shù)值模擬的方法對(duì)不同包角下的離心泵性能進(jìn)行了試驗(yàn)和模擬對(duì)比,并分析了性能的差異.Yang等[7]采用試驗(yàn)和數(shù)值模擬的方法研究了葉片包角對(duì)離心泵的性能影響,通過(guò)對(duì)比發(fā)現(xiàn),在不同范圍內(nèi),包角大小的變化對(duì)離心泵的性能會(huì)有不同的影響:在相對(duì)小角度下,隨著包角的增加,泵性能提升;而在相對(duì)大角度下,卻出現(xiàn)相反的結(jié)論.張翔等[8]基于CFD數(shù)值計(jì)算方法分析了不同葉片包角與離心泵性能的關(guān)系,得出了離心泵揚(yáng)程和功率特性發(fā)生變化的原因在于葉輪出口相對(duì)速度液流角發(fā)生改變.王勇等[9]采用試驗(yàn)的方法分析了不同葉片包角對(duì)離心泵空化振動(dòng)流動(dòng)誘導(dǎo)噪聲的影響.

    葉片包角問(wèn)題引起了許多學(xué)者的高度關(guān)注,但是,大多數(shù)學(xué)者的研究只關(guān)注葉片包角的變化對(duì)離心泵整機(jī)性能的影響,而對(duì)葉片包角變化影響葉輪出口“射流-尾跡”及壓力脈動(dòng)的研究相對(duì)較少.作者通過(guò)數(shù)值求解方法,分析不同葉片包角對(duì)離心泵內(nèi)部流動(dòng)特征及壓力脈動(dòng)及葉輪徑向力的影響,為離心泵葉輪葉片包角的設(shè)計(jì)提供一定的參考依據(jù).

    1 計(jì)算模型和數(shù)值模擬方法

    1.1計(jì)算模型及網(wǎng)格

    研究的原型泵型號(hào)為IS100-80-160單級(jí)單吸離心泵,泵的設(shè)計(jì)流量QS=100 m3/h,設(shè)計(jì)揚(yáng)程H=32 m,轉(zhuǎn)速n=2 900 r/min,模型的基本參數(shù)如表1所示.整個(gè)離心泵的計(jì)算域如圖1所示.為了提高計(jì)算精度和收斂性,蝸殼和葉輪流體域全部用六面體網(wǎng)格生成,并且在葉片和蝸殼的壁面均進(jìn)行邊界層的加密,蝸殼壁面第一層網(wǎng)格尺度為0.012 mm,保證壓力脈動(dòng)分析的準(zhǔn)確性[10].

    表1 離心泵主要幾何參數(shù)Tab.1 Main geometric parameters of the centrifugal pump

    圖1 離心泵計(jì)算域模型Fig.1 Computational domain of the centrifugal pump

    1.2包角的選定

    葉片包角是離心泵葉片的一個(gè)重要幾何參數(shù),其值的大小實(shí)質(zhì)表征了葉片流道的擴(kuò)散程度[11].葉片包角φ定義為葉片入口邊與圓心的連線和出口邊與圓心連線間的夾角[12],如圖2所示(見(jiàn)下頁(yè)).在對(duì)許多具有較好水利性能的葉輪進(jìn)行研究的基礎(chǔ)上,有學(xué)者給出了Zφ/360隨著比轉(zhuǎn)速ns變化的最優(yōu)規(guī)律,如表2所示[13].原始離心泵模型的比轉(zhuǎn)速為131.12,由此根據(jù)表2確定了模型泵葉片包角的取值范圍為90~108 °,保證離心泵的主要基本參數(shù)不變,選用葉片包角依次為95,100,105,108 °這4副葉輪,通過(guò)數(shù)值求解的方法來(lái)研究包角變化對(duì)離心泵內(nèi)部非定常流動(dòng)壓力脈動(dòng)特性和葉輪徑向力的影響,從而為離心泵的設(shè)計(jì)與優(yōu)化提供參考.

    1.3大渦模擬控制方程及邊界條件

    圖2 葉片包角示意圖Fig.2 Schematic diagram of the blade wrag angle

    表2 葉片數(shù)Z與包角φ的乘積與比轉(zhuǎn)速ns的關(guān)系

    離心泵內(nèi)部的流動(dòng)是十分復(fù)雜的三維紊流流動(dòng),而大渦模擬能夠捕捉到離心泵葉輪流道內(nèi)復(fù)雜的湍流流動(dòng)現(xiàn)象[14].大渦模擬是計(jì)算流體力學(xué)中的一種十分重要的數(shù)值模擬方法,它是基于一種數(shù)學(xué)濾波函數(shù),通過(guò)在大渦流場(chǎng)中引入應(yīng)力項(xiàng)來(lái)過(guò)濾掉小渦對(duì)大渦的影響,在流體機(jī)械非定常流動(dòng)模擬計(jì)算中具有很大優(yōu)勢(shì)[15].本文以模型泵為基礎(chǔ),通過(guò)求解三維不可壓非定常流動(dòng)N-S方程獲得離心泵內(nèi)非定常流動(dòng)的數(shù)值解,初始流場(chǎng)由離心泵流場(chǎng)的定常數(shù)值解近似給出,其控制方程為

    (1)

    (2)

    式中:帶上劃線的量是濾波后的場(chǎng)變量;ρ為流體的密度;μ為流體的動(dòng)力黏性系數(shù);τij為亞格子尺度應(yīng)力;t為時(shí)間;ui,uj為i,j方向上的雷諾時(shí)均速度.

    2 計(jì)算結(jié)果分析

    2.1外特性曲線分析

    從流量-揚(yáng)程功率曲線圖3中可以看出,隨著葉片包角的增大,離心泵的揚(yáng)程H在相應(yīng)流量下明顯下降。Q/QS為流量比.軸功率P隨著葉片包角的增大而減小,包角增大,相應(yīng)流量下消耗的功率變小,葉片包角從95 °增加到108 °,設(shè)計(jì)工況下?lián)P程下降了3.83%,軸功率下降了2.92%.

    圖3 流量-揚(yáng)程曲線Fig.3 Lift against flow coefficient

    從流量-效率曲線圖4可以發(fā)現(xiàn),在小流量情況下,隨著葉片包角的增大,泵的效率η有升高的趨勢(shì),而在大流量下,離心泵的效率因包角的增大而顯著下降.4副葉輪的最高效率在設(shè)計(jì)工況下分別為81.95%,81.84%,81.77%,81.19%,效率下降了0.76%.包角為95 °時(shí),效率最高,這是因?yàn)槿~片包角增大后,葉片變長(zhǎng),當(dāng)量擴(kuò)散角變小,葉片對(duì)流體的約束能力增強(qiáng)的同時(shí)摩擦損失也越來(lái)越大,從理論上來(lái)講,必定存在最優(yōu)的葉片包角.同時(shí),隨著葉片包角的增大,泵最高效率點(diǎn)的流量也相對(duì)減小,根據(jù)葉輪和蝸殼的耦合關(guān)系[16],葉片包角的增大使得葉片出口安放角和葉輪出口相對(duì)液流角減小,從而會(huì)引起泵最高效率點(diǎn)流量偏向小流量.由4副葉輪的性能計(jì)算結(jié)果可見(jiàn),隨著包角的增大,離心泵的揚(yáng)程下降且效率減小.95 °葉片包角的性能要明顯優(yōu)于其他幾副葉輪.在葉輪外尺寸確定后,造成這些性能差異的根本原因是包角變化引起葉片出口角、流道長(zhǎng)度和葉型的改變,而不同的葉型使得葉輪內(nèi)部的流場(chǎng)結(jié)構(gòu)不同,包角的改變導(dǎo)致葉輪內(nèi)部流動(dòng)損失差異和葉片出口液體的滑移,最終使得葉片包角改變了離心泵的外特性.

    2.2葉輪出口流動(dòng)狀態(tài)

    為了研究包角的改變對(duì)葉輪出口流動(dòng)狀態(tài)的影響,圖5給出了葉輪出口徑向速度Cr沿周向的分布.其中:X=0,表示葉片壓力面位置;X=1,表示葉片吸力面位置.葉輪出口20%葉高為靠盤側(cè),50%葉高為跨盤蓋流道中央,80%葉高處為靠蓋側(cè).U2為周向速度.可以看出,在葉輪出口處徑向速度沿周向分布非常的不均勻,葉片的壓力面?zhèn)却嬖谥鴱较蛩俣鹊母咚賲^(qū),而在葉片的吸力面?zhèn)葎t對(duì)應(yīng)著徑向速度的低速區(qū).這實(shí)際上反映的就是葉輪出口的“射流-尾跡”流動(dòng)結(jié)構(gòu);由于尾跡區(qū)實(shí)質(zhì)上是低速流體微團(tuán)的聚集區(qū),其存在將導(dǎo)致流體在葉輪流道的出口處產(chǎn)生“堵塞”,減小有效通流面積,增加速度分布的不均勻性,影響葉輪出口流動(dòng)狀態(tài).從圖5中可以看出,靠盤側(cè)葉輪出口徑向速度隨包角變化不明顯,而從跨盤蓋中央與靠蓋側(cè)整體上來(lái)看,在X=0.5處,徑向速度“虧損”最大,徑向速度隨著包角的增大先呈現(xiàn)一定的平緩趨勢(shì)而后波動(dòng)變大.這說(shuō)明過(guò)大或者過(guò)小的包角加強(qiáng)了葉輪出口射流-尾跡流動(dòng)結(jié)構(gòu),使葉輪內(nèi)部由邊界層分離、漩渦二次流和尾跡等引起的渦流損失變大,而在一定范圍內(nèi)相對(duì)大葉片包角減弱了葉輪出口的射流-尾跡,使葉輪內(nèi)部的渦流損失變小.

    圖4 流量-效率曲線Fig.4 Efficiency against flow coefficient

    圖5 不同包角下葉輪出口徑向速度Cr沿周向分布Fig.5 Radial velocity along circumferential direction at different wrap angles

    2.3壓力脈動(dòng)分析

    由于流體的黏性作用及葉輪出口與蝸殼進(jìn)口動(dòng)靜部件的相互干涉作用,使得離心泵內(nèi)的流場(chǎng)呈現(xiàn)非定常流動(dòng)特征,從而引起流場(chǎng)的壓力脈動(dòng).這種周期性的壓力脈動(dòng)會(huì)引起系統(tǒng)及設(shè)備的振動(dòng)和噪聲,嚴(yán)重時(shí)會(huì)導(dǎo)致設(shè)備損壞.

    為了監(jiān)測(cè)離心泵在設(shè)計(jì)工況下內(nèi)部壓力脈動(dòng),引入了無(wú)量綱壓力系數(shù)Cp=Δp/(0.5ρU2),其中,Δp為壓力與平均壓力的差值,ρ為水的密度.在離心泵蝸殼壁面的中截面上及沿葉輪流道中間位置設(shè)置了監(jiān)測(cè)點(diǎn)y1,y2,y3,y4,y5和p1,p2,p3,p4,p5,如圖6所示.

    圖7給出了不同包角下蝸殼螺旋段測(cè)點(diǎn)的壓力脈動(dòng)頻譜特性,f為頻率.從圖7中可以看出,不同包角下的離心泵壓力脈動(dòng)信號(hào)均呈現(xiàn)出明顯的離散特性,并且各測(cè)點(diǎn)壓力脈動(dòng)的主頻均在葉頻處.壓力脈動(dòng)最大的地方出現(xiàn)在蝸舌附近的p1點(diǎn)處,在此位置處的壓力脈動(dòng)低頻信號(hào)也比較豐富.從p1點(diǎn)起,流體沿流動(dòng)方向各測(cè)點(diǎn)在主頻上的壓力脈動(dòng)強(qiáng)度逐漸減弱,高頻成分逐漸消失.這是因?yàn)槭艿饺~輪出口處的射流-尾跡作用及蝸殼斷面幾何形狀影響,隔舌處徑向間隙小,葉輪與蝸殼的動(dòng)靜干涉作用最強(qiáng)烈,沿著流動(dòng)方向徑間隙變大,動(dòng)量摻混作用也就隨之變?nèi)?

    圖6 監(jiān)測(cè)點(diǎn)位置示意圖Fig.6 Schematic diagram of monitoring points

    圖7 不同包角蝸殼螺旋段測(cè)點(diǎn)的壓力脈動(dòng)頻譜圖Fig.7 Frequency spectrum of the pressure pulsation in the spiral case with different wrap angles

    蝸舌附近是離心泵內(nèi)部距離葉輪最近的區(qū)域,也是葉輪通道氣流流出直接沖擊的部位,該部位最能體現(xiàn)壓力的非定常變化,選擇蝸舌處監(jiān)測(cè)點(diǎn)p1進(jìn)行研究.圖8給出了不同包角下蝸舌處測(cè)點(diǎn)壓力脈動(dòng)頻譜特性,從圖8中可以看出,在p1處監(jiān)測(cè)點(diǎn)壓力脈動(dòng)的低頻特性十分復(fù)雜,并且不同葉片包角在葉頻處壓力脈動(dòng)最為強(qiáng)烈,由于蝸舌的存在和葉輪出口氣流的擾動(dòng)作用使得蝸舌附近的流場(chǎng)較為復(fù)雜,表明了葉頻脈動(dòng)是壓力脈動(dòng)的主要貢獻(xiàn)因素.而除了葉頻外,在2倍葉頻和3倍葉頻處也存在脈動(dòng)相對(duì)明顯的峰值,隨后幅值逐漸減小,這是因?yàn)樵谌~輪旋轉(zhuǎn)過(guò)程中監(jiān)測(cè)點(diǎn)受前一個(gè)葉片產(chǎn)生的壓力脈動(dòng)的影響還未消失時(shí),后一個(gè)葉片又從該位置經(jīng)過(guò).

    圖8 不同包角蝸舌處測(cè)點(diǎn)壓力脈動(dòng)頻譜圖Fig.8 Frequency spectrum of the pressure pulsation at the measuring point on the volute tongue

    同時(shí),為了更加清楚地呈現(xiàn)不同包角對(duì)壓力脈動(dòng)強(qiáng)度的影響,對(duì)于p1點(diǎn)在葉頻下的壓力脈動(dòng)強(qiáng)度用直方圖表示,如圖9所示.可以清晰地看出,隨著包角的增大,壓力脈動(dòng)幅值先增大后減小,95 °包角時(shí)壓力脈動(dòng)最弱,105 °包角時(shí)壓力脈動(dòng)最強(qiáng),葉片包角從95 °增加到105 °時(shí),脈動(dòng)幅值增加了14.7%.該現(xiàn)象與前面分析葉輪出口處射流-尾跡結(jié)構(gòu)現(xiàn)象不一致,這可能是由于此時(shí)葉輪與蝸舌處的動(dòng)靜干涉作用強(qiáng)于射流-尾跡結(jié)構(gòu)帶來(lái)的壓力脈動(dòng)強(qiáng)度.

    圖10 (見(jiàn)下頁(yè))給出了設(shè)計(jì)工況下不同包角葉輪流道內(nèi)測(cè)點(diǎn)的壓力脈動(dòng)頻譜特性,可以發(fā)現(xiàn),在葉輪流道內(nèi)沿著流體流動(dòng)方向,各測(cè)點(diǎn)的壓力脈動(dòng)幅值在葉頻及諧頻處均增大,位于葉輪進(jìn)口處的測(cè)點(diǎn)y1的壓力脈動(dòng)最弱,而距離葉輪出口處最近測(cè)點(diǎn)y5的壓力脈動(dòng)則最為強(qiáng)烈.由此說(shuō)明葉輪內(nèi)壓力脈動(dòng)沿著流體流動(dòng)的方向逐漸增強(qiáng),在葉輪出口處壓力脈動(dòng)表現(xiàn)得最為強(qiáng)烈,這是因?yàn)槿~輪流道的出口處射流-尾跡結(jié)構(gòu)作用明顯且距離蝸殼壁更近,動(dòng)靜耦合干擾作用造成出口處湍流十分強(qiáng)烈,從而導(dǎo)致脈動(dòng)幅值很大.同時(shí)可以發(fā)現(xiàn),隨著包角的增大,葉輪出口y5壓力脈動(dòng)整體呈現(xiàn)遞減,這與包角增大導(dǎo)致葉輪出口處射流-尾跡減弱的規(guī)律一致.

    圖9 不同包角蝸舌測(cè)點(diǎn)壓力脈動(dòng)強(qiáng)度直方圖Fig.9 Frequency histogram of the pressure pulsation at the measuring point on the volute tongue

    2.4葉輪徑向力分析

    從圖11中可以看出,離心泵葉輪所受徑向力呈現(xiàn)明顯的周期性和脈動(dòng)性,波峰和波谷的數(shù)目與葉片數(shù)呈現(xiàn)對(duì)應(yīng)關(guān)系,并且隨著葉片包角的增大,徑向力幅值呈現(xiàn)減小趨勢(shì).圖12 (見(jiàn)下頁(yè))為葉輪所受徑向力的頻譜分析,徑向力呈現(xiàn)出明顯的離散特性,不同包角下葉輪徑向力最大幅值均出現(xiàn)在葉頻處,同時(shí)在二階和三階葉頻及諧頻處都出現(xiàn)了較大諧波.為了能夠更加清楚地呈現(xiàn)不同包角下徑向力的變化情況,圖13 (見(jiàn)下頁(yè))給出了葉頻下不同包角的徑向力大小直方圖,可以清楚地看到,葉頻下徑向力的幅值隨著包角的增大而逐漸減小.在葉頻下葉片包角從95 °增加到108 °時(shí),徑向力下降了16.0%,由此可見(jiàn),包角增大,增強(qiáng)了離心泵運(yùn)行的可靠性.壓力沿葉輪出口分布不均導(dǎo)致離心泵葉輪受到徑向作用力,包角的改變對(duì)葉輪徑向力的規(guī)律同對(duì)葉輪出口處射流-尾跡的規(guī)律一致,證明了之前包角對(duì)葉輪出口流場(chǎng)影響結(jié)論的準(zhǔn)確性,也說(shuō)明包角改變了葉輪出口的壓力分布,從而影響了葉輪徑向力的產(chǎn)生.

    圖10 不同包角葉輪流道內(nèi)測(cè)點(diǎn)壓力脈動(dòng)頻譜圖Fig.10 Frequency spectrum of the pressure pulsation at the measuring point in the impeller with different wrap angles

    圖11 不同包角葉輪徑向力差時(shí)域圖Fig.11 Time-domain diagram of the radial force difference of the impeller with different wrap angles

    圖12 不同包角葉輪徑向力差大小頻域圖Fig.12 Frequency spectrum of the radial force difference of the impeller with different wrap angles

    圖13 不同包角徑向力差大小直方圖Fig.13 Frequency histogram of the radial force difference of the impeller with different wrap angles

    3 結(jié) 論

    a. 隨著葉片包角的增大,葉輪出口相對(duì)液流角減小,進(jìn)而改變了離心泵的外特性,具體表現(xiàn)為,隨著葉片包角的增加,離心泵的揚(yáng)程下降,離心泵的效率降低,且最大值向小流量方向偏移.

    b. 包角的變化改變了葉輪流道的結(jié)構(gòu).從而使得葉輪內(nèi)部的流場(chǎng)發(fā)生變化,產(chǎn)生了射流-尾跡結(jié)構(gòu),隨著包角的增大,葉輪內(nèi)部射-流尾跡作用先減弱后加強(qiáng).

    c. 蝸殼螺旋段測(cè)點(diǎn)沿流動(dòng)方向,尾跡流的動(dòng)量摻混作用隨距離尾緣越遠(yuǎn)而越弱;沿葉輪流道流動(dòng)方向,壓力脈動(dòng)逐漸增強(qiáng),葉輪出口處壓力脈動(dòng)最為強(qiáng)烈;隨著包角的增大,蝸舌處壓力脈動(dòng)強(qiáng)度先增大后減小,而葉輪出口處的壓力脈動(dòng)卻先減小后增大,這是射流尾跡結(jié)構(gòu)作用和葉輪與蝸殼動(dòng)量摻混作用相互耦合導(dǎo)致的.

    d. 葉輪所受徑向力表現(xiàn)出明顯的周期性,且徑向力脈動(dòng)頻率在葉頻處強(qiáng)度最大;徑向力隨著葉片包角的增大而減小,因此,適當(dāng)?shù)卦龃笕~片包角,可以提高離心泵運(yùn)行的可靠性.

    [1] TAN L,CAO S L,GUI S B.Hydraulic design and pre-whirl regulation law of inlet guide vane for centrifugal pump[J].Science China Technological Sciences,2010,53(8):2142-2151.

    [2] KAUPERT K A,STAUBLI T.The unsteady pressure field in a high specific speed centrifugal pump impeller-part I:influence of the volute[J].Journal of Fluids Engineering,1999,121(3):621-626.

    [3] GUELICH J F,BOLLETER U.Pressure pulsations in centrifugal pumps[J].Journal of Vibration and Acoustics,1992,114(2):272-279.

    [4] DONG R,CHU S,KATZ J.Effect of modification to tongue and impeller geometry on unsteady flow,pressure fluctuations,and noise in a centrifugal pump[J].Journal of Turbomachinery,1997,119(3):506-515.

    [5] 張啟華,李月仙,林建忠,等.基于紙漿泵流場(chǎng)數(shù)值模擬的葉輪設(shè)計(jì)準(zhǔn)則及其應(yīng)用[J].江蘇大學(xué)學(xué)報(bào)(自然科學(xué)版),2011,32(2):185-189.

    [6] 楊華,劉超,湯方平,等.不同葉片包角的離心泵試驗(yàn)與數(shù)值模擬[J].機(jī)械工程學(xué)報(bào),2007,43(10):166-169.

    [7] YANG S S,KONG F Y,CHEN H,et al.Effects of blade wrap angle influencing a pump as turbine[J].Journal of Fluids Engineering,2012,134(6):061102.

    [8] 張翔,王洋,徐小敏,等.葉片包角對(duì)離心泵性能的影響[J].農(nóng)業(yè)機(jī)械學(xué)報(bào),2010,44(11):38-42.

    [9] 王勇,劉厚林,劉東喜,等.葉片包角對(duì)離心泵流動(dòng)誘導(dǎo)振動(dòng)噪聲的影響[J].農(nóng)業(yè)工程學(xué)報(bào),2013,29(1):72-77.

    [10] 馬尊領(lǐng),郭艷麗,陳二云,等.離心泵壓力脈動(dòng)的數(shù)值模擬及敏感性分析[J].工程熱物理學(xué)報(bào),2015,36(7):1461-1465.

    [11] 嚴(yán)敬,張江源,何敏,等.離心泵可控包角圓柱形葉片型線方程[J].排灌機(jī)械,2008,26(5):46-49.

    [12] 離心泵設(shè)計(jì)基礎(chǔ)編寫組.離心泵設(shè)計(jì)基礎(chǔ)[M].北京:機(jī)械工業(yè)出版社,1974:134-162.

    [13] 陳乃祥,吳玉林.離心泵[M].北京:機(jī)械工業(yè)出版社,2003:105-117.

    [14] BYSKOV R K,JACOBSEN C B,PEDERSEN N.Flow in a centrifugal pump impeller at design and off-design conditions-Part II:large eddy simulations[J].Journal of Fluids Engineering,2003,125(1):73-83.

    [15] 王福軍.計(jì)算流體動(dòng)力學(xué)分析:CFD軟件原理與應(yīng)用[M].北京:清華大學(xué)出版社,2004:15-35.

    [16] 張翔.低比速離心泵蝸殼對(duì)水力性能的影響[D].鎮(zhèn)江:江蘇大學(xué),2008.

    EffectofBladeWrapAngleontheInternalFlowFieldandPressurePulsationinCentrifugalPumps

    XIE Zhibin1,CHENEryun1,2,DINGJingtao1,YANGAiling1,LIGuoping3

    (1.SchoolofEnergyandPowerEngineering,UniversityofShanghaiforScienceandTechnology,Shanghai200093,China; 2.ShanghaiKeyLaboratroyofMultiphaseFlowandHeatTransferinPowerEngineer,UniversityofShanghaiforScienceandTechnology,Shanghai200093,China; 3.ShanghaiMarineEquipmentResearchInstitute,Shanghai200031,China)

    For investigating the pressure pulsation characteristics of the unsteady flow in centrifugal pumps,a single-stage single-suction centrifugal pump was taken as an example to perform a numerical simulation by using the technology of large eddy simulation with sliding grid.The three dimensional flow field was simulated numerically with the blade wrap angle 95 °,100 °,105° and 108° respectively.The influences of the blade wrap angle on the centrifugal pump hydraulic performance,impeller exit jet-wake,pressure pulsation frequency spectrum characteristics at each measuring point and radial force of the impeller were analyzed.The results show that as the angle increases,the centrifugal pump hydraulic performance degradates.The angle increase will lead to the impeller outlet jet-wake flow structure.Under design conditions,near the volute tongue the pressure pulsation at measuring points is the strongest,at the spiral screw section along the flow direction,the pressure pulsation intensity gradually weaken,and at the impeller passage point along the flow direction,the pressure pulsation gradually enhances and reaches the maximum at the impeller outlet.The radial force of the centrifugal pump impeller decreases with the increase of angle.Therefore,appropriately increasing the wrap angle can improve the reliability of the centrifugal pump running.

    bladewrapangle;jet-wake;pressurepulsation;radialforce

    1007-6735(2017)05-0430-08

    10.13255/j.cnki.jusst.2017.05.005

    2017-03-14

    國(guó)家自然科學(xué)基金資助項(xiàng)目(51106099,50976072);上海市科委科研計(jì)劃項(xiàng)目(13DZ2260900)

    謝志賓(1990-),男,碩士研究生.研究方向:葉輪機(jī)械氣動(dòng)力學(xué).E-mail:xiezhibin2015@163.com

    陳二云(1977-),男,副教授.研究方向:葉輪機(jī)械氣動(dòng)力學(xué)和聲學(xué).E-mail:cheneryun@usst.edu.cn

    TH311

    A

    (編輯:石 瑛)

    猜你喜歡
    包角尾跡蝸殼
    塔式太陽(yáng)能發(fā)電系統(tǒng)S-CO2向心透平氣動(dòng)性能優(yōu)化
    一種基于Radon 變換和尾跡模型的尾跡檢測(cè)算法
    焊接蝸殼泵模型開(kāi)發(fā)
    葉片包角對(duì)水輪機(jī)模式多級(jí)液力透平性能的影響
    汽車后包角拉延成形數(shù)值模擬及驗(yàn)證
    污水泵蝸殼斷裂原因分析
    某柴油發(fā)動(dòng)機(jī)輪系優(yōu)化設(shè)計(jì)分析
    汽車科技(2019年6期)2019-12-31 06:13:41
    基于EEMD-Hilbert譜的渦街流量計(jì)尾跡振蕩特性
    蝸殼差壓法在印尼ASAHAN一級(jí)水電站的應(yīng)用
    高功率低燃油耗的新一代蝸殼式增壓器
    精品欧美国产一区二区三| av.在线天堂| 精品久久久久久久人妻蜜臀av| 免费看日本二区| 91久久精品国产一区二区成人| 综合色av麻豆| 亚洲最大成人av| 成年女人永久免费观看视频| 国产亚洲91精品色在线| 精品午夜福利在线看| 日本wwww免费看| 七月丁香在线播放| 国产欧美另类精品又又久久亚洲欧美| 亚洲国产欧美在线一区| 久久精品夜夜夜夜夜久久蜜豆| 亚洲国产欧美人成| 97热精品久久久久久| 国产一级毛片在线| 精品久久久久久成人av| 午夜久久久久精精品| 天堂影院成人在线观看| 亚洲在线观看片| 高清在线视频一区二区三区 | 免费观看精品视频网站| 十八禁国产超污无遮挡网站| 久久久久久久亚洲中文字幕| 亚洲在线观看片| 看片在线看免费视频| 国产亚洲精品av在线| 少妇熟女aⅴ在线视频| av天堂中文字幕网| 国产高潮美女av| 日本午夜av视频| 国产精品一区二区三区四区久久| 久久这里只有精品中国| 亚洲精品成人久久久久久| 九九在线视频观看精品| 色网站视频免费| 看十八女毛片水多多多| 两性午夜刺激爽爽歪歪视频在线观看| 国产成人freesex在线| www.av在线官网国产| 亚洲经典国产精华液单| 欧美区成人在线视频| av.在线天堂| 国产爱豆传媒在线观看| 国产乱人视频| 全区人妻精品视频| 国产精品永久免费网站| 看十八女毛片水多多多| 成人国产麻豆网| 欧美一区二区亚洲| 国产亚洲一区二区精品| 久久人妻av系列| 精品久久久久久久久久久久久| 成人欧美大片| 色播亚洲综合网| 午夜精品一区二区三区免费看| 非洲黑人性xxxx精品又粗又长| 午夜激情欧美在线| 久久99热6这里只有精品| 国产69精品久久久久777片| 如何舔出高潮| 国产高清三级在线| 赤兔流量卡办理| 特级一级黄色大片| 亚洲经典国产精华液单| 国产精品一区二区性色av| av国产免费在线观看| 少妇的逼好多水| 男人的好看免费观看在线视频| 内地一区二区视频在线| 国产精品久久久久久久电影| 在线免费观看的www视频| 六月丁香七月| av线在线观看网站| 老司机影院毛片| 国产精品久久久久久精品电影| 欧美zozozo另类| 欧美成人精品欧美一级黄| 特级一级黄色大片| 精品久久久久久电影网 | 美女被艹到高潮喷水动态| 亚洲中文字幕一区二区三区有码在线看| 亚洲一级一片aⅴ在线观看| 亚洲va在线va天堂va国产| 亚洲天堂国产精品一区在线| 1000部很黄的大片| 国产亚洲最大av| 午夜激情欧美在线| 亚洲欧美日韩东京热| 自拍偷自拍亚洲精品老妇| 乱人视频在线观看| 久久精品久久久久久久性| 亚洲,欧美,日韩| 哪个播放器可以免费观看大片| 一夜夜www| 女的被弄到高潮叫床怎么办| 在线观看66精品国产| 亚洲伊人久久精品综合 | 亚洲在久久综合| 婷婷色综合大香蕉| 97热精品久久久久久| 观看免费一级毛片| 麻豆久久精品国产亚洲av| 国产亚洲91精品色在线| 少妇人妻精品综合一区二区| 国产亚洲精品av在线| 国产午夜精品久久久久久一区二区三区| 欧美性猛交黑人性爽| 51国产日韩欧美| 国产激情偷乱视频一区二区| 国产亚洲5aaaaa淫片| 一级毛片久久久久久久久女| 国产91av在线免费观看| 小说图片视频综合网站| 国产黄a三级三级三级人| 国产亚洲精品久久久com| 国产色爽女视频免费观看| 欧美区成人在线视频| 亚洲一级一片aⅴ在线观看| 欧美色视频一区免费| 国产精品电影一区二区三区| 亚洲精品,欧美精品| 欧美一区二区亚洲| 久久久久久久久久久丰满| 三级经典国产精品| 国产精品,欧美在线| 午夜免费激情av| 成人美女网站在线观看视频| 国产女主播在线喷水免费视频网站 | 免费观看a级毛片全部| 最近中文字幕高清免费大全6| 极品教师在线视频| 精品不卡国产一区二区三区| 亚洲内射少妇av| 51国产日韩欧美| 国产伦理片在线播放av一区| 久久精品久久久久久噜噜老黄 | 国产免费一级a男人的天堂| 少妇的逼好多水| 三级国产精品片| 国内精品一区二区在线观看| 久久人人爽人人爽人人片va| 成人二区视频| 午夜福利在线观看免费完整高清在| 亚洲国产精品专区欧美| 亚洲最大成人av| 尤物成人国产欧美一区二区三区| 国产人妻一区二区三区在| 日本一二三区视频观看| 嫩草影院新地址| 美女内射精品一级片tv| 少妇人妻精品综合一区二区| 在线播放国产精品三级| 伦精品一区二区三区| 亚洲精品aⅴ在线观看| 婷婷六月久久综合丁香| 欧美性猛交黑人性爽| 18禁裸乳无遮挡免费网站照片| 如何舔出高潮| 国产又色又爽无遮挡免| 永久网站在线| 我要看日韩黄色一级片| 免费黄网站久久成人精品| 亚洲中文字幕日韩| 欧美性感艳星| 中文在线观看免费www的网站| 日本-黄色视频高清免费观看| 国产亚洲午夜精品一区二区久久 | 久久久久久伊人网av| 国产在视频线精品| 联通29元200g的流量卡| 欧美区成人在线视频| 欧美激情国产日韩精品一区| 久久久久久久午夜电影| 超碰av人人做人人爽久久| 99久国产av精品国产电影| 久久久久久久久中文| 午夜a级毛片| 最近手机中文字幕大全| 日韩在线高清观看一区二区三区| 亚洲综合色惰| 国产69精品久久久久777片| 精品国产一区二区三区久久久樱花 | 成年女人永久免费观看视频| 精品99又大又爽又粗少妇毛片| 亚洲国产高清在线一区二区三| 我要看日韩黄色一级片| 高清在线视频一区二区三区 | 在现免费观看毛片| 中文资源天堂在线| 中文精品一卡2卡3卡4更新| 国产视频内射| 日本-黄色视频高清免费观看| 99九九线精品视频在线观看视频| 大又大粗又爽又黄少妇毛片口| 亚洲自拍偷在线| 国产精华一区二区三区| 国产精品.久久久| 国产黄a三级三级三级人| 亚洲欧美精品综合久久99| 国产成人精品久久久久久| 亚洲av一区综合| 中文在线观看免费www的网站| 久久99精品国语久久久| 禁无遮挡网站| 免费看光身美女| 高清日韩中文字幕在线| 国内精品宾馆在线| 麻豆乱淫一区二区| h日本视频在线播放| 男女边吃奶边做爰视频| 亚洲人成网站在线播| 日韩人妻高清精品专区| 免费观看的影片在线观看| 欧美+日韩+精品| 亚洲av福利一区| 毛片女人毛片| 亚洲最大成人中文| 亚洲熟妇中文字幕五十中出| 久久久国产成人精品二区| 伦精品一区二区三区| 五月玫瑰六月丁香| 晚上一个人看的免费电影| 色综合色国产| 久久久精品94久久精品| 桃色一区二区三区在线观看| 精品国产一区二区三区久久久樱花 | 我的女老师完整版在线观看| 在现免费观看毛片| 午夜视频国产福利| 亚洲av免费在线观看| 国语对白做爰xxxⅹ性视频网站| 日韩欧美在线乱码| 又粗又硬又长又爽又黄的视频| 成人二区视频| 午夜福利在线在线| 欧美色视频一区免费| 国产精品.久久久| av免费在线看不卡| 啦啦啦啦在线视频资源| 人妻夜夜爽99麻豆av| 变态另类丝袜制服| 我要搜黄色片| 又爽又黄a免费视频| 亚洲电影在线观看av| 国产黄色小视频在线观看| 如何舔出高潮| 久久99精品国语久久久| 精品久久久久久久末码| 国产伦理片在线播放av一区| 久久久成人免费电影| a级毛色黄片| 高清视频免费观看一区二区 | 热99在线观看视频| 人妻少妇偷人精品九色| 在线播放国产精品三级| 国产伦精品一区二区三区四那| 全区人妻精品视频| 精品人妻视频免费看| 99热这里只有是精品50| 久久精品国产亚洲av涩爱| 国产极品精品免费视频能看的| 婷婷色麻豆天堂久久 | 成年av动漫网址| 午夜视频国产福利| 亚洲欧美中文字幕日韩二区| 国产精品一二三区在线看| 国产在视频线在精品| 国产伦一二天堂av在线观看| 你懂的网址亚洲精品在线观看 | 高清毛片免费看| 精品99又大又爽又粗少妇毛片| 国产精品熟女久久久久浪| 韩国高清视频一区二区三区| 看片在线看免费视频| av在线蜜桃| 1000部很黄的大片| 男人舔奶头视频| 天堂中文最新版在线下载 | 淫秽高清视频在线观看| 免费看美女性在线毛片视频| 国产免费视频播放在线视频 | 偷拍熟女少妇极品色| 97热精品久久久久久| 桃色一区二区三区在线观看| 伦精品一区二区三区| 日本免费在线观看一区| 精品久久久久久久久av| 国产色爽女视频免费观看| 又粗又爽又猛毛片免费看| 男女视频在线观看网站免费| 三级国产精品片| 又黄又爽又刺激的免费视频.| 国产淫语在线视频| 久久久久久久久久成人| 高清午夜精品一区二区三区| 黄片无遮挡物在线观看| 天堂√8在线中文| 国产精品熟女久久久久浪| 国产91av在线免费观看| 免费av观看视频| 久久久久久久久久久免费av| 亚洲国产高清在线一区二区三| 亚洲av中文av极速乱| 久久久成人免费电影| av在线天堂中文字幕| 日韩高清综合在线| 日韩一区二区视频免费看| 在线a可以看的网站| 国产不卡一卡二| 日韩高清综合在线| 午夜免费男女啪啪视频观看| 亚洲国产精品sss在线观看| 日韩欧美精品免费久久| 免费人成在线观看视频色| 国产成人免费观看mmmm| 免费在线观看成人毛片| 国产探花极品一区二区| 亚洲精品乱久久久久久| 亚洲一区高清亚洲精品| 熟女电影av网| or卡值多少钱| 欧美另类亚洲清纯唯美| 亚洲三级黄色毛片| 精品不卡国产一区二区三区| 国产不卡一卡二| 亚洲四区av| 精品久久久久久久人妻蜜臀av| 日本色播在线视频| 69av精品久久久久久| 春色校园在线视频观看| 亚洲精华国产精华液的使用体验| 亚洲美女视频黄频| 国产麻豆成人av免费视频| 最近中文字幕2019免费版| 国产成年人精品一区二区| 亚洲国产色片| 国产免费男女视频| 国产高清视频在线观看网站| 成人美女网站在线观看视频| 91aial.com中文字幕在线观看| 久久6这里有精品| 亚洲人成网站高清观看| 成人毛片60女人毛片免费| 一本一本综合久久| av在线亚洲专区| 亚洲国产精品专区欧美| 欧美3d第一页| 国产午夜精品久久久久久一区二区三区| 两个人视频免费观看高清| 久久人人爽人人片av| 亚洲精品自拍成人| 中文乱码字字幕精品一区二区三区 | 97超视频在线观看视频| 国产av一区在线观看免费| 亚洲精品亚洲一区二区| 少妇熟女欧美另类| 我的老师免费观看完整版| 国产亚洲精品久久久com| 亚洲在久久综合| 久久久成人免费电影| 久久久久久久久久久丰满| 国产爱豆传媒在线观看| 韩国av在线不卡| 亚洲电影在线观看av| 国产精华一区二区三区| 精品酒店卫生间| 日韩av在线大香蕉| 天堂av国产一区二区熟女人妻| 男女边吃奶边做爰视频| 午夜久久久久精精品| 色吧在线观看| 蜜桃亚洲精品一区二区三区| 久久精品夜夜夜夜夜久久蜜豆| 国模一区二区三区四区视频| 永久网站在线| 国产精品久久久久久精品电影小说 | 白带黄色成豆腐渣| 岛国在线免费视频观看| 日韩一区二区视频免费看| 亚洲内射少妇av| 99久久成人亚洲精品观看| av在线蜜桃| 午夜爱爱视频在线播放| 色网站视频免费| 亚洲精品国产av成人精品| 国产伦在线观看视频一区| 精品国产三级普通话版| 久久人妻av系列| 自拍偷自拍亚洲精品老妇| 干丝袜人妻中文字幕| 寂寞人妻少妇视频99o| 久久久久九九精品影院| 大话2 男鬼变身卡| 丰满少妇做爰视频| 国产欧美日韩精品一区二区| 国产精品,欧美在线| 一卡2卡三卡四卡精品乱码亚洲| 日韩av不卡免费在线播放| 深爱激情五月婷婷| 老师上课跳d突然被开到最大视频| 精品酒店卫生间| 久久久久精品久久久久真实原创| 国产亚洲精品久久久com| av在线老鸭窝| 精品熟女少妇av免费看| 国产精品爽爽va在线观看网站| 精品一区二区三区视频在线| 国产成人a区在线观看| 日产精品乱码卡一卡2卡三| 欧美一区二区亚洲| 中文字幕精品亚洲无线码一区| 久久精品国产亚洲av涩爱| 少妇被粗大猛烈的视频| 嫩草影院精品99| 亚洲欧美一区二区三区国产| 国产日韩欧美在线精品| 精品熟女少妇av免费看| 国产高清视频在线观看网站| 午夜精品在线福利| 国产综合懂色| 你懂的网址亚洲精品在线观看 | 一级毛片我不卡| 乱码一卡2卡4卡精品| 蜜臀久久99精品久久宅男| 日本av手机在线免费观看| 欧美97在线视频| 欧美一区二区精品小视频在线| 长腿黑丝高跟| 欧美3d第一页| 3wmmmm亚洲av在线观看| 最近2019中文字幕mv第一页| 日本与韩国留学比较| 国产精品美女特级片免费视频播放器| 国产午夜福利久久久久久| 免费观看人在逋| 久久久久久久久久久免费av| 欧美性感艳星| 国产精品一区二区在线观看99 | 亚洲成人中文字幕在线播放| 直男gayav资源| 久久精品人妻少妇| 国产一级毛片在线| 精品国产一区二区三区久久久樱花 | 99热这里只有精品一区| 网址你懂的国产日韩在线| 爱豆传媒免费全集在线观看| 日本爱情动作片www.在线观看| 精品无人区乱码1区二区| 久久亚洲国产成人精品v| 久久99精品国语久久久| 亚洲av一区综合| 国产在线一区二区三区精 | 国产午夜精品论理片| 男人舔女人下体高潮全视频| 两个人视频免费观看高清| 内地一区二区视频在线| 熟女电影av网| 成人漫画全彩无遮挡| 男人狂女人下面高潮的视频| 免费无遮挡裸体视频| 男人舔女人下体高潮全视频| 久久6这里有精品| 天堂影院成人在线观看| 亚洲第一区二区三区不卡| av在线播放精品| 中文字幕久久专区| 久久99热这里只有精品18| 中文字幕久久专区| 精品午夜福利在线看| kizo精华| 国产亚洲最大av| 欧美激情久久久久久爽电影| 最近视频中文字幕2019在线8| 伦精品一区二区三区| or卡值多少钱| 亚洲欧美日韩东京热| 深爱激情五月婷婷| 久久亚洲国产成人精品v| 国产免费视频播放在线视频 | 日韩成人av中文字幕在线观看| 久久精品熟女亚洲av麻豆精品 | 国产老妇女一区| 中文乱码字字幕精品一区二区三区 | 国产激情偷乱视频一区二区| 国产精品久久久久久精品电影| 久久久国产成人精品二区| 99热网站在线观看| 国产精品久久电影中文字幕| 国产亚洲91精品色在线| 欧美一级a爱片免费观看看| 亚洲av电影在线观看一区二区三区 | 日韩人妻高清精品专区| 成人特级av手机在线观看| 老女人水多毛片| 噜噜噜噜噜久久久久久91| 欧美潮喷喷水| av视频在线观看入口| 国产探花极品一区二区| 男插女下体视频免费在线播放| 国语自产精品视频在线第100页| 青春草亚洲视频在线观看| 免费av不卡在线播放| 国产精品久久久久久av不卡| 日韩欧美 国产精品| 能在线免费观看的黄片| 欧美一区二区精品小视频在线| 久热久热在线精品观看| 亚洲精品成人久久久久久| 最近的中文字幕免费完整| 久久热精品热| 国产视频首页在线观看| 人人妻人人澡人人爽人人夜夜 | 噜噜噜噜噜久久久久久91| 秋霞伦理黄片| 国产精品一区二区三区四区久久| 不卡视频在线观看欧美| 狠狠狠狠99中文字幕| 身体一侧抽搐| 日本欧美国产在线视频| 中文欧美无线码| 久热久热在线精品观看| 又粗又硬又长又爽又黄的视频| 亚洲在线自拍视频| 精品免费久久久久久久清纯| 身体一侧抽搐| 国产美女午夜福利| 人人妻人人澡人人爽人人夜夜 | 国产精品熟女久久久久浪| 国产探花在线观看一区二区| 精品不卡国产一区二区三区| 亚洲激情五月婷婷啪啪| 免费看美女性在线毛片视频| 女的被弄到高潮叫床怎么办| 亚洲欧美日韩无卡精品| 简卡轻食公司| 边亲边吃奶的免费视频| 精品免费久久久久久久清纯| 欧美潮喷喷水| 97在线视频观看| 欧美另类亚洲清纯唯美| 国产爱豆传媒在线观看| 美女内射精品一级片tv| 国产精品人妻久久久久久| 国产精品久久久久久精品电影| 亚洲最大成人手机在线| 久久精品熟女亚洲av麻豆精品 | 日韩欧美精品免费久久| 国产精品99久久久久久久久| 国产男人的电影天堂91| 久久久久久久久久久丰满| 狂野欧美激情性xxxx在线观看| 精品久久久噜噜| 亚州av有码| 日韩,欧美,国产一区二区三区 | 亚洲第一区二区三区不卡| 免费看日本二区| 免费看a级黄色片| 亚洲人成网站高清观看| 免费不卡的大黄色大毛片视频在线观看 | 最近2019中文字幕mv第一页| 在线观看一区二区三区| 丝袜喷水一区| 国产高清有码在线观看视频| 国产午夜精品论理片| 美女cb高潮喷水在线观看| 男的添女的下面高潮视频| 别揉我奶头 嗯啊视频| 国产亚洲最大av| 一个人看视频在线观看www免费| 国产片特级美女逼逼视频| 男女啪啪激烈高潮av片| 精品久久久久久久久久久久久| av免费在线看不卡| 亚洲性久久影院| 欧美性猛交黑人性爽| 中文亚洲av片在线观看爽| 久久精品久久久久久久性| av免费观看日本| 麻豆一二三区av精品| 偷拍熟女少妇极品色| 成人毛片60女人毛片免费| 少妇人妻精品综合一区二区| 日本色播在线视频| 欧美另类亚洲清纯唯美| 嫩草影院精品99| 亚洲国产日韩欧美精品在线观看| 丰满少妇做爰视频| 身体一侧抽搐| 男人和女人高潮做爰伦理| 少妇被粗大猛烈的视频| 国产av一区在线观看免费| 日本免费a在线| 精品少妇黑人巨大在线播放 | 亚洲国产精品久久男人天堂| 七月丁香在线播放| 色尼玛亚洲综合影院| 观看美女的网站| 少妇熟女aⅴ在线视频| 热99re8久久精品国产| 日韩欧美精品v在线| 国产激情偷乱视频一区二区| 99九九线精品视频在线观看视频| 我的老师免费观看完整版| 秋霞在线观看毛片| 亚洲精品乱码久久久久久按摩| 一卡2卡三卡四卡精品乱码亚洲| 成人国产麻豆网| 国产成人一区二区在线| 免费搜索国产男女视频| 亚洲欧美精品专区久久| 内地一区二区视频在线| 免费观看人在逋|