• <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.

    国产伦精品一区二区三区四那| 18禁在线播放成人免费| 国产成人freesex在线| 亚洲国产色片| 久久中文看片网| 黄色配什么色好看| av在线老鸭窝| 99久久九九国产精品国产免费| 日韩av不卡免费在线播放| 99热网站在线观看| 亚洲av一区综合| av女优亚洲男人天堂| 九九爱精品视频在线观看| 少妇的逼水好多| 久久久久久久久久久免费av| 边亲边吃奶的免费视频| 国产一级毛片七仙女欲春2| 中文欧美无线码| 91久久精品电影网| 99热网站在线观看| 免费人成视频x8x8入口观看| 国产日本99.免费观看| 亚洲在线自拍视频| 小说图片视频综合网站| 精华霜和精华液先用哪个| 亚洲欧洲日产国产| 久久九九热精品免费| 亚洲av电影不卡..在线观看| 少妇的逼水好多| 在线观看av片永久免费下载| 99在线人妻在线中文字幕| 久久久久九九精品影院| 别揉我奶头 嗯啊视频| 欧美日韩国产亚洲二区| 国产三级在线视频| 偷拍熟女少妇极品色| 三级毛片av免费| 国产精品电影一区二区三区| 99久久久亚洲精品蜜臀av| 亚洲自偷自拍三级| 最后的刺客免费高清国语| 久久精品国产亚洲av天美| 少妇熟女aⅴ在线视频| 人人妻人人澡欧美一区二区| 中国国产av一级| 国产麻豆成人av免费视频| 欧美潮喷喷水| 91久久精品电影网| 国产一区二区亚洲精品在线观看| 午夜福利视频1000在线观看| 国产亚洲av嫩草精品影院| 一级黄片播放器| 亚洲人成网站高清观看| 亚洲欧美精品自产自拍| 国产精品一区二区三区四区久久| 日韩亚洲欧美综合| 国产日韩欧美在线精品| 精品一区二区三区人妻视频| 久久人人爽人人爽人人片va| 麻豆成人午夜福利视频| 亚洲丝袜综合中文字幕| 精品人妻偷拍中文字幕| 色综合站精品国产| 午夜激情福利司机影院| 久久这里有精品视频免费| 日本免费一区二区三区高清不卡| 26uuu在线亚洲综合色| 久久婷婷人人爽人人干人人爱| 美女国产视频在线观看| avwww免费| 欧美日韩精品成人综合77777| 成人国产麻豆网| 在线观看一区二区三区| 久久久午夜欧美精品| 欧美+亚洲+日韩+国产| 日韩一本色道免费dvd| 日韩精品青青久久久久久| 日本五十路高清| 男女啪啪激烈高潮av片| 国产老妇女一区| 国产精品不卡视频一区二区| 精品久久久噜噜| 波多野结衣巨乳人妻| 看黄色毛片网站| 看免费成人av毛片| 亚洲成人av在线免费| 亚洲精品日韩在线中文字幕 | 人妻少妇偷人精品九色| 国产成人午夜福利电影在线观看| 亚洲真实伦在线观看| 国产精品精品国产色婷婷| 麻豆av噜噜一区二区三区| 国产亚洲精品av在线| 久久人人精品亚洲av| 亚洲av一区综合| 国产精品蜜桃在线观看 | 丰满乱子伦码专区| 人妻系列 视频| 中国美白少妇内射xxxbb| 自拍偷自拍亚洲精品老妇| 少妇裸体淫交视频免费看高清| 亚洲精品456在线播放app| 欧美激情久久久久久爽电影| 免费av毛片视频| 禁无遮挡网站| 成人午夜精彩视频在线观看| 神马国产精品三级电影在线观看| 偷拍熟女少妇极品色| 如何舔出高潮| 亚洲成人av在线免费| 国产成年人精品一区二区| 麻豆国产av国片精品| 亚洲人成网站在线观看播放| 黄片无遮挡物在线观看| 一区二区三区四区激情视频 | 亚洲内射少妇av| 99热6这里只有精品| 久久久久免费精品人妻一区二区| 日韩中字成人| 亚洲精品国产av成人精品| 亚洲在线观看片| 久久久久久九九精品二区国产| 淫秽高清视频在线观看| 99久国产av精品| 在线观看66精品国产| 国产精品国产三级国产av玫瑰| 亚洲av中文字字幕乱码综合| 久久精品91蜜桃| 亚洲国产欧美在线一区| 免费在线观看成人毛片| 伊人久久精品亚洲午夜| 美女国产视频在线观看| 国产美女午夜福利| 97超碰精品成人国产| 久久国产乱子免费精品| 97热精品久久久久久| 精品久久久久久久久久免费视频| ponron亚洲| 看非洲黑人一级黄片| 久久久久久久午夜电影| 干丝袜人妻中文字幕| 欧美激情久久久久久爽电影| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 精品少妇黑人巨大在线播放 | 精品久久久久久久末码| 观看免费一级毛片| 国产三级中文精品| 亚洲七黄色美女视频| 亚洲国产精品国产精品| 边亲边吃奶的免费视频| 欧美日韩乱码在线| 国产综合懂色| 又粗又爽又猛毛片免费看| 一个人免费在线观看电影| 日本三级黄在线观看| 不卡一级毛片| 成人漫画全彩无遮挡| 亚洲av电影不卡..在线观看| 国内精品美女久久久久久| 国产精品爽爽va在线观看网站| 久久草成人影院| 精品不卡国产一区二区三区| 亚洲无线在线观看| 精品日产1卡2卡| 99视频精品全部免费 在线| 精品久久久久久久久av| 欧美成人精品欧美一级黄| 女同久久另类99精品国产91| 国产精品爽爽va在线观看网站| 午夜福利成人在线免费观看| 国产女主播在线喷水免费视频网站 | 美女被艹到高潮喷水动态| 国产真实伦视频高清在线观看| 日韩av在线大香蕉| 亚洲欧美日韩高清在线视频| 久久人人精品亚洲av| videossex国产| 欧美潮喷喷水| 精华霜和精华液先用哪个| 免费观看的影片在线观看| 久久久欧美国产精品| 久久精品国产自在天天线| 亚洲欧美清纯卡通| 国产伦理片在线播放av一区 | 日本免费a在线| 插阴视频在线观看视频| 成年免费大片在线观看| 国产av一区在线观看免费| 午夜福利成人在线免费观看| 中出人妻视频一区二区| 一个人看视频在线观看www免费| 日本五十路高清| 色哟哟·www| 欧美极品一区二区三区四区| 色哟哟哟哟哟哟| 免费观看人在逋| 真实男女啪啪啪动态图| 日本与韩国留学比较| 欧美成人一区二区免费高清观看| 在线国产一区二区在线| 国产精品久久久久久精品电影| 日韩 亚洲 欧美在线| 国产精品.久久久| 99在线人妻在线中文字幕| 国产v大片淫在线免费观看| 亚洲电影在线观看av| 国产精品女同一区二区软件| 老师上课跳d突然被开到最大视频| 亚洲精品自拍成人| 97热精品久久久久久| 免费av观看视频| 国内揄拍国产精品人妻在线| 欧美xxxx性猛交bbbb| 伦理电影大哥的女人| 少妇的逼好多水| 校园春色视频在线观看| 亚洲激情五月婷婷啪啪| 亚洲五月天丁香| 欧美变态另类bdsm刘玥| 丰满的人妻完整版| 成年免费大片在线观看| av在线蜜桃| а√天堂www在线а√下载| 一本精品99久久精品77| 深夜a级毛片| 一级毛片我不卡| 91精品一卡2卡3卡4卡| 久久久久久伊人网av| 久久久久性生活片| 亚洲国产高清在线一区二区三| av在线天堂中文字幕| 国产大屁股一区二区在线视频| 亚洲精品日韩在线中文字幕 | 不卡一级毛片| 99久久无色码亚洲精品果冻| 欧美高清成人免费视频www| 精品一区二区免费观看| 女人被狂操c到高潮| 不卡一级毛片| 亚洲精品影视一区二区三区av| 小蜜桃在线观看免费完整版高清| 国产精品久久久久久精品电影| 麻豆精品久久久久久蜜桃| 在线天堂最新版资源| 青青草视频在线视频观看| 精品久久久噜噜| 3wmmmm亚洲av在线观看| 国产黄色视频一区二区在线观看 | 菩萨蛮人人尽说江南好唐韦庄 | 亚洲美女搞黄在线观看| 日本与韩国留学比较| 男人和女人高潮做爰伦理| 不卡视频在线观看欧美| 亚洲国产欧美人成| 午夜久久久久精精品| 亚洲成人精品中文字幕电影| 成人性生交大片免费视频hd| 欧美最新免费一区二区三区| 超碰av人人做人人爽久久| 黄色欧美视频在线观看| 搡女人真爽免费视频火全软件| 国产中年淑女户外野战色| 免费观看在线日韩| 六月丁香七月| 日本熟妇午夜| 国产视频首页在线观看| 大又大粗又爽又黄少妇毛片口| 男女下面进入的视频免费午夜| 最近的中文字幕免费完整| 在线观看免费视频日本深夜| 内地一区二区视频在线| 99视频精品全部免费 在线| 校园人妻丝袜中文字幕| 欧美潮喷喷水| 亚洲av免费在线观看| 中文字幕制服av| 亚洲精品日韩在线中文字幕 | 久久人人精品亚洲av| 国产淫片久久久久久久久| 亚洲经典国产精华液单| 搡老妇女老女人老熟妇| 中文字幕av在线有码专区| 高清日韩中文字幕在线| 中文字幕制服av| 日产精品乱码卡一卡2卡三| 国产精品嫩草影院av在线观看| 国语自产精品视频在线第100页| 舔av片在线| 国产色婷婷99| 国产蜜桃级精品一区二区三区| 婷婷色av中文字幕| 久久久午夜欧美精品| 麻豆国产av国片精品| 亚洲成人av在线免费| 91在线精品国自产拍蜜月| 国产在线男女| 少妇的逼好多水| 中文在线观看免费www的网站| 少妇被粗大猛烈的视频| 日本黄色片子视频| 欧美+亚洲+日韩+国产| 男人和女人高潮做爰伦理| 欧美最新免费一区二区三区| 丝袜美腿在线中文| 亚洲欧美日韩东京热| 久久精品综合一区二区三区| 日韩欧美国产在线观看| 国产高清激情床上av| 久久精品国产99精品国产亚洲性色| 天天躁夜夜躁狠狠久久av| 免费观看的影片在线观看| 亚洲经典国产精华液单| 91久久精品国产一区二区成人| 亚洲精华国产精华液的使用体验 | 特级一级黄色大片| 嫩草影院精品99| 亚洲丝袜综合中文字幕| 日韩三级伦理在线观看| 日韩欧美国产在线观看| 国产成人91sexporn| 亚洲av不卡在线观看| 久久99精品国语久久久| 99热6这里只有精品| 男的添女的下面高潮视频| 亚洲va在线va天堂va国产| 国产黄片美女视频| 99国产精品一区二区蜜桃av| 成人高潮视频无遮挡免费网站| 中文字幕熟女人妻在线| 日韩欧美在线乱码| a级毛片免费高清观看在线播放| 精品久久久久久久久久免费视频| 99热只有精品国产| 在线播放国产精品三级| 亚洲精品自拍成人| 国产一级毛片七仙女欲春2| 久久久久国产网址| 大型黄色视频在线免费观看| 五月玫瑰六月丁香| 久久久久性生活片| 一个人看的www免费观看视频| 久久久精品94久久精品| 国产亚洲精品久久久com| 人人妻人人看人人澡| 变态另类丝袜制服| 啦啦啦韩国在线观看视频| 亚洲av免费在线观看| 国产亚洲av片在线观看秒播厂 | 免费人成在线观看视频色| 亚洲精品色激情综合| 国产大屁股一区二区在线视频| 日韩欧美一区二区三区在线观看| 高清毛片免费看| 国产 一区精品| 成人鲁丝片一二三区免费| 波多野结衣高清无吗| 99国产精品一区二区蜜桃av| 岛国在线免费视频观看| 日本免费a在线| 午夜久久久久精精品| 亚洲激情五月婷婷啪啪| 一级毛片久久久久久久久女| 国产黄色视频一区二区在线观看 | 秋霞在线观看毛片| 国产日韩欧美在线精品| 精品一区二区三区人妻视频| 成人漫画全彩无遮挡| 九九爱精品视频在线观看| 国产午夜福利久久久久久| a级毛片免费高清观看在线播放| 精品无人区乱码1区二区| 在线国产一区二区在线| 国产精品综合久久久久久久免费| 亚洲自偷自拍三级| 麻豆av噜噜一区二区三区| 亚洲无线观看免费| 久久久久久伊人网av| 亚洲aⅴ乱码一区二区在线播放| 成人一区二区视频在线观看| 国产精品一及| 女的被弄到高潮叫床怎么办| 人妻夜夜爽99麻豆av| 久久午夜亚洲精品久久| 日韩欧美 国产精品| 日韩亚洲欧美综合| 日韩av不卡免费在线播放| 在线免费观看不下载黄p国产| 亚洲精品日韩在线中文字幕 | 菩萨蛮人人尽说江南好唐韦庄 | 人人妻人人看人人澡| 激情 狠狠 欧美| 91狼人影院| 观看美女的网站| 日韩国内少妇激情av| 在线观看美女被高潮喷水网站| av在线亚洲专区| 一夜夜www| 欧美不卡视频在线免费观看| 国产精品久久视频播放| 亚洲国产高清在线一区二区三| 国内精品宾馆在线| 欧美+亚洲+日韩+国产| 麻豆精品久久久久久蜜桃| 亚洲成人精品中文字幕电影| 成年av动漫网址| 久久久成人免费电影| 国产综合懂色| 精品久久久久久久久av| 亚洲国产日韩欧美精品在线观看| av在线播放精品| 国产精品.久久久| 亚洲精品久久国产高清桃花| 亚洲欧美精品自产自拍| 三级毛片av免费| 亚洲乱码一区二区免费版| 欧美成人精品欧美一级黄| 久久亚洲精品不卡| 久久久久久久亚洲中文字幕| 有码 亚洲区| 国产成人午夜福利电影在线观看| 国产黄片视频在线免费观看| 一级黄色大片毛片| 国内精品宾馆在线| 神马国产精品三级电影在线观看| 国产精品蜜桃在线观看 | 男人狂女人下面高潮的视频| 亚洲欧洲国产日韩| 国产精品麻豆人妻色哟哟久久 | 久久久成人免费电影| 男女那种视频在线观看| 长腿黑丝高跟| 中文资源天堂在线| av视频在线观看入口| 国产高清不卡午夜福利| av女优亚洲男人天堂| 午夜福利在线在线| 精品日产1卡2卡| 久久精品国产亚洲网站| 亚洲欧美日韩东京热| 日韩强制内射视频| 丝袜喷水一区| 在线观看av片永久免费下载| 最后的刺客免费高清国语| 天天一区二区日本电影三级| 一个人观看的视频www高清免费观看| 人体艺术视频欧美日本| 成人性生交大片免费视频hd| 黄片无遮挡物在线观看| 国产精品永久免费网站| 一边亲一边摸免费视频| 中文资源天堂在线| 亚洲精品影视一区二区三区av| av.在线天堂| 日韩 亚洲 欧美在线| 黄片无遮挡物在线观看| 国产乱人偷精品视频| 长腿黑丝高跟| 国产精品99久久久久久久久| 夜夜爽天天搞| 悠悠久久av| 亚洲成a人片在线一区二区| 婷婷色综合大香蕉| 乱码一卡2卡4卡精品| 亚洲精品粉嫩美女一区| 特级一级黄色大片| 99riav亚洲国产免费| 久久久久久九九精品二区国产| 美女高潮的动态| 一区福利在线观看| 欧美一级a爱片免费观看看| 亚洲最大成人av| 国产毛片a区久久久久| 亚洲欧美清纯卡通| 久久久午夜欧美精品| 午夜精品国产一区二区电影 | 免费看日本二区| 久久久久久伊人网av| 日韩一区二区视频免费看| 欧美日韩乱码在线| 全区人妻精品视频| 国产日本99.免费观看| 一区福利在线观看| 国产精品电影一区二区三区| 最近中文字幕高清免费大全6| 免费av不卡在线播放| 久久久久网色| 国产精品无大码| 精品人妻一区二区三区麻豆| a级毛片免费高清观看在线播放| 久久鲁丝午夜福利片| 国产精华一区二区三区| 日本五十路高清| 高清毛片免费看| 一边亲一边摸免费视频| 久久鲁丝午夜福利片| 亚洲最大成人手机在线| 日产精品乱码卡一卡2卡三| 亚洲精品乱码久久久v下载方式| 国产成人精品婷婷| 内射极品少妇av片p| 日韩欧美精品免费久久| 欧美又色又爽又黄视频| 蜜桃久久精品国产亚洲av| 天堂√8在线中文| 内射极品少妇av片p| 亚洲精品国产av成人精品| 欧美又色又爽又黄视频| 国产av一区在线观看免费| 国国产精品蜜臀av免费| 亚洲色图av天堂| 成人av在线播放网站| 婷婷色综合大香蕉| 联通29元200g的流量卡| 岛国在线免费视频观看| 两个人的视频大全免费| 91精品国产九色| 美女黄网站色视频| 亚洲成av人片在线播放无| 夜夜爽天天搞| 亚洲激情五月婷婷啪啪| 国内久久婷婷六月综合欲色啪| 国产精品日韩av在线免费观看| 不卡视频在线观看欧美| 亚洲国产精品成人综合色| h日本视频在线播放| 国产亚洲精品久久久com| 久久99蜜桃精品久久| 精品日产1卡2卡| 草草在线视频免费看| 最近最新中文字幕大全电影3| 久久久久久大精品| 夫妻性生交免费视频一级片| 美女 人体艺术 gogo| 男女做爰动态图高潮gif福利片| 成人一区二区视频在线观看| 卡戴珊不雅视频在线播放| 久久久久久九九精品二区国产| 成人亚洲精品av一区二区| 国产成人福利小说| 久久久久久久久久久免费av| 99久国产av精品国产电影| 亚洲精品日韩在线中文字幕 | a级毛片免费高清观看在线播放| 网址你懂的国产日韩在线| 国产麻豆成人av免费视频| 男女做爰动态图高潮gif福利片| 赤兔流量卡办理| 亚洲婷婷狠狠爱综合网| 亚洲精品日韩在线中文字幕 | 国产精品爽爽va在线观看网站| av福利片在线观看| 久久久久性生活片| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 久久久久久久久大av| 日韩在线高清观看一区二区三区| 日韩高清综合在线| 日韩欧美精品v在线| 国产成人精品一,二区 | 亚洲成av人片在线播放无| 插阴视频在线观看视频| 又粗又爽又猛毛片免费看| 尤物成人国产欧美一区二区三区| 中国国产av一级| 久久精品91蜜桃| 如何舔出高潮| 久久久久久国产a免费观看| 夫妻性生交免费视频一级片| 亚洲成a人片在线一区二区| 日韩三级伦理在线观看| 人妻制服诱惑在线中文字幕| 国产美女午夜福利| 国产成人91sexporn| 日日摸夜夜添夜夜添av毛片| 男人的好看免费观看在线视频| 村上凉子中文字幕在线| 国产在线精品亚洲第一网站| 国产高清有码在线观看视频| 国产爱豆传媒在线观看| av黄色大香蕉| 国产黄色小视频在线观看| 人妻夜夜爽99麻豆av| 蜜桃久久精品国产亚洲av| 亚洲最大成人av| 午夜福利在线观看免费完整高清在 | 亚洲第一区二区三区不卡| 晚上一个人看的免费电影| 色综合站精品国产| 国产高清不卡午夜福利| 干丝袜人妻中文字幕| 少妇的逼好多水| 一级黄色大片毛片| 久久精品国产99精品国产亚洲性色| 亚洲内射少妇av| 欧美最新免费一区二区三区| 成人鲁丝片一二三区免费| 一边亲一边摸免费视频| 成人亚洲欧美一区二区av| 男人狂女人下面高潮的视频| 亚洲国产欧美人成| 熟女电影av网| 人妻久久中文字幕网| or卡值多少钱| 色哟哟哟哟哟哟| 直男gayav资源| 午夜亚洲福利在线播放| 日韩人妻高清精品专区| 午夜久久久久精精品| av在线播放精品| 国产午夜福利久久久久久| 夜夜夜夜夜久久久久| 又粗又硬又长又爽又黄的视频 | 国产黄片视频在线免费观看| av天堂中文字幕网| 免费看日本二区| 国产视频首页在线观看|