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

    葉片鋸齒前緣控制流動分離的數(shù)值研究

    2021-08-08 16:36:10孫貴洋陳二云楊愛玲
    能源研究與信息 2021年4期

    孫貴洋 陳二云 楊愛玲

    摘要:針對鋸齒前緣結構調控葉片近壁面流場特性,以 NACA0018葉片為對象,采用大渦模擬方法研究不同鋸齒前緣結構對葉片近壁面流場的影響機制。獲得了來流速度為30 m· s?1、雷諾數(shù)為513440、0°攻角下葉片近壁面流場分布特性。分析了鋸齒前緣和葉片前緣和尾緣處壓力脈動及分離渦的影響。數(shù)值結果表明:對正弦波齒而言,隨著振幅的增大,在波谷處的小渦開始向前緣移動,整體上小尺度渦增多,前緣近壁面壓力脈動增大,尾緣近壁面壓力脈動減小;對疊加波形齒而言,尾跡渦進一步破碎,厚度變薄,葉片表面出現(xiàn)破碎的小尺度渦,在尾緣處葉片壓力脈動幅值下降最為明顯,且未出現(xiàn)明顯的窄帶尖峰。

    關鍵詞:尾跡渦;壓力脈動;鋸齒前緣;疊加波形

    中圖分類號: TB17??? 文獻標志碼: A

    Numerical study on the flow separation controlling by sawtooth leading edge of the blade

    SUN Guiyang ,CHEN Eryun,YANG Ailing

    (School of Energy and Power Engineering/Shanghai Key Laboratory of Multiphase Flow and Heat Transfer inPower Engineering, University of Shanghai for Science and Technology, Shanghai 200093, China)

    Abstract:Aiming at the flow field characteristics of the near-wall blade surface controlled by sawtooth leading edge, taking the NACA0018 blade as the object, the influence of sawtooth leadingedge on the near-wall flow field of the blade was studied by large eddy simulation method. Thedistribution characteristics of flow field near the blade wall at a flow velocity of 30 m · s? ,? a Reynolds number of 513440, and an angle of attack of 0° were obtained. The influence of the sawtooth leading edge on the pressure fluctuation and the separation vortex at the leading and trailing edges of the blade was analyzed. The numerical results show that for the sine-shaped sawtooth, as the amplitude increases, the small-scale vortex at the wave trough moves to the frontedge, which increases generally. The pressure fluctuation near the front wall increases while thepressure fluctuation near the wall of the trailing edge decreases. For the superimposed sawtooth, the wake vortex is further broken, whose thickness decreases. And the small-scale broken vortices appear on the blade surface. The pressure fluctuation amplitude at the blade trailing edge decreases obviously, and there is no obvious narrow-band spike.

    Keywords:wake vortex; pressure fluctuation; sawtooth leading edge; superimposed wave

    隨著人民生活的日益提高,噪聲問題慢慢進入人們的視野,風機、飛機、風力機及空調壓縮機的噪聲給人們的日常生活帶來了諸多不便。葉片氣動噪聲是這類機械最主要的噪聲源。國內外學者在研究葉片自噪聲時,對葉片尾緣自噪聲關注度頗高,也采取了許多降噪措施,例如:尾緣鋸齒、多孔介質和刷式尾緣等,并取得了顯著成果[1]。Moreau 等[2]以平板為對象,研究雷諾數(shù)在1.6×105~4.2×105范圍內鋸齒對平板噪聲的影響,發(fā)現(xiàn)尾緣鋸齒可以削弱低頻范圍內的寬頻噪聲約3 dB,同時削弱尾渦脫落噪聲約13 dB。 Chong 等[3]探究了層流葉片窄帶單峰噪聲產(chǎn)生機理與尾緣鋸齒對葉片窄帶單峰噪聲的影響,并指出原始葉片的聲學反饋回路,在尾緣處發(fā)生了衍射,出現(xiàn)了 T?S 波,從而解釋了寬帶噪聲中“駝峰”形成的原因。仝帆等[4]以 SD2030葉片為對象,經(jīng)研究發(fā)現(xiàn)鋸齒尾緣可以明顯降低葉片中低頻范圍內的噪聲,在4000 Hz 以下,窄帶噪聲最多可降低約16 dB。與葉片鋸齒尾緣相對的是葉片前緣鋸齒,王國付[5]通過粒子圖像測速(PIV )流場測試給出了鋸齒前緣葉片和光滑葉片的流場變化形態(tài),得出凹谷和凸包的存在使得流場產(chǎn)生流向渦對,通過流向渦對的作用引起凹谷從前緣的層流邊界層轉涙為湍流邊界層再附,同時凸包處層流邊界層附體流動。馬群毅等[6]將鋸齒前緣應用于透平靜葉柵中,研究其在不同攻角下的流動特性,結果表明在大負攻角工況下,結節(jié)葉片可將均勻來流沿展向分割成幾部分,每部分氣流由兩個旋向相反的小渦組成,并沿流向逐漸融合,避免了常規(guī)葉片在壓力面?zhèn)妊卣麄€展向大渦的形成與發(fā)展,降低了流動損失,提高了流動效率。陳偉杰等[7]研究了雷諾數(shù)在2×105~8×105范圍內前緣鋸齒結構對葉片層流邊界層不穩(wěn)定噪聲的影響,得出前緣鋸齒可以減弱甚至完全抑制邊界層不穩(wěn)定噪聲,降噪效果對鋸齒振幅和鋸齒波長均比較敏感,鋸齒振幅越大,波長越小,降噪效果越好,降噪量可達30 dB,且前緣鋸齒結構可以誘導產(chǎn)生流向渦,影響葉片下游邊界層流動,破壞聲學反饋回路。

    前緣鋸齒靈感來源于仿座頭鯨鰭肢前緣的不規(guī)則波狀凸起。目前,大部分學者研究的前緣鋸齒為正弦波或三角鋸齒結構。本文主要在正弦波鋸齒前緣結構基礎上,研究疊加波形鋸齒調控葉片近壁面的壓力脈動及漩渦流場分布特性,為后期研究疊加波形降噪設計提供理論支持。

    1數(shù)值計算方法

    1.1計算模型

    本文中計算模型為 NACA0018葉片,弦長 C =250 mm。為保持原始葉片幾何特征不變,僅在距前緣30%內的部分設計前緣鋸齒結構,前緣鋸齒設計公式為

    y = ai sin (bix)???????????? (1)

    式中:ai = hi ,ai 和 hi 均為振幅;bi =2πC/λi,bi 為波形的角頻率;λi為波長;x 為自變量,表示橫坐標; y 為函數(shù)值,表示縱坐標。

    正弦波葉片的相對波長λ/C =25%,相對振幅變化范圍為hi/C =5.0%~7.5%。每種疊加波形葉片前緣鋸齒形狀均由兩個正弦波疊加而成,正弦波振幅均為h/C =2.5%,其中每組疊加波形的正弦波都有一個相同的波長,λ/C =25%。計算模型為原始 NACA0018葉片、兩組正弦波前緣鋸齒葉片和兩組疊加波前緣鋸齒葉片,共計5組。鋸齒前緣結構葉片及模型幾何參數(shù)分別如圖1、表1所示,其中圖1中單波形葉片上1、2和3分別代表波峰、基線和波谷,疊加波形葉片上1、2、3和4分別代表波峰1、波谷2、波峰2和波谷1。

    1.2網(wǎng)格劃分

    本文中網(wǎng)格拓撲結構采用“O”型結構。圖2為數(shù)值計算區(qū)域結構尺寸示意圖,采用分區(qū)方式劃分網(wǎng)格。圖3為鋸齒前緣葉片網(wǎng)格。經(jīng)過無量綱化的壁面距離y+計算,近壁面第一層網(wǎng)格高度為0.01 mm,滿足?y+≤1的要求,增長率為1.05。A 域繞葉片一圈網(wǎng)格節(jié)點數(shù)為484,葉片展向節(jié)點數(shù)為72,B域 X方向節(jié)點數(shù)為141,Y 方向節(jié)點數(shù)為101,Z方向節(jié)點數(shù)為30,A 域網(wǎng)格數(shù)為6.5×106,B 域網(wǎng)格數(shù)為1.5×106。由于鋸齒前緣葉片的前緣特殊結構,使得波峰與波谷的曲率最大,且最易發(fā)生分離流動,故此處網(wǎng)格應加密處理,且為保證計算精度與計算時間的平衡,故網(wǎng)格尺度從波峰與波谷處逐漸向兩側增長來劃分網(wǎng)格。鋸齒前緣葉片網(wǎng)格的劃分方式如圖3(a)所示。

    1.3流場邊界條件

    本文中計算流場時使用Fluent商業(yè)軟件,借助 k?ω SST 湍流模型進行定常計算。以定常的收斂值為初始值,運用大渦模擬動力Smagorin? sky?Lilly 亞格子模型進行瞬態(tài)計算。邊界條件為:入口為速度入口,進口來流速度U =30 m · s?1,葉片攻角α=0?;流場出口為自由出流;葉片展向方向為平移周期性邊界條件,左右與速度入口邊界條件一致,葉片表面為無滑移壁面條件。時間步長為1.0×10?5 s ,流場計算約2500個時間步長后收斂,流場收斂后繼續(xù)計算2500個時間步長獲得流場統(tǒng)計平均結果。

    1.4網(wǎng)格無關性驗證

    為驗證仿真的準確性,建立與文獻[8]中實驗模型一致的仿真模型進行網(wǎng)格驗證,選擇k ?! SST兩方程湍流模型進行網(wǎng)格驗證。網(wǎng)格無關性參數(shù)比較如表2所示。葉片不同網(wǎng)格數(shù)時表面靜壓分布如圖4所示,其中Cp為靜壓系數(shù)。從圖中可以看出,計算結果與實驗值吻合很好,靜壓分布大體一致,但在吸力面20%~30%內,仿真值高于實驗值,而在壓力面60%~100%內,仿真值低于實驗值。從表2中可以看出,隨著網(wǎng)格數(shù)增加,分離點的位置先增大后減小并逐漸趨于穩(wěn)定,升阻比一直增大,并逐漸趨于穩(wěn)定。故可以確定網(wǎng)格數(shù)在500萬~700萬范圍內時可忽略網(wǎng)格數(shù)對數(shù)值計算結果的影響。

    靜壓系數(shù)計算式為

    p ? p0

    式中: p為葉片壁面靜壓, Pa ;p0為入口靜壓, Pa;ρ為空氣密度, kg ·m?3。

    1.5湍流模型驗證

    本文中選擇S? A 、RNG k ?"及 k?ω SST 三種湍流模型進行驗證。不同湍流模型葉片表面靜壓分布如圖5所示。由圖中可以看出:采用 k ?! SST湍流模型計算的仿真值與實驗值在前緣處吻合較好,而且兩者整體靜壓基本吻合;采用 k ?" RNG 和S? A湍流模型計算的仿真值與實驗值在前緣處偏差較大。因此,本文中選擇k?! SST兩方程湍流模型作為定常計算模型。該湍流模型能較好地捕捉更多的流動分離現(xiàn)象,并且對葉片前緣靜壓捕捉與實驗值接近,可以用于模擬葉片繞流流動。

    2計算結果分析

    2.1鋸齒前緣結構對葉片尾跡的影響

    根據(jù)文獻[4, 9–10]引入Q準則渦核等值面瞬時速度云圖,分析葉片尾跡渦核變化規(guī)律。圖6為 Q 準則渦核等值面圖瞬時速度云圖。由圖中可知,鋸齒前緣結構使葉片的卡門渦街尾跡發(fā)生明顯的變化,不再呈現(xiàn)規(guī)則的周期性渦交錯現(xiàn)象,尾跡渦厚度變薄,尾跡渦發(fā)生擠壓、拉伸,從而使大尺度渦破碎成小尺度渦。緊接著,發(fā)現(xiàn)葉片表面渦結構也發(fā)生了顯著的改變。與 NACA0018葉片相比,鋸齒前緣結構葉片在接近尾緣處,葉片表面開始出現(xiàn)破碎的小尺度渦,渦系結構更加緊致且沿流向不斷擴大。

    隨著h/C的不斷增大,正弦波形葉片在波谷區(qū)的小渦開始向前緣移動,整體上小尺度渦增多,而疊加波形鋸齒前緣葉片與正弦波鋸齒前緣葉片相比,尾跡渦厚度更薄,渦系結構相對寬松且小尺度渦區(qū)域更廣。在低馬赫數(shù)下,可以認為流動噪聲來源于渦。根據(jù) Powell 渦聲方程可知,渦量的減少有利于降低流動噪聲,大渦破碎生成小渦的過程,伴隨著能量的衰減,但小渦在高頻段對噪聲有較高的貢獻量[11]。

    為了進一步剖析葉片展向不同位置尾跡渦結構,分別在模型 B1的波谷、基線及波峰處做二維切面,分析二維切面處鋸齒前緣結構對尾跡渦的影響。

    圖7為二維切面渦量分布。從圖中可以發(fā)現(xiàn),模型 B1的波峰、基線及波谷三個位置處的切面均未出現(xiàn)與模型 A 葉片尾跡一樣的周期性交錯渦現(xiàn)象,并且尾跡渦量相對較小。鋸齒前緣結構葉片從波峰到波谷有效弦長逐漸減小,渦逐漸向前緣移動,有效弦長的減小使得氣流附體流動的時間減小,流動更加劇烈,葉片表面渦量從波峰到波谷逐漸增大,但整體小于模型 A 葉片表面渦量。葉片在波谷處的氣流受到剪切、擠壓、卷吸作用,在主流挾帶下向尾緣移動,此切面氣流受到的剪切力最大且流動最為復雜,故可解釋渦量相對較大的原因。鋸齒前緣結構不僅僅改變了葉片尾跡流動結構,從葉片前緣至尾緣的渦量分布可以看出流動結構也發(fā)生了巨大改變。通過控制全局影響尾跡的變化,減小了尾跡及葉片表面渦量大小,達到了流動控制的目的。

    2.2鋸齒前緣結構對葉片前緣和尾緣處壓力脈動的影響

    為了探究葉片前緣和尾緣處壓力脈動變化,分別在距葉片前緣和尾緣法向近壁面設置測點監(jiān)測此處的壓力脈動,監(jiān)測點位置分布如圖8所示。從整體上看,改型后的葉片前緣和尾緣處壓力脈動幅值明顯下降,尖峰均出現(xiàn)在1500 Hz 以內,且尾緣處壓力脈動幅值相對較大。

    圖9為葉片前緣和尾緣近壁面測點壓力脈動。從圖中可以看出, NACA0018葉片前緣和尾緣處壓力脈動幅值明顯高于鋸齒前緣結構葉片,其中模型 C1與模型 C2尾緣處壓力脈動未出現(xiàn)明顯峰值,且壓力脈動最大下降值為50 Pa,這表明此組合方式形成的疊加波形鋸齒前緣結構葉片降低壓力脈動幅值表現(xiàn)較好。壓力脈動的大小不僅與脈動幅值有關,還與壓力脈動相位角有關。假如兩個相位角相反的脈動源,盡管壓力脈動幅值較大,也會由于相位角的差異,最終疊加后的壓力脈動幅值可能不大。原始葉片展向一致,則壓力脈動相位角在展向上相差無幾,葉片展向結構越不規(guī)則,壓力脈動相位角出現(xiàn)不同的可能性越大,所以正弦波鋸齒前緣葉片與疊加波形葉片壓力脈動幅值低于原始葉片壓力脈動幅值。

    從圖9中可知,壓力脈動主峰頻率在改型后發(fā)生頻移現(xiàn)象,而壓力脈動主峰頻率與分離渦脫落頻率有關,壓力脈動主峰頻率發(fā)生了改變,則說明鋸齒前緣結構改變了原始葉片的分離渦脫落的頻率,從而能進一步影響渦的發(fā)展。根據(jù) Lockard[12]和Casalino等[13]的研究,選擇葉片表面作為積分面,關注偶極子聲源而忽略單極子與四極子聲源,可知壓力脈動幅值的降低有利于降低偶極子源噪聲。鋸齒前緣結構降低了葉片前緣和尾緣處壓力脈動幅值,這對削弱偶極子源噪聲有利。

    3結論

    (1)鋸齒前緣結構使葉片的卡門渦街尾跡發(fā)生了明顯的變化,不再呈現(xiàn)規(guī)則的周期性渦交錯現(xiàn)象,尾跡渦厚度變薄。與 NACA0018葉片相比,鋸齒前緣葉片在接近尾緣處,葉片表面開始出現(xiàn)破碎的小尺度渦,渦系結構更加緊致且沿流向不斷擴大。隨著h/C的不斷增大,葉片表面小渦區(qū)域開始向前緣移動,整體上小尺度渦增多。其中疊加波形鋸齒前緣葉片與正弦波鋸齒前緣葉片相比,尾跡渦厚度更薄,渦系結構相對寬松且小尺度渦區(qū)域更大。

    (2)疊加波形鋸齒前緣葉片壓力脈動幅值下降最為明顯,在尾緣處基本未出現(xiàn)明顯的窄帶峰值。鋸齒前緣葉片展向不同位置處,波谷切面所在尾緣位置壓力波動最大,其與原始葉片相比壓力脈動主峰頻率發(fā)生頻移現(xiàn)象,壓力脈動幅值最大下降50 Pa。

    (3)鋸齒前緣結構改變了原始葉片壓力脈動主峰頻率,壓力脈動主峰頻率與渦脫落的頻率有關,從而可知鋸齒前緣結構也改變了原始葉片渦脫落的頻率。

    參考文獻:

    [1]謝志賓, 陳二云, 丁京濤, 等.葉片包角對離心泵流場及脈動特性的影響[J].上海理工大學學報 , 2017, 39(5):430-437.

    [2] MOREAU? D? J,? DOOLAN? C? J. Tonal? noise? fromtrailing? edge? serrations? at? low? Reynolds? number[R]. AIAA 2013?2010, Reston: AIAA, 2013.

    [3] CHONG? T? P,? JOSEPH? P,? GRUBER? M. Anexperimental? study? of? airfoil? instability? noise? with trailing? edge? serrations[C]//16th? AIAA/CEASAeroacoustics Conference. Stockholm, Sweden: AIAA, 2010.

    [4] 仝帆, 喬渭陽, 王良鋒, 等.仿生學翼型尾緣鋸齒降噪機理[J].航空學報, 2015, 36(9):2911-2922.

    [5] 王國付, 張明明, 徐建中.仿鯨魚鰭翼段氣動控制實驗研究[J].工程熱物理學報, 2013, 34(10):1842-1846.

    [6] 馬群毅, 鄧清華, 豐鎮(zhèn)平.基于座頭鯨胸鰭結節(jié)的仿生透平葉片內流特性研究[C]//高等學校工程熱物理第二十屆全國學術會議論文集——熱流體力學專輯.青島, 中國:中國高等教育學會工程熱物理專業(yè)委員會,2014.

    [7] 陳偉杰, 喬渭陽, 仝帆, 等.前緣鋸齒對邊界層不穩(wěn)定噪聲的影響[J].航空學報, 2016, 37(12):3634-3645.

    [8] JIANG? M,? LI? X? D,? BAI? B? H,? et? al. Numericalsimulation? on? the? NACA0018 airfoil? self-noise

    [9] 楊景茹, 楊愛玲, 陳二云, 等.鋸齒尾緣葉片氣動特性和繞流流場的數(shù)值研究[J].航空動力學報 , 2017, 32(4):900-908.

    [10]戴麗萍, 焦建東, 李新凱, 等.渦發(fā)生器安裝位置對風力機翼型流場影響的數(shù)值研究[J].太陽能學報, 2016,37(2):276-281.

    [11]孫曉峰, 周盛.氣動聲學[M].北京:國防工業(yè)出版社,1994:150?156.

    [12] LOCKARD D P. A comparison of ffowcs Williamshawkings? solvers for airframe noise applications[C]// Proceedings? of? the 8th? AIAA/CEAS? Aeroacoustics Conference & Exhibit. Breckenridge, Colorado: AIAA,2002:2580.

    [13] CASALINO D, JACOB M, ROGER M. Prediction ofrod-airfoil interaction noiseusing the ffowcs-Williams- hawkings analogy[J]. AIAAJournal, 2003, 41(2):182-191.

    成人二区视频| 亚洲国产欧美日韩在线播放| 亚洲国产av影院在线观看| av网站免费在线观看视频| 国产精品.久久久| 国产在线视频一区二区| 日韩电影二区| 久久久久久久亚洲中文字幕| 少妇 在线观看| 日韩欧美精品免费久久| 午夜福利影视在线免费观看| 国精品久久久久久国模美| 一区二区av电影网| 日韩不卡一区二区三区视频在线| 人人妻人人澡人人看| 制服诱惑二区| 国产精品无大码| 久久久午夜欧美精品| 成人毛片60女人毛片免费| 插逼视频在线观看| 夫妻性生交免费视频一级片| 99精国产麻豆久久婷婷| 久久人人爽人人爽人人片va| 麻豆成人av视频| 少妇精品久久久久久久| 国产精品国产三级国产专区5o| videossex国产| 午夜福利视频精品| 美女xxoo啪啪120秒动态图| 成人漫画全彩无遮挡| 自拍欧美九色日韩亚洲蝌蚪91| 这个男人来自地球电影免费观看 | 永久网站在线| 精品一品国产午夜福利视频| 大香蕉97超碰在线| 日本av手机在线免费观看| 99久久中文字幕三级久久日本| 九草在线视频观看| 久久狼人影院| 精品亚洲成a人片在线观看| 九色亚洲精品在线播放| 一本色道久久久久久精品综合| 国产成人精品久久久久久| 高清午夜精品一区二区三区| 亚洲中文av在线| 国产欧美亚洲国产| 国产精品国产av在线观看| 久久狼人影院| av视频免费观看在线观看| 亚洲图色成人| 国产精品久久久久久久久免| 曰老女人黄片| 人人妻人人添人人爽欧美一区卜| 在线观看三级黄色| 成人毛片a级毛片在线播放| 人妻一区二区av| 三上悠亚av全集在线观看| 我要看黄色一级片免费的| 国产男女内射视频| 欧美bdsm另类| a级毛片在线看网站| 久久99一区二区三区| 免费av不卡在线播放| 欧美三级亚洲精品| 九九在线视频观看精品| 大话2 男鬼变身卡| av网站免费在线观看视频| 99视频精品全部免费 在线| 欧美日韩精品成人综合77777| 性色avwww在线观看| 久久青草综合色| 国产伦精品一区二区三区视频9| 亚洲av国产av综合av卡| 制服丝袜香蕉在线| 美女脱内裤让男人舔精品视频| 中国美白少妇内射xxxbb| 精品亚洲成国产av| 特大巨黑吊av在线直播| 麻豆乱淫一区二区| 好男人视频免费观看在线| 美女xxoo啪啪120秒动态图| 国产日韩一区二区三区精品不卡 | 毛片一级片免费看久久久久| 91久久精品电影网| 永久免费av网站大全| 人成视频在线观看免费观看| 成人亚洲欧美一区二区av| 在线精品无人区一区二区三| 伊人久久国产一区二区| 国产亚洲一区二区精品| 夜夜骑夜夜射夜夜干| 免费日韩欧美在线观看| 妹子高潮喷水视频| 高清在线视频一区二区三区| 在线观看一区二区三区激情| 一区二区日韩欧美中文字幕 | 国产视频首页在线观看| 亚洲欧洲日产国产| 18禁裸乳无遮挡动漫免费视频| 久久人人爽人人爽人人片va| 国产视频首页在线观看| 午夜福利视频精品| 国产精品久久久久成人av| 最后的刺客免费高清国语| 日韩亚洲欧美综合| 国产精品熟女久久久久浪| 中文字幕av电影在线播放| 久久久国产一区二区| 国产成人91sexporn| 亚洲天堂av无毛| 制服诱惑二区| 日日撸夜夜添| 日韩精品有码人妻一区| 草草在线视频免费看| 91aial.com中文字幕在线观看| 菩萨蛮人人尽说江南好唐韦庄| 99久久中文字幕三级久久日本| 美女国产视频在线观看| 国产 一区精品| 欧美日韩成人在线一区二区| 日韩伦理黄色片| 亚洲精品日韩av片在线观看| 又大又黄又爽视频免费| 国产成人精品久久久久久| 成人毛片60女人毛片免费| 精品国产一区二区三区久久久樱花| 美女脱内裤让男人舔精品视频| 亚洲精品久久久久久婷婷小说| 精品一区二区免费观看| 国产一区二区三区av在线| 制服诱惑二区| 免费播放大片免费观看视频在线观看| 精品一区二区免费观看| 久久久精品94久久精品| 99九九线精品视频在线观看视频| 啦啦啦视频在线资源免费观看| 男人爽女人下面视频在线观看| 国产精品国产三级国产专区5o| 欧美成人精品欧美一级黄| 美女视频免费永久观看网站| 中文字幕免费在线视频6| 老司机亚洲免费影院| 免费日韩欧美在线观看| 男人爽女人下面视频在线观看| 国产精品嫩草影院av在线观看| 99热6这里只有精品| 亚洲欧美一区二区三区国产| 国产精品久久久久久精品古装| 91成人精品电影| 免费人成在线观看视频色| 日本爱情动作片www.在线观看| 国产精品欧美亚洲77777| 国产成人免费观看mmmm| 啦啦啦在线观看免费高清www| 少妇的逼好多水| 精品一区在线观看国产| 制服人妻中文乱码| 久久午夜综合久久蜜桃| 国产精品免费大片| 色哟哟·www| 夫妻午夜视频| 搡老乐熟女国产| 中国国产av一级| 卡戴珊不雅视频在线播放| xxxhd国产人妻xxx| 国产免费视频播放在线视频| 日韩成人av中文字幕在线观看| 黑人高潮一二区| 在线观看美女被高潮喷水网站| 国产精品熟女久久久久浪| 三上悠亚av全集在线观看| 免费人成在线观看视频色| 欧美激情国产日韩精品一区| 欧美日韩一区二区视频在线观看视频在线| 亚洲欧美一区二区三区黑人 | 国产高清国产精品国产三级| 午夜免费鲁丝| 欧美另类一区| √禁漫天堂资源中文www| 精品卡一卡二卡四卡免费| 一本大道久久a久久精品| 精品一区二区三区视频在线| 下体分泌物呈黄色| 亚洲人与动物交配视频| 国产爽快片一区二区三区| av视频免费观看在线观看| 久久久久久久久久久免费av| 精品熟女少妇av免费看| 亚洲精品日韩av片在线观看| 69精品国产乱码久久久| 成人国语在线视频| 亚洲精品国产色婷婷电影| 一区二区三区免费毛片| 在线 av 中文字幕| 亚洲精品一区蜜桃| 五月开心婷婷网| av播播在线观看一区| 亚洲精品国产色婷婷电影| 下体分泌物呈黄色| av黄色大香蕉| 夫妻性生交免费视频一级片| 久久精品国产自在天天线| 观看美女的网站| 99热全是精品| 在线观看免费视频网站a站| 夜夜看夜夜爽夜夜摸| 亚洲国产av新网站| 卡戴珊不雅视频在线播放| 亚洲精品中文字幕在线视频| 成人毛片60女人毛片免费| 久久久久久久精品精品| 久久久精品免费免费高清| 久久久久精品久久久久真实原创| 少妇人妻精品综合一区二区| 亚洲第一区二区三区不卡| 一级二级三级毛片免费看| 乱人伦中国视频| 日韩人妻高清精品专区| 热re99久久国产66热| 在线观看免费视频网站a站| 26uuu在线亚洲综合色| 国产片特级美女逼逼视频| 久久鲁丝午夜福利片| 日韩av免费高清视频| 18禁观看日本| 亚洲精品视频女| 男男h啪啪无遮挡| 成人午夜精彩视频在线观看| 精品久久蜜臀av无| 999精品在线视频| 高清黄色对白视频在线免费看| 啦啦啦中文免费视频观看日本| 国产免费又黄又爽又色| 美女主播在线视频| 亚洲成色77777| 午夜激情久久久久久久| 亚洲精品乱码久久久v下载方式| 日韩av不卡免费在线播放| 国产精品久久久久久av不卡| 十八禁网站网址无遮挡| av天堂久久9| 高清在线视频一区二区三区| 日本与韩国留学比较| 国产亚洲精品久久久com| 国产成人一区二区在线| 丰满乱子伦码专区| 青春草视频在线免费观看| 久久久久精品久久久久真实原创| 97精品久久久久久久久久精品| 国产高清不卡午夜福利| 美女xxoo啪啪120秒动态图| 九九爱精品视频在线观看| 母亲3免费完整高清在线观看 | 亚洲一区二区三区欧美精品| 丰满少妇做爰视频| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 国产欧美亚洲国产| 亚洲av福利一区| 欧美成人精品欧美一级黄| 久久久久久久久久久久大奶| 亚洲av免费高清在线观看| 亚洲精品美女久久av网站| 亚洲精品色激情综合| 亚洲国产色片| 亚洲精品久久成人aⅴ小说 | 26uuu在线亚洲综合色| 黑人猛操日本美女一级片| 下体分泌物呈黄色| 国产av一区二区精品久久| 国产永久视频网站| 成人18禁高潮啪啪吃奶动态图 | 欧美日本中文国产一区发布| 亚洲欧洲国产日韩| 亚洲,一卡二卡三卡| 九九爱精品视频在线观看| 大陆偷拍与自拍| 日韩中文字幕视频在线看片| 一级毛片我不卡| 国产日韩欧美在线精品| 满18在线观看网站| 日韩av免费高清视频| 国产69精品久久久久777片| 欧美激情极品国产一区二区三区 | 国产精品一区二区在线不卡| 国产亚洲一区二区精品| 国产白丝娇喘喷水9色精品| 国产免费又黄又爽又色| 高清黄色对白视频在线免费看| 国产视频首页在线观看| 欧美3d第一页| 久久精品国产亚洲av涩爱| 亚洲美女视频黄频| 日韩免费高清中文字幕av| 伊人久久精品亚洲午夜| 亚洲精品乱久久久久久| 少妇人妻精品综合一区二区| 午夜福利影视在线免费观看| 满18在线观看网站| 亚洲欧洲日产国产| 久久 成人 亚洲| 亚洲内射少妇av| 午夜91福利影院| 免费黄网站久久成人精品| 久久国产精品大桥未久av| 亚洲av男天堂| 秋霞在线观看毛片| 成人亚洲欧美一区二区av| 免费久久久久久久精品成人欧美视频 | 黄色毛片三级朝国网站| 三上悠亚av全集在线观看| 久久久久久久久久人人人人人人| 啦啦啦在线观看免费高清www| 少妇熟女欧美另类| 少妇高潮的动态图| av免费在线看不卡| 国产成人aa在线观看| kizo精华| 看非洲黑人一级黄片| 亚洲精品国产av成人精品| 另类亚洲欧美激情| 哪个播放器可以免费观看大片| 王馨瑶露胸无遮挡在线观看| 一级毛片aaaaaa免费看小| 国产精品成人在线| 日日摸夜夜添夜夜爱| 久久亚洲国产成人精品v| 国产精品人妻久久久影院| 亚洲图色成人| 亚洲综合色惰| videos熟女内射| 久久精品久久久久久噜噜老黄| 婷婷成人精品国产| 亚洲精品日韩在线中文字幕| 国产乱来视频区| 亚洲一区二区三区欧美精品| 熟女人妻精品中文字幕| 久久免费观看电影| 97精品久久久久久久久久精品| 两个人的视频大全免费| 老司机亚洲免费影院| 肉色欧美久久久久久久蜜桃| 午夜福利影视在线免费观看| 亚洲第一区二区三区不卡| av电影中文网址| 色94色欧美一区二区| 中文字幕av电影在线播放| 伊人久久精品亚洲午夜| 啦啦啦视频在线资源免费观看| 精品少妇黑人巨大在线播放| 啦啦啦视频在线资源免费观看| 午夜福利,免费看| 丰满乱子伦码专区| 91精品一卡2卡3卡4卡| 国产亚洲最大av| av专区在线播放| 精品久久久久久久久av| 欧美变态另类bdsm刘玥| 一级毛片电影观看| 亚洲欧美成人精品一区二区| 视频中文字幕在线观看| 精品人妻熟女毛片av久久网站| 久久久久人妻精品一区果冻| 久久久久久久久久成人| 日日爽夜夜爽网站| 免费高清在线观看日韩| 丰满少妇做爰视频| 天堂中文最新版在线下载| 国产亚洲精品久久久com| 天堂中文最新版在线下载| 男男h啪啪无遮挡| 久久久久久久久久久免费av| 精品久久久精品久久久| 日本黄大片高清| 99久国产av精品国产电影| 国产免费福利视频在线观看| 午夜激情av网站| 日韩一本色道免费dvd| 免费黄色在线免费观看| 91aial.com中文字幕在线观看| 国产无遮挡羞羞视频在线观看| 日本与韩国留学比较| 妹子高潮喷水视频| 精品国产乱码久久久久久小说| 18禁动态无遮挡网站| 26uuu在线亚洲综合色| 狠狠精品人妻久久久久久综合| 国产国语露脸激情在线看| 丝瓜视频免费看黄片| 观看av在线不卡| 国产 精品1| 熟女电影av网| 亚洲中文av在线| 免费人妻精品一区二区三区视频| 国产高清有码在线观看视频| 又黄又爽又刺激的免费视频.| 汤姆久久久久久久影院中文字幕| 男男h啪啪无遮挡| 在线播放无遮挡| 一本—道久久a久久精品蜜桃钙片| 久久99热6这里只有精品| 一级毛片黄色毛片免费观看视频| 亚洲国产精品一区三区| 又黄又爽又刺激的免费视频.| 欧美人与性动交α欧美精品济南到 | 最后的刺客免费高清国语| 伦精品一区二区三区| 大码成人一级视频| 中文字幕久久专区| 久久久久人妻精品一区果冻| 看免费成人av毛片| videossex国产| 免费观看av网站的网址| 精品一区二区三区视频在线| 蜜臀久久99精品久久宅男| 青春草国产在线视频| 欧美精品一区二区免费开放| 国产精品蜜桃在线观看| 女人精品久久久久毛片| 国产色婷婷99| 中文字幕最新亚洲高清| 精品少妇黑人巨大在线播放| 精品卡一卡二卡四卡免费| 夜夜看夜夜爽夜夜摸| 下体分泌物呈黄色| 亚洲美女黄色视频免费看| 啦啦啦啦在线视频资源| 草草在线视频免费看| tube8黄色片| 日本猛色少妇xxxxx猛交久久| 午夜视频国产福利| 18禁在线播放成人免费| 久久久国产精品麻豆| 人妻人人澡人人爽人人| 日韩人妻高清精品专区| 久热久热在线精品观看| av免费在线看不卡| 两个人免费观看高清视频| 亚洲情色 制服丝袜| 国产在线视频一区二区| 国语对白做爰xxxⅹ性视频网站| 国产成人免费无遮挡视频| 欧美激情 高清一区二区三区| av国产精品久久久久影院| 久久女婷五月综合色啪小说| 久久久久网色| 超色免费av| 亚洲av不卡在线观看| 色婷婷av一区二区三区视频| 美女中出高潮动态图| 欧美少妇被猛烈插入视频| 亚洲性久久影院| 中国国产av一级| 亚洲在久久综合| 久久午夜福利片| 熟女电影av网| 免费黄色在线免费观看| 少妇丰满av| 亚洲国产毛片av蜜桃av| 久久午夜综合久久蜜桃| 成人亚洲精品一区在线观看| 国产片特级美女逼逼视频| 2022亚洲国产成人精品| 自线自在国产av| 大码成人一级视频| 色5月婷婷丁香| 男男h啪啪无遮挡| 寂寞人妻少妇视频99o| 秋霞在线观看毛片| 亚洲欧美成人综合另类久久久| 精品亚洲成a人片在线观看| 国产精品免费大片| 国产探花极品一区二区| 婷婷色综合www| 亚洲精品自拍成人| 熟妇人妻不卡中文字幕| 五月玫瑰六月丁香| 国产成人一区二区在线| 能在线免费看毛片的网站| av网站免费在线观看视频| 秋霞在线观看毛片| 草草在线视频免费看| 久久狼人影院| 国产黄片视频在线免费观看| 肉色欧美久久久久久久蜜桃| 三级国产精品片| 国产国语露脸激情在线看| 精品久久久噜噜| 制服人妻中文乱码| 亚洲人成77777在线视频| 久久久久精品性色| 成人黄色视频免费在线看| 新久久久久国产一级毛片| 日韩av免费高清视频| 一个人看视频在线观看www免费| 日日爽夜夜爽网站| 成人毛片60女人毛片免费| 亚洲成人一二三区av| av免费在线看不卡| a级片在线免费高清观看视频| 观看美女的网站| 在线天堂最新版资源| 免费黄色在线免费观看| 国产精品秋霞免费鲁丝片| 成人综合一区亚洲| 91在线精品国自产拍蜜月| 免费观看在线日韩| 性色av一级| 我的女老师完整版在线观看| 国产又色又爽无遮挡免| 亚洲成人手机| 一级,二级,三级黄色视频| 精品99又大又爽又粗少妇毛片| 精品午夜福利在线看| 欧美激情国产日韩精品一区| 国产精品久久久久成人av| 亚洲精品,欧美精品| 新久久久久国产一级毛片| 亚洲欧洲精品一区二区精品久久久 | 亚洲美女黄色视频免费看| 国产日韩欧美亚洲二区| 成人18禁高潮啪啪吃奶动态图 | 久久精品久久精品一区二区三区| 国产成人精品久久久久久| 超碰97精品在线观看| 日韩三级伦理在线观看| 一个人看视频在线观看www免费| 大香蕉97超碰在线| 哪个播放器可以免费观看大片| 国产成人精品一,二区| 人人妻人人添人人爽欧美一区卜| 欧美xxⅹ黑人| 如何舔出高潮| 2018国产大陆天天弄谢| 亚洲国产日韩一区二区| 久久久久精品性色| 亚洲av成人精品一二三区| 80岁老熟妇乱子伦牲交| 少妇猛男粗大的猛烈进出视频| 母亲3免费完整高清在线观看 | 精品久久国产蜜桃| 男人添女人高潮全过程视频| 久久免费观看电影| 狂野欧美激情性bbbbbb| 寂寞人妻少妇视频99o| 在线 av 中文字幕| 天天躁夜夜躁狠狠久久av| 在线 av 中文字幕| av网站免费在线观看视频| 国产精品国产三级国产专区5o| 午夜福利网站1000一区二区三区| 欧美 亚洲 国产 日韩一| 制服丝袜香蕉在线| av国产久精品久网站免费入址| 岛国毛片在线播放| 18禁在线无遮挡免费观看视频| 肉色欧美久久久久久久蜜桃| 国产精品嫩草影院av在线观看| 久久久亚洲精品成人影院| 丰满饥渴人妻一区二区三| 2022亚洲国产成人精品| 久久久久国产精品人妻一区二区| 天堂中文最新版在线下载| 下体分泌物呈黄色| 2022亚洲国产成人精品| 国产日韩欧美视频二区| 一区二区三区精品91| 女的被弄到高潮叫床怎么办| 欧美精品国产亚洲| 日韩熟女老妇一区二区性免费视频| 亚洲人与动物交配视频| 久久久久久久久大av| 亚洲国产精品999| 热99国产精品久久久久久7| 亚洲熟女精品中文字幕| 亚洲三级黄色毛片| 卡戴珊不雅视频在线播放| 亚洲av国产av综合av卡| 人妻系列 视频| 国产精品一国产av| 婷婷成人精品国产| 视频区图区小说| 精品亚洲成a人片在线观看| 各种免费的搞黄视频| 国产免费又黄又爽又色| 国产视频首页在线观看| 国产在线视频一区二区| 一区二区三区免费毛片| 国产成人精品无人区| 日韩一区二区三区影片| 观看av在线不卡| 乱人伦中国视频| 啦啦啦在线观看免费高清www| kizo精华| 男男h啪啪无遮挡| 三级国产精品欧美在线观看| 亚洲国产av新网站| 欧美一级a爱片免费观看看| 狂野欧美激情性xxxx在线观看| 日韩,欧美,国产一区二区三区| 精品人妻熟女av久视频| 久久久久久伊人网av| 99国产精品免费福利视频| 国产老妇伦熟女老妇高清| 九九久久精品国产亚洲av麻豆| 亚洲欧美成人精品一区二区| 免费日韩欧美在线观看| 草草在线视频免费看| 人妻一区二区av| 下体分泌物呈黄色| 女人精品久久久久毛片| 国产欧美亚洲国产| 国产日韩欧美视频二区| 大片电影免费在线观看免费| 久久久国产精品麻豆| 热re99久久精品国产66热6|