• <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在线观看| 高清欧美精品videossex| 男女高潮啪啪啪动态图| 亚洲精品aⅴ在线观看| 性高湖久久久久久久久免费观看| 午夜福利免费观看在线| 男人舔女人的私密视频| 欧美日韩av久久| 国产av国产精品国产| 色吧在线观看| 亚洲一区二区三区欧美精品| 亚洲av电影在线进入| 黄片小视频在线播放| 搡老岳熟女国产| 青春草国产在线视频| 少妇猛男粗大的猛烈进出视频| 成人国语在线视频| 日本色播在线视频| 天天躁夜夜躁狠狠久久av| 伦理电影大哥的女人| 免费看不卡的av| 久久精品国产a三级三级三级| 少妇人妻精品综合一区二区| 免费黄网站久久成人精品| 欧美成人精品欧美一级黄| 97在线人人人人妻| 考比视频在线观看| 看十八女毛片水多多多| 国产精品 国内视频| 国产一区二区三区av在线| 黄片无遮挡物在线观看| 人人妻人人澡人人看| 日韩伦理黄色片| 色吧在线观看| 黄色怎么调成土黄色| 亚洲,一卡二卡三卡| 亚洲国产精品成人久久小说| 狠狠精品人妻久久久久久综合| 国产精品香港三级国产av潘金莲 | 飞空精品影院首页| 亚洲精品国产av蜜桃| 精品一品国产午夜福利视频| 亚洲精品美女久久av网站| 人人妻人人澡人人爽人人夜夜| 亚洲,欧美,日韩| 欧美在线黄色| 人人妻,人人澡人人爽秒播 | 国产男人的电影天堂91| 亚洲一卡2卡3卡4卡5卡精品中文| 人成视频在线观看免费观看| www.av在线官网国产| 久热这里只有精品99| 嫩草影视91久久| 国产精品一区二区精品视频观看| 午夜日韩欧美国产| av天堂久久9| 桃花免费在线播放| 男人操女人黄网站| 成人毛片60女人毛片免费| 国产精品.久久久| 国产精品一区二区在线观看99| 黄片小视频在线播放| 伊人久久国产一区二区| 国产亚洲av高清不卡| 亚洲成人手机| 成人18禁高潮啪啪吃奶动态图| 777米奇影视久久| 国产极品粉嫩免费观看在线| 一级,二级,三级黄色视频| 国产男女超爽视频在线观看| 人妻一区二区av| 一区福利在线观看| 十八禁高潮呻吟视频| 亚洲av综合色区一区| 亚洲欧美中文字幕日韩二区| 精品国产国语对白av| 看免费av毛片| 男女边吃奶边做爰视频| 一本—道久久a久久精品蜜桃钙片| 男人操女人黄网站| 亚洲精品视频女| 欧美日韩一区二区视频在线观看视频在线| 赤兔流量卡办理| 无遮挡黄片免费观看| 纯流量卡能插随身wifi吗| 午夜激情av网站| 日韩免费高清中文字幕av| 在线观看免费高清a一片| 男女高潮啪啪啪动态图| 可以免费在线观看a视频的电影网站 | 国产av一区二区精品久久| 国产野战对白在线观看| 国产爽快片一区二区三区| 韩国av在线不卡| 久久久久网色| 丰满少妇做爰视频| 男女床上黄色一级片免费看| 国产在线一区二区三区精| 日韩中文字幕欧美一区二区 | 黄片播放在线免费| 国产精品 欧美亚洲| 亚洲精品乱久久久久久| 国产欧美日韩综合在线一区二区| 免费高清在线观看视频在线观看| 中文天堂在线官网| 免费女性裸体啪啪无遮挡网站| 国产免费现黄频在线看| 亚洲男人天堂网一区| 80岁老熟妇乱子伦牲交| 国产淫语在线视频| 少妇精品久久久久久久| 亚洲av在线观看美女高潮| 欧美亚洲日本最大视频资源| 精品少妇一区二区三区视频日本电影 | 黑人欧美特级aaaaaa片| 国产精品欧美亚洲77777| xxx大片免费视频| av片东京热男人的天堂| 亚洲欧美成人综合另类久久久| 一区在线观看完整版| 中文欧美无线码| 亚洲人成电影观看| 日韩人妻精品一区2区三区| 精品国产露脸久久av麻豆| 大陆偷拍与自拍| 亚洲中文av在线| 国产欧美亚洲国产| 大香蕉久久网| 国产一区二区三区综合在线观看| 国产又色又爽无遮挡免| 午夜福利乱码中文字幕| 亚洲第一av免费看| 丝袜喷水一区| 91成人精品电影| 国产亚洲精品第一综合不卡| 99久久人妻综合| 丝瓜视频免费看黄片| 亚洲伊人色综图| 久久人人97超碰香蕉20202| 自线自在国产av| 999久久久国产精品视频| 国产精品蜜桃在线观看| 晚上一个人看的免费电影| 亚洲av综合色区一区| 99热全是精品| 99热全是精品| 超碰97精品在线观看| 欧美亚洲日本最大视频资源| 亚洲国产欧美日韩在线播放| 女的被弄到高潮叫床怎么办| 老司机靠b影院| 国产黄色视频一区二区在线观看| 国产乱人偷精品视频| 99久久人妻综合| 亚洲精品国产区一区二| 国产无遮挡羞羞视频在线观看| 亚洲熟女毛片儿| 久久久精品区二区三区| 激情视频va一区二区三区| 999精品在线视频| 在线观看www视频免费| 婷婷成人精品国产| 欧美亚洲日本最大视频资源| 国产精品久久久久久精品古装| 国产精品国产三级国产专区5o| 国产无遮挡羞羞视频在线观看| 性色av一级| 亚洲国产欧美日韩在线播放| 久久久久久人人人人人| 午夜精品国产一区二区电影| 亚洲图色成人| 精品亚洲乱码少妇综合久久| 在线观看人妻少妇| 国产又爽黄色视频| 巨乳人妻的诱惑在线观看| 国产乱人偷精品视频| 久久久久精品性色| 天天躁日日躁夜夜躁夜夜| 成人18禁高潮啪啪吃奶动态图| 国产精品一区二区在线不卡| 人人妻人人添人人爽欧美一区卜| 最近中文字幕高清免费大全6| 9191精品国产免费久久| 亚洲一区中文字幕在线| 亚洲欧美激情在线| 精品人妻一区二区三区麻豆| 国产精品麻豆人妻色哟哟久久| 日韩电影二区| 少妇被粗大的猛进出69影院| 高清视频免费观看一区二区| 女人久久www免费人成看片| 国产精品av久久久久免费| 亚洲国产精品一区二区三区在线| 国产野战对白在线观看| 丝袜喷水一区| 考比视频在线观看| 日本黄色日本黄色录像| 免费高清在线观看日韩| 欧美精品一区二区大全| 街头女战士在线观看网站| 免费不卡黄色视频| 18禁裸乳无遮挡动漫免费视频| 国产精品免费大片| 国产高清国产精品国产三级| 久久99热这里只频精品6学生| 嫩草影院入口| 超碰成人久久| 亚洲精品第二区| 999久久久国产精品视频| 国产亚洲av片在线观看秒播厂| 国产精品香港三级国产av潘金莲 | 亚洲专区中文字幕在线 | 欧美日韩一级在线毛片| 国产男人的电影天堂91| 久久久久精品国产欧美久久久 | 国产精品免费视频内射| 性色av一级| 97在线人人人人妻| 亚洲人成网站在线观看播放| 亚洲精品美女久久av网站| 免费观看av网站的网址| 18在线观看网站| 亚洲美女视频黄频| 国产精品 欧美亚洲| 中文字幕av电影在线播放| 一本大道久久a久久精品| 在线免费观看不下载黄p国产| 国产精品久久久久久精品古装| av在线app专区| 熟女av电影| 夫妻性生交免费视频一级片| 天天操日日干夜夜撸| 男人操女人黄网站| 日韩中文字幕视频在线看片| 久久热在线av| 欧美人与性动交α欧美软件| 99热全是精品| av线在线观看网站| 久久精品aⅴ一区二区三区四区| 十八禁网站网址无遮挡| 成年美女黄网站色视频大全免费| 十八禁人妻一区二区| 视频区图区小说| 午夜久久久在线观看| 国产成人免费无遮挡视频| 丝袜脚勾引网站| 中文字幕制服av| 天天影视国产精品| 国产毛片在线视频| 免费观看人在逋| 久久天躁狠狠躁夜夜2o2o | 国产有黄有色有爽视频| 欧美精品一区二区大全| 免费在线观看黄色视频的| 免费观看性生交大片5| 欧美日韩精品网址| 亚洲精品中文字幕在线视频| 日韩精品免费视频一区二区三区| 亚洲中文av在线| 久久免费观看电影| 久久久精品国产亚洲av高清涩受| 久久久精品94久久精品| 丝袜美足系列| 亚洲五月色婷婷综合| 国产爽快片一区二区三区| 精品一区二区三卡| 伦理电影大哥的女人| 日韩av不卡免费在线播放| 久久这里只有精品19| 制服人妻中文乱码| 久久久久网色| 久久久久久人妻| 热99久久久久精品小说推荐| 国产 一区精品| 成人国语在线视频| 久久青草综合色| 97人妻天天添夜夜摸| 男人舔女人的私密视频| 成人午夜精彩视频在线观看| 午夜福利在线免费观看网站| 18在线观看网站| 可以免费在线观看a视频的电影网站 | 久久午夜综合久久蜜桃| 91精品伊人久久大香线蕉| 少妇 在线观看| 国产免费又黄又爽又色| 女性生殖器流出的白浆| 久久99精品国语久久久| 国产亚洲欧美精品永久| 亚洲成人手机| 天堂俺去俺来也www色官网| 欧美av亚洲av综合av国产av | 国产激情久久老熟女| 丝袜美足系列| 如何舔出高潮| 亚洲,欧美,日韩| 国产免费福利视频在线观看| 久久久国产一区二区| 色播在线永久视频| 777久久人妻少妇嫩草av网站| 国产乱人偷精品视频| 精品少妇黑人巨大在线播放| 欧美精品人与动牲交sv欧美| 一区二区av电影网| 久久人人爽av亚洲精品天堂| 九色亚洲精品在线播放| avwww免费| 一级毛片黄色毛片免费观看视频| 成人黄色视频免费在线看| 高清av免费在线| 午夜福利影视在线免费观看| 久久久久久人人人人人| 欧美激情 高清一区二区三区| 中文天堂在线官网| 一级毛片 在线播放| 啦啦啦视频在线资源免费观看| 日韩一区二区三区影片| 久久精品aⅴ一区二区三区四区| 精品国产一区二区三区久久久樱花| 日本色播在线视频| 国产片内射在线| 高清欧美精品videossex| 国产深夜福利视频在线观看| 久久久久精品人妻al黑| 亚洲精品日韩在线中文字幕| 中文精品一卡2卡3卡4更新| 夜夜骑夜夜射夜夜干| 久久人人97超碰香蕉20202| 国产在线免费精品| 成人手机av| 黄片播放在线免费| 亚洲欧美激情在线| 亚洲第一青青草原| 日韩av在线免费看完整版不卡| 我的亚洲天堂| 麻豆精品久久久久久蜜桃| 人妻一区二区av| 18在线观看网站| 精品一区在线观看国产| 丰满饥渴人妻一区二区三| av视频免费观看在线观看| 在线天堂最新版资源| 亚洲av福利一区| 免费在线观看视频国产中文字幕亚洲 | 亚洲国产精品999| 免费高清在线观看视频在线观看| 国产一区二区 视频在线| 黄色视频在线播放观看不卡| 国产精品国产av在线观看| 久久性视频一级片| 黄色怎么调成土黄色| 国产一级毛片在线| 久久女婷五月综合色啪小说| 日日撸夜夜添| 91精品伊人久久大香线蕉| 天天躁夜夜躁狠狠躁躁| 女人久久www免费人成看片| 少妇人妻 视频| 国产成人a∨麻豆精品| 男人舔女人的私密视频| 日韩熟女老妇一区二区性免费视频| 欧美精品一区二区免费开放| 欧美日韩视频高清一区二区三区二| 国产精品一区二区精品视频观看| 青春草视频在线免费观看| 久久人妻熟女aⅴ| 狂野欧美激情性xxxx| 亚洲自偷自拍图片 自拍| 国产又爽黄色视频| 亚洲一码二码三码区别大吗| 国产一区亚洲一区在线观看| 美国免费a级毛片| 国产xxxxx性猛交| 这个男人来自地球电影免费观看 | 亚洲五月色婷婷综合| 欧美黑人精品巨大| 国产精品 国内视频| 欧美精品av麻豆av| 亚洲熟女毛片儿| 亚洲久久久国产精品| 又大又黄又爽视频免费| 亚洲婷婷狠狠爱综合网| 精品亚洲乱码少妇综合久久| 国产欧美亚洲国产| 一区在线观看完整版| 久久精品久久精品一区二区三区| 韩国高清视频一区二区三区| 国产不卡av网站在线观看| 国产精品偷伦视频观看了| 永久免费av网站大全| 精品国产乱码久久久久久男人| 一级毛片 在线播放| av卡一久久| 高清在线视频一区二区三区| 免费观看人在逋| 国产精品三级大全| 久久精品久久久久久久性| 国产伦理片在线播放av一区| 美女主播在线视频| 亚洲国产最新在线播放| 天天影视国产精品| 午夜福利网站1000一区二区三区| 久久久精品94久久精品| 免费看av在线观看网站| 99国产精品免费福利视频| 99热国产这里只有精品6| 欧美乱码精品一区二区三区| tube8黄色片| 久久天躁狠狠躁夜夜2o2o | 国产国语露脸激情在线看| 亚洲国产最新在线播放| 人妻 亚洲 视频| 国产精品国产三级专区第一集| 日本爱情动作片www.在线观看| 国产爽快片一区二区三区| av.在线天堂| 国产精品欧美亚洲77777| 晚上一个人看的免费电影| 精品国产露脸久久av麻豆| 美女扒开内裤让男人捅视频| 韩国av在线不卡| 亚洲中文av在线| 国产一区亚洲一区在线观看| 欧美成人午夜精品| 亚洲av综合色区一区| 国产精品一区二区在线观看99| 亚洲综合色网址| 免费看av在线观看网站| 叶爱在线成人免费视频播放| 欧美少妇被猛烈插入视频| 老司机靠b影院| 中文字幕另类日韩欧美亚洲嫩草| 国产成人欧美在线观看 | 成年美女黄网站色视频大全免费| 搡老乐熟女国产| 久久久久久免费高清国产稀缺| 99re6热这里在线精品视频| e午夜精品久久久久久久| 在现免费观看毛片| 欧美黑人欧美精品刺激| 亚洲国产欧美在线一区| 黄网站色视频无遮挡免费观看| 香蕉丝袜av| 久久人人爽av亚洲精品天堂| 一个人免费看片子| 大香蕉久久成人网| 中文字幕精品免费在线观看视频| 18禁国产床啪视频网站| 视频区图区小说| a级毛片黄视频| 亚洲精品美女久久av网站| 国产一区二区在线观看av| 国产精品久久久久久精品电影小说| 久久亚洲国产成人精品v| 亚洲av男天堂| 成人亚洲精品一区在线观看| 丝袜人妻中文字幕| 欧美日韩国产mv在线观看视频| 在线天堂中文资源库| 2018国产大陆天天弄谢| 晚上一个人看的免费电影| 制服丝袜香蕉在线| 国产av国产精品国产| 精品一区二区三卡| 免费观看a级毛片全部| 下体分泌物呈黄色| 波野结衣二区三区在线| av又黄又爽大尺度在线免费看| 久久久精品区二区三区| 欧美日韩一区二区视频在线观看视频在线| 青春草亚洲视频在线观看| 九九爱精品视频在线观看| av网站免费在线观看视频| 91老司机精品| 如何舔出高潮| 99re6热这里在线精品视频| 777米奇影视久久| 亚洲国产欧美日韩在线播放| 母亲3免费完整高清在线观看| 如何舔出高潮| 男女免费视频国产| 国产精品女同一区二区软件| 国产精品麻豆人妻色哟哟久久| 视频在线观看一区二区三区| 99九九在线精品视频| 亚洲精品一二三| av片东京热男人的天堂| 国产精品久久久久久久久免| 制服丝袜香蕉在线| 我的亚洲天堂| 国产成人系列免费观看| 亚洲精品,欧美精品| 一本大道久久a久久精品| 女性生殖器流出的白浆| 国产精品女同一区二区软件| 国产色婷婷99| 亚洲国产av影院在线观看| 丝袜在线中文字幕| 亚洲av国产av综合av卡| 在现免费观看毛片| 欧美最新免费一区二区三区| 亚洲情色 制服丝袜| 午夜福利乱码中文字幕| 亚洲av成人精品一二三区| 在线亚洲精品国产二区图片欧美| 欧美精品一区二区大全| 丝袜美足系列| 亚洲男人天堂网一区| 男女边吃奶边做爰视频| 久久狼人影院| 菩萨蛮人人尽说江南好唐韦庄| 建设人人有责人人尽责人人享有的| 人体艺术视频欧美日本| 久久久久久久久免费视频了| 国产淫语在线视频| av片东京热男人的天堂| 人人妻,人人澡人人爽秒播 | 热re99久久精品国产66热6| 波野结衣二区三区在线| 成人漫画全彩无遮挡| 国产淫语在线视频| 18禁裸乳无遮挡动漫免费视频| 晚上一个人看的免费电影| 午夜av观看不卡| 九色亚洲精品在线播放| 肉色欧美久久久久久久蜜桃| 欧美日韩视频高清一区二区三区二| 国产精品久久久人人做人人爽| 亚洲免费av在线视频| 国产深夜福利视频在线观看| 人人妻人人澡人人看| 美女视频免费永久观看网站| 国产精品 国内视频| 亚洲欧美清纯卡通| 欧美最新免费一区二区三区| 日韩成人av中文字幕在线观看| 亚洲成av片中文字幕在线观看| 伦理电影免费视频| 免费少妇av软件| 成人国产av品久久久| 亚洲成人手机| 捣出白浆h1v1| 午夜免费鲁丝| 欧美老熟妇乱子伦牲交| 大片电影免费在线观看免费| 黑人巨大精品欧美一区二区蜜桃| 成人黄色视频免费在线看| 亚洲精品成人av观看孕妇| 大话2 男鬼变身卡| 人妻人人澡人人爽人人| 亚洲av日韩精品久久久久久密 | 国产成人精品无人区| 天天躁日日躁夜夜躁夜夜| 人人妻人人澡人人爽人人夜夜| 午夜免费鲁丝| 国产一卡二卡三卡精品 | 我要看黄色一级片免费的| 中文字幕人妻熟女乱码| 欧美成人午夜精品| 男女边摸边吃奶| av女优亚洲男人天堂| 成年美女黄网站色视频大全免费| 一本久久精品| 在线天堂最新版资源| 少妇猛男粗大的猛烈进出视频| 久久久久精品人妻al黑| 精品第一国产精品| 亚洲久久久国产精品| 亚洲av电影在线观看一区二区三区| 免费看不卡的av| 伦理电影免费视频| 国产又色又爽无遮挡免| 国产视频首页在线观看| 午夜免费观看性视频| 美女扒开内裤让男人捅视频| 巨乳人妻的诱惑在线观看| 色吧在线观看| 新久久久久国产一级毛片| 国产精品秋霞免费鲁丝片| 亚洲精品一区蜜桃| 亚洲国产av影院在线观看| 久久久久久久国产电影| 日韩不卡一区二区三区视频在线| 狠狠精品人妻久久久久久综合| 欧美日韩一级在线毛片| 久久久久久久大尺度免费视频| 看非洲黑人一级黄片| 老熟女久久久| 人人妻人人澡人人看| 一级a爱视频在线免费观看| 中文字幕另类日韩欧美亚洲嫩草| 日韩欧美精品免费久久| 国产高清不卡午夜福利| 亚洲国产最新在线播放| 午夜福利在线免费观看网站| 中文字幕高清在线视频| videosex国产| 亚洲国产最新在线播放| 夫妻性生交免费视频一级片| 久久精品国产综合久久久| 天天添夜夜摸| 欧美日韩国产mv在线观看视频| 大陆偷拍与自拍| 中文字幕精品免费在线观看视频| 狠狠婷婷综合久久久久久88av| 亚洲四区av| 亚洲av国产av综合av卡| 亚洲欧美激情在线| 久久性视频一级片| 飞空精品影院首页| 精品国产乱码久久久久久小说| 亚洲美女搞黄在线观看| 日日撸夜夜添| 人妻 亚洲 视频|