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

    圓柱繞流氣動噪聲數(shù)值分析

    2016-09-07 06:24:32顧信忠李舜酩
    聲學技術 2016年2期
    關鍵詞:渦的測量點聲壓級

    顧信忠,李舜酩

    ?

    圓柱繞流氣動噪聲數(shù)值分析

    顧信忠1,2,李舜酩1

    (1. 南京航空航天大學能源與動力學院,江蘇南京210016; 2. 南京航空航天大學金城學院車輛系,江蘇南京211156)

    為了快速預測剛性圓柱繞流的氣動噪聲,研究了一種將離散渦方法(DVM)和渦聲理論結合起來計算低馬赫數(shù)、高雷諾數(shù)流場氣動噪聲的方法。首先用Oseen粘性渦模型改進了離散渦方法并模擬了圓柱繞流,分析結果與實際情況相符。根據(jù)流場計算的結果,應用渦聲理論進一步計算了遠場的聲壓。測點的總聲壓級與實驗值及其他數(shù)值計算結果都比較吻合。最后繪制了聲場的指向性特性曲線,表明圓柱繞流聲場明顯的偶極子特性。

    氣動噪聲;圓柱繞流;離散渦方法;渦聲理論;Oseen渦

    0 引言

    圓柱體繞流產生聲音是一種非常普遍的現(xiàn)象。產生的機理包括:流體與圓柱壁面的相互作用使圓柱產生振動,振動的圓柱排開流體就形成了單極子聲源;當旋渦從圓柱表面交替脫落時,會在圓柱表面產生正負環(huán)量,柱面附近區(qū)域因此產生周期性壓力波動,從而形成了偶極子聲源;圓柱表面邊界層中的渦旋和脫落的渦旋會因湍流應力而不斷衰減,此湍流應力就是四極子聲源[1]。

    圓柱繞流產聲問題在航空航天、航海、風工程等實際工程中非常具有代表性,因此眾多的學者對此進行過深入研究,所采用的方法大致可以分為聲類比方法和直接數(shù)值模擬方法兩大類。

    文獻[2]采用有限體積法分別求解二維、三維可壓縮平均雷諾數(shù)N-S方程,對馬赫數(shù)為0.2、雷諾數(shù)為100~5000000的圓柱繞流進行了數(shù)值分析,并用Lighthill聲類比方法獲得了遠場聲信息。文獻[3]用大渦模擬方法模擬了雷諾數(shù)為50000~90000的圓柱繞流湍流流場,提取聲源的強度和分布信息后,用FW-H方程的積分外推法計算了遠場氣動噪聲,研究了圓柱繞流氣動噪聲的聲場特性。文獻[4]則采用高精度的空間和時間離散算法直接求解N-S方程(DNS方法),計算結果可以清晰顯示聲波的產生和傳播過程而不受馬赫數(shù)和雷諾數(shù)的限制。文獻[5]同樣采用DNS方法分析了馬赫數(shù)為0.2、雷諾數(shù)為150時圓柱繞流的氣動噪聲,著重研究了渦與聲壓之間的內在關系。

    上述兩類方法都需要采用非常細密的網格和高精度算法對流場進行準確的計算,因此需要大量的存儲空間和很長的計算時間。本文則是基于離散渦方法和渦聲理論,采用Oseen粘性渦模型,研究了可快速預測繞流產生氣動噪聲的方法。

    1 流場模擬

    1.1 離散渦方法

    離散渦方法通過計算渦量,求解Navier-Stokes方程,以拉格朗日方式跟蹤脫落渦的運動軌跡,恰能滿足流場渦聲計算的需要,因此本文采用離散渦方法模擬流場。

    離散渦方法的基本思路是將流場分成有旋和無旋兩部分,有旋部分主要集中在柱面的邊界層和分離后的剪切層,并認為兩者都是薄薄的連續(xù)渦層,可用離散的渦元來模擬,并將其置于無旋流場中,通過計算離散渦元之間相互作用和演化實現(xiàn)對流場的快速模擬。當柱面的軸向尺寸遠大于其斷面尺寸時,則可用二維離散渦方法模擬流場[6-7],如圖1所示。

    1.2 渦模型選擇

    當兩個點渦接近時產生誘導速度非常大,與實際情況不符,因此很多文獻[8-9]采用Rankin渦模型進行光滑處理。但是渦聲與渦量的變化率有關,而Rankin渦的渦量是恒定的,因此本文采用的是Oseen粘性渦模型,其周向速度為[10]

    其中:為渦的環(huán)量;為離渦核中心的距離;為流體的運動粘性系數(shù);為計算時間。

    1.3 離散渦方法計算流程

    1.3.1 邊界層渦量的確定

    將圓柱面用若干個等長首尾相連的直線段來離散,直線段的端點布置適當強度的渦元,稱為渦點,直線段的中點稱為控制點。根椐Dirichet邊界條件——物面的切向速度為零,要求任意控制點處的流速與渦元產生的誘導速度相應[11],即

    其中:(s,s)為相關系數(shù),表示s處單位強度點渦在s控制點處產生的誘導速度;(s)為s渦點處渦的強度;u、v分別為來流的速度分量;β為控制點s的位矢與水平方向的夾角。寫成矩陣形式為

    (3)

    求解該線性方程組即可得到各渦點處渦元的強度。

    1.3.2 渦的脫落

    邊界層中的渦會不斷從柱面分離進入尾流,本文采用固定分離點的方式產生脫落渦,即預先在物面上確定幾個渦點產生脫落渦。根據(jù)庫塔(Kutta)條件可確定初生脫落渦的環(huán)量為[12]

    其中:u為相應渦點處邊界層外的流速;為計算的時間步長。

    初生脫落渦的平均對流速度為[13]

    1.3.3 渦的對流

    已知二維不可壓縮無粘流場中所有渦元(包括邊界層中的渦元和脫落渦)的強度和位置,則可用Biot-Savart公式計算出流場中任意點處的誘導速度為

    式中:為渦元在復平面中的位矢;xy分別為其水平和豎直分量。

    將誘導速度與來流速度疊加就可以得到流場中任意位置的實際速度為

    假設在時刻渦元的位置為()、(),則可得到渦元對流的計算公式為

    (8)

    其中:和為渦元因粘性擴散產生的運動,可采用隨機走位估算[14]:

    式中:和分別是(0,1)和(0,2π)區(qū)間內兩個相互獨立的隨機數(shù)。97

    1.3.4 渦的合并與鏡像

    實際情況下,當兩渦元之間距離足夠小時,就會因渦核被破壞而發(fā)生合并。于是Spalart提出了渦的融合機制[15],即認為當兩個渦元滿足一定條件時,就被合并成一個渦元。

    式中:0和0都是控制參數(shù);0控制圓柱附近的渦元數(shù)目;0控制流場中渦元總數(shù)。合并后新渦元的位置和環(huán)量分別為:

    , (11)

    當渦元運動到圓柱內部時,文中采用鏡像的方法將其強制移到圓柱外部。

    1.3.5 作用在圓柱上的力

    流體流經圓柱表面時會產生相互作用,圓柱上的作用力可用布拉休斯公式計算,將該作用力在豎直方向和水平方向進行分解,并除以0.5∞2進行無量綱化,即可得到圓柱上的升力和阻力系數(shù)[12]:

    式中:Г、u、v、β分別為第渦點處的渦的環(huán)量、柱面的切向速度、法向速度和方位角;∞為平均來流速度;為圓柱直徑。

    1.4 流場仿真結果

    本文以實驗驗證的模型[16]為研究對象,圓柱直徑為19 mm,空氣流速為69.2 m/s,此時雷諾數(shù)為Re≈90000,馬赫數(shù)約為0.2。圓柱面上均勻布置了40個渦點,其中5個為預定的脫落點,如圖1所示。設置計算時間步長為10-5 s,圖2是=0.02 s時的渦元分布圖,此時流場中有3600多個脫落渦元,可見從圓柱脫落的渦并沒有形成穩(wěn)定的卡門渦街,而是擴散到尾流區(qū)很寬的范圍內形成了紊流,與實際情況相符。

    升力系數(shù)和阻力系數(shù)如圖3和圖4所示,因為在計算過程中存在大量的渦元合并與鏡像,使得升力系數(shù)的變化范圍非常大,流場穩(wěn)定后升力系數(shù)在0值的上下等幅振蕩。阻力系數(shù)的時間平均值約為1.52,比實驗值1.32[16]大15.2%,比文獻[17]中的大渦模擬結果1.47大3.4%,如表1所示。說明文中采用的離散渦方法能夠較為準確地模擬圓柱繞流。

    表1 圓柱繞流平均阻力系數(shù)(Re≈90000,馬赫數(shù)≈0.2)

    2 聲場計算

    2.1 渦聲理論

    根據(jù)氣動聲學理論[18],在大雷諾數(shù)、低馬赫數(shù)無粘三維流場中,運動渦元產生的遠場聲壓為

    式中:0為流體密度;為渦量;為渦的運動速度;為聲源發(fā)射聲波的時刻;為緊致格林函數(shù);1為觀測點的位矢;2為聲源即渦的位矢。

    若柱體的軸向尺寸遠大于其斷面尺寸,可在二維中計算渦聲[18],此時緊致格林函數(shù)為

    其中:2為基爾霍夫矢量;0為聲波傳播速度;H(.)為Heaviside廣義函數(shù)。將式(15)代入式(14)即可得到二維流場的聲壓為

    (16)

    分別將=1、2代入上式,就可以得到柱體上脈動升力和阻力產生的遠場聲壓。若聲場中的柱體為圓柱,則基爾霍夫矢量為

    2.2 聲場計算結果

    2.2.1 聲壓

    為了與實驗值比較,與文獻[16]相同,在遠場取三個測量點,坐標分別為:測量點1(0.665,0)、測量點2(0,-0.665)、測量點3(0,2.432),如圖5所示,各測點的聲壓如圖6所示。

    對測點的聲壓信號進行快速傅里葉變換,得到頻域下的聲壓規(guī)律,利用式(18)計算各測點的聲壓級為

    (a)

    (b)

    圖6 各測點聲壓

    Fig.6 The time-histories of sound pressures at different testing points

    測量點聲壓級的頻譜特性曲線如圖7所示,噪聲頻率分布范圍很寬,且測點的聲壓級在750 Hz左右有明顯的峰值,與實際流場中渦從柱面脫落的周期相應。

    還可進一步計算各測點的總聲壓級:

    得到測量點1總聲壓級為118.52 dB,而實驗測量值為117 dB[16],誤差僅為1.3%,比大渦模擬的結果115.72[3]高2.4%,比115.16 dB[17]高2.9%;測量點3 總聲壓級的仿真值為113.74 dB,實驗測量值為100 dB,誤差為13.7%,比大渦模擬的結果104.52 dB[3]高8.8%,比103.89 dB[17]高9.5%,如表2所示。

    表2 圓柱繞流總聲壓級(Re≈90000,馬赫數(shù)≈0.2)

    Tabie 2 Overall sound pressure levels at different testing points (Re≈90000, Mach number≈0.2)

    研究方法測量點1 測量點3 總聲壓 誤差總聲壓 誤差 本文中方法118.52/ 113.74/ 二維大渦模擬[3]115.722.4% 104.528.8% 二維大渦模擬[17]115.162.9% 103.899.5% 實驗結果[16]1171.3% 10013.7%

    測點1、3的總聲壓級均大于實測值及大渦模擬的結果,可能產生誤差的原因包括以下幾個方面:第一,在仿真過程中用大量的離散渦元模擬邊界層和剪切層中的連續(xù)渦層,使得聲壓信號中的高頻成分偏高;第二,為了編程方便,文中將邊界層的離散渦元直接布置在物面上,而當聲源靠近物面時,對外輻射的聲壓較大;第三,在計算渦元對流時采用的是精度較低的一階差分格式,可能對計算結果存在一定的影響??傮w來看,仿真預測結果是較為可信的。

    2.2.2 指向性特性

    圖8是聲場在=0.019 s、=0.0195 s和=0.02 s時的輻射特性曲線,聲壓最大輻射值都在垂直于來流的軸線上,指向性與偶極子聲源的指向性相近,說明圓柱繞流的氣動噪聲主要是由偶極子聲源引起的。在平行于來流方向仍有一定的聲輻射,這是由作用在圓柱上的阻力周期性波動以及渦元粘性衰減的四極子聲源共同作用的結果。

    =0.019

    =0.0195

    3 結論

    (1) 采用離散渦方法對二維圓柱繞流流場的模擬結果顯示,流場圖與實際情況相符,升力系數(shù)和阻力系數(shù)的計算值都與實驗結果較為吻合,模擬結果是可信的。

    (2) 圓柱繞流產生的氣動噪聲的峰值頻率與渦的脫落頻率相應,總聲壓級與實驗值及大渦模擬仿真結果比較吻合,聲場的輻射特性與實際情況相符,聲場的計算結果是比較準確的。

    (3) 由于采用大量的離散渦元模擬邊界層和剪切層的連續(xù)渦層,獲得的聲信號中的高頻成分偏高。

    (4) 采用離散渦方法和渦聲理論相結合的方法,對聲場進行快速預測的方法是可行的。

    [1] 張強. 氣動聲學基礎[M]. 北京: 國防工業(yè)出版社, 2012.

    ZHANG Qiang. Foundation of Acoustics[M]. Beijing: National Defense Industry Press, 2012.

    [2] COX J S, Brentner K S, Rumsey C L. Computation of vortex shedding and radiated sound for a circular cylinder: subcritical to transonic Reynolds numbers[J]. Theoret. Comput. Fluid Dyn, 1988, 12(4): 233-253.

    [3] 龍雙麗, 聶宏, 許鑫. 不同雷諾數(shù)下圓柱繞流氣動噪聲數(shù)值模擬[J]. 聲學技術, 2011, 30(2): 111-116.

    LONG Shuangli, NIE Hong, XU Xin. Numerical simulation of noise induced by flow around a cylinder at different Reynolds number[J]. Technical Acoustics, 2011, 30(2): 111-116.

    [4] Mitchell B E, Lele S K, Moin P. Direct computation of the sound generated by vortex pairing in an axisymmetric jet[J]. J. Fluid Mech, 1999, 383(6): 113-142.

    [5] Inoue O, Hatakeyama N. Sound generation by a two-dimensional circular cylinder in a uniform flow[J]. Journal of Fluid Mechanics, 2002, 471(01): 285-314.

    [6] MUSTTO A A, Hirata M H, Bodstein G C R. Discrete vortex method simulation of the flow around a circular cylinder with and without rotation[J]. AIAA paper, 1998, 98-2409: 59-69.

    [7] Etoh F, Kamemoto K, Matsumoto H, et al. Numerical simulation of flow around a rotary oscillating foil with constant amplitude angle by use of the vortex method[C]//Proc. 11th Symp. on CFD, Tokyo, 1997: 385-386.

    [8] 潘巖松. 高層建筑二維流場的離散渦方法數(shù)值模擬[D]. 武漢: 華中科技大學, 2004.

    PAN Yansong. Numerical simulation of two-dimensional flow past a tall building using a discrete vortex method[D]. Wuhan: Huazhong University of Science & Technology, 2004.

    [9] 徐曉亮, 黃海明, 章梓茂. 燒蝕環(huán)境下的圓柱繞流計算模型[J]. 空氣動力學學報, 2011, 29(2): 170-176.

    XU Xiaoliang, Huang Haiming, Zhang Zimao. Numerical models for flow past a cylinder in the ablation environment[J]. Acta Aerodynam Ica Sinica, 2011, 29(2): 170-176.

    [10] Schaefer J W, Eskinazi S. An analysis of the vortex street generated in a viscous fluid[J]. J. Fluid Mech.,1959, 6(2): 241-260.

    [11] Lewis R I, Vortex element methods for fluid dynamic analysis of engineering systems[M]. Cambridge University Press, 2005.

    [12] Sarpkaya T, Shoaff R L. An inviscid model of two dimensional vortex shedding for transient and asymptotically-steady separated flow over a cylinder [J]. AIAA Paper, 1979, 79-0281: 1-9.

    [13] 童秉綱, 夏南, 李潛.物體繞流的離散渦方法[J]. 力學進展, 1985, 15(3): 318-328.

    TONG Bingang, XIA Nan, LI Qian. Discrete Vortex Method in external flow problems [J].Advances in Mechanics, 1985, 15(3): 318-328.

    [14] 陳偉, 宗智. 二維圓柱繞流的離散渦數(shù)值模擬[J]. 艦船科學技術, 2010, 32(5): 111-115.

    CHEN Wei, ZONG Zhi. Numerical simulation of two- dimensional flow around circular cylinder using discrete vortex method[J]. Ship SciencE and Technology, 2010, 32(5): 111-115.

    [15] SPALARTR S, LEONARD A, BAGANOFF D. Numerical simulation of separated flows[D]. Stanford University, California, 1983.

    [16] Revell J D, Prydz R A, Hays A P. Experimental study of aerody-namic noise vs drag relationships for circular cylinders[J]. Anz Journal of Surgery, 1978, 16(6): 889-897.

    [17] YANG Dangguo, LI Jianqiang, LIU Jun. Analysis on physical mechanism of sound generation inside cavities based on acoustic analogy method[J]. Journal of Fluid Dynamics, 2013, 3(1): 23-31.

    [18] HOWE M S. Theory of vortex sound [M]. Cambridge Texts in Applied Mathematics, 2002.

    Numerical analysis of noise induced by flow around a cylinder

    GU Xin-zhong1,2, LI Shun-ming1

    (1. College of Energy and Power Engineering,Nanjing University of Aeronautics &Astronautics, Nanjing 210016,Jiangsu, China;2. Department of Vehicle Engineering, Nanhang Jincheng College, Nanjing211156, Jiangsu,China)

    The method of combination of Discrete Vortex Method (DVM) and the theory of vortex sound are used to predict two-dimensional far field noise induced by flow around a cylinder. Firstly, the field of flow around a cylinder is simulated using the Oseen vortex model. The simulation results are in agreement with experimental data. Based on the simulation results, the characteristic of the acoustic field is studied. The overall sound pressure levels of receivers also accord wellwith both experimental data and other numerical simulation results. The polar diagrams of the aeroacoutic pressure are plotted. And the dipolar nature of the generated sound is confirmed.

    aerodynamic noise; flow around a cylinder; Discrete Vortex Method(DVM); theory of vortex sound; oseen vortex

    TB533

    A

    1000-3630(2016)-02-0095-06

    10.16300/j.cnki.1000-3630.2016.02.003

    2015-04-10;

    2015-06-10

    顧信忠(1980-), 男, 江蘇儀征人, 博士研究生, 研究方向為氣動噪聲預測與控制技術。

    顧信忠, E-mail: guxingzhong007@126.com

    猜你喜歡
    渦的測量點聲壓級
    機器噪聲平均聲壓級計算方法差異性實證研究
    電動工具(2024年1期)2024-02-29 01:40:24
    飛機部件數(shù)字化調姿定位測量點的優(yōu)選與構造算法
    航空學報(2022年5期)2022-07-04 02:24:32
    南海中尺度渦的形轉、內轉及平移運動研究
    海洋通報(2020年5期)2021-01-14 09:26:52
    一種計算消聲室聲壓級的新方法
    計量學報(2020年2期)2020-04-11 04:33:22
    淺析沖壓件測量點的規(guī)劃
    模具制造(2019年10期)2020-01-06 09:13:08
    全新DXR mkll有源揚聲器
    演藝科技(2019年4期)2019-03-30 03:21:46
    基于CAD模型的三坐標測量機測量點分布規(guī)劃
    PM2.5空中探測器的設計
    溝槽對湍流邊界層中展向渦影響的實驗研究
    開縫圓柱縫隙傾斜角對脫落渦的影響
    美女cb高潮喷水在线观看| 热99re8久久精品国产| 啦啦啦啦在线视频资源| 国产91av在线免费观看| 一本精品99久久精品77| 全区人妻精品视频| 亚洲精品乱码久久久久久按摩| av在线老鸭窝| 人人妻人人看人人澡| av在线蜜桃| 1000部很黄的大片| 国产国拍精品亚洲av在线观看| 一个人免费在线观看电影| 国产成人a∨麻豆精品| 女人被狂操c到高潮| 久久久久久久久久久丰满| 欧美激情久久久久久爽电影| 欧美bdsm另类| 卡戴珊不雅视频在线播放| 久久九九热精品免费| 精品人妻一区二区三区麻豆| 国产熟女欧美一区二区| 长腿黑丝高跟| 最近2019中文字幕mv第一页| 精品国内亚洲2022精品成人| 日韩av不卡免费在线播放| 波多野结衣高清作品| 亚洲欧美清纯卡通| 日韩视频在线欧美| 天天躁夜夜躁狠狠久久av| 国产av不卡久久| 国产精华一区二区三区| 免费人成在线观看视频色| 久久精品国产亚洲av涩爱 | 欧美一区二区精品小视频在线| 美女大奶头视频| 亚洲精品国产成人久久av| 久久午夜福利片| 99久久精品国产国产毛片| 国产精品精品国产色婷婷| 亚洲欧美日韩卡通动漫| 欧美性猛交黑人性爽| 97人妻精品一区二区三区麻豆| 精品一区二区三区人妻视频| 两个人的视频大全免费| 神马国产精品三级电影在线观看| 波多野结衣高清无吗| 3wmmmm亚洲av在线观看| 国产精品免费一区二区三区在线| 欧美精品国产亚洲| 国产精品精品国产色婷婷| 好男人视频免费观看在线| 久99久视频精品免费| 免费看日本二区| 乱码一卡2卡4卡精品| 欧美日韩国产亚洲二区| 国产午夜福利久久久久久| 国产真实伦视频高清在线观看| 亚洲国产欧美人成| 最近2019中文字幕mv第一页| 18+在线观看网站| 亚洲精品久久久久久婷婷小说 | 久久久久久国产a免费观看| 能在线免费观看的黄片| 熟妇人妻久久中文字幕3abv| 婷婷六月久久综合丁香| 午夜老司机福利剧场| 国产在视频线在精品| 精品不卡国产一区二区三区| 日韩欧美三级三区| 欧美xxxx性猛交bbbb| 免费一级毛片在线播放高清视频| 国产三级中文精品| 一夜夜www| 成熟少妇高潮喷水视频| 在线观看66精品国产| 九九久久精品国产亚洲av麻豆| 日本爱情动作片www.在线观看| 久久人人爽人人片av| 日韩在线高清观看一区二区三区| 欧美日本亚洲视频在线播放| 亚洲欧美日韩卡通动漫| 亚洲av成人精品一区久久| 日本黄大片高清| 国产亚洲91精品色在线| 桃色一区二区三区在线观看| 黄色日韩在线| 哪里可以看免费的av片| 又爽又黄无遮挡网站| 伊人久久精品亚洲午夜| 国内揄拍国产精品人妻在线| 国内精品宾馆在线| 久久久久免费精品人妻一区二区| 久久99蜜桃精品久久| 我要搜黄色片| 又爽又黄a免费视频| 国产熟女欧美一区二区| 中文字幕熟女人妻在线| 中出人妻视频一区二区| 在线天堂最新版资源| 亚洲自偷自拍三级| 桃色一区二区三区在线观看| 色综合色国产| АⅤ资源中文在线天堂| 亚洲成a人片在线一区二区| 国产精品99久久久久久久久| 国产一区二区在线观看日韩| 毛片女人毛片| 成人永久免费在线观看视频| 美女国产视频在线观看| 国产精品美女特级片免费视频播放器| av在线亚洲专区| 免费观看a级毛片全部| 亚洲国产精品成人综合色| 又爽又黄a免费视频| 亚洲一区二区三区色噜噜| 国产一区二区激情短视频| 人妻系列 视频| 两性午夜刺激爽爽歪歪视频在线观看| 午夜a级毛片| 欧美xxxx黑人xx丫x性爽| 男女那种视频在线观看| 久久人人精品亚洲av| 亚洲av熟女| 亚洲国产精品国产精品| 欧美+亚洲+日韩+国产| 我要看日韩黄色一级片| av黄色大香蕉| 神马国产精品三级电影在线观看| 少妇猛男粗大的猛烈进出视频 | 国产成人福利小说| 日本黄色片子视频| 亚洲国产高清在线一区二区三| 日日啪夜夜撸| 免费大片18禁| 嘟嘟电影网在线观看| 精品国内亚洲2022精品成人| av国产免费在线观看| 一区二区三区免费毛片| 午夜视频国产福利| 亚洲美女视频黄频| 亚洲人与动物交配视频| 三级男女做爰猛烈吃奶摸视频| www.色视频.com| 久久久久网色| 黄色配什么色好看| 麻豆国产97在线/欧美| 亚洲成a人片在线一区二区| 少妇高潮的动态图| 久久精品综合一区二区三区| 亚洲国产精品久久男人天堂| 欧美精品一区二区大全| 日本免费a在线| 日本在线视频免费播放| 国产乱人视频| 欧美成人a在线观看| 日韩视频在线欧美| 自拍偷自拍亚洲精品老妇| 日韩av在线大香蕉| 少妇的逼水好多| 3wmmmm亚洲av在线观看| 日韩 亚洲 欧美在线| 久久久久久久久久成人| 亚洲国产欧洲综合997久久,| 精品日产1卡2卡| av.在线天堂| a级毛片免费高清观看在线播放| 又爽又黄无遮挡网站| 波野结衣二区三区在线| 舔av片在线| 内地一区二区视频在线| av在线蜜桃| 国产成人福利小说| 日本一二三区视频观看| 免费看日本二区| 亚洲性久久影院| 成人特级黄色片久久久久久久| 国产精品一区二区三区四区免费观看| 成年女人看的毛片在线观看| 黄片wwwwww| 观看美女的网站| 人妻系列 视频| 免费av观看视频| 在线观看av片永久免费下载| 国产精品一及| 一区二区三区免费毛片| 亚洲不卡免费看| 婷婷亚洲欧美| av国产免费在线观看| 免费大片18禁| 性欧美人与动物交配| 国产精品不卡视频一区二区| 蜜臀久久99精品久久宅男| 亚洲最大成人av| 欧美色欧美亚洲另类二区| 黄色日韩在线| 一级毛片久久久久久久久女| 大又大粗又爽又黄少妇毛片口| 午夜福利在线在线| 高清毛片免费观看视频网站| 搞女人的毛片| 国产一区二区激情短视频| 欧美极品一区二区三区四区| 亚洲精品亚洲一区二区| 国内精品宾馆在线| 国产精品一区二区在线观看99 | 99久久久亚洲精品蜜臀av| 在线观看美女被高潮喷水网站| 免费看日本二区| 国内精品美女久久久久久| 深爱激情五月婷婷| 亚洲在线自拍视频| 九九久久精品国产亚洲av麻豆| 国产片特级美女逼逼视频| 日韩国内少妇激情av| 深夜精品福利| 91精品国产九色| 亚洲欧美精品自产自拍| 亚洲国产精品成人综合色| 天堂av国产一区二区熟女人妻| 精华霜和精华液先用哪个| 亚洲av免费高清在线观看| 亚洲性久久影院| av又黄又爽大尺度在线免费看 | 日本在线视频免费播放| 只有这里有精品99| 国产淫片久久久久久久久| 大香蕉久久网| 国产视频首页在线观看| 精品久久久久久久久久久久久| ponron亚洲| 欧美zozozo另类| 一级毛片aaaaaa免费看小| 国产不卡一卡二| 国产精品美女特级片免费视频播放器| 成人亚洲欧美一区二区av| 乱人视频在线观看| 国产激情偷乱视频一区二区| 黑人高潮一二区| 国产单亲对白刺激| 天天躁夜夜躁狠狠久久av| 97热精品久久久久久| 国产成人午夜福利电影在线观看| 国产精品,欧美在线| 夜夜夜夜夜久久久久| 国产精品久久久久久久电影| 卡戴珊不雅视频在线播放| 日韩欧美在线乱码| 亚洲av成人av| 丝袜喷水一区| 婷婷六月久久综合丁香| 人妻久久中文字幕网| 午夜福利高清视频| 在现免费观看毛片| 亚洲人成网站在线播放欧美日韩| 免费看美女性在线毛片视频| 亚洲人成网站高清观看| 国产一区二区在线观看日韩| 免费看av在线观看网站| 日韩av不卡免费在线播放| 亚洲va在线va天堂va国产| 成人永久免费在线观看视频| 久久久精品欧美日韩精品| 午夜精品一区二区三区免费看| 美女大奶头视频| 成人性生交大片免费视频hd| 伦理电影大哥的女人| 日韩三级伦理在线观看| 搡老妇女老女人老熟妇| 久久这里只有精品中国| 性插视频无遮挡在线免费观看| 99热全是精品| 国产精品一区二区三区四区免费观看| 高清日韩中文字幕在线| 欧美xxxx性猛交bbbb| 亚洲最大成人中文| 久久亚洲精品不卡| 简卡轻食公司| 九草在线视频观看| 美女黄网站色视频| 国产精品嫩草影院av在线观看| 中文在线观看免费www的网站| 国产毛片a区久久久久| 草草在线视频免费看| 99riav亚洲国产免费| 国产精品乱码一区二三区的特点| 久久精品综合一区二区三区| 欧美潮喷喷水| 国产成人精品婷婷| 国产成人精品一,二区 | 婷婷色av中文字幕| av在线观看视频网站免费| 国产精品久久久久久精品电影小说 | 啦啦啦韩国在线观看视频| 国产高清激情床上av| 婷婷色av中文字幕| 男人舔奶头视频| 在现免费观看毛片| 精品一区二区三区人妻视频| 午夜精品在线福利| 不卡一级毛片| 国产精华一区二区三区| 精品熟女少妇av免费看| 男人的好看免费观看在线视频| 欧美+日韩+精品| 久久久久久国产a免费观看| 亚洲av成人精品一区久久| 国产精品久久久久久精品电影| 国产69精品久久久久777片| 欧美一区二区精品小视频在线| 精品久久久久久久末码| 一个人看的www免费观看视频| 亚洲无线在线观看| 欧美一级a爱片免费观看看| 久久精品国产亚洲av香蕉五月| 国产视频内射| 中文欧美无线码| 97超视频在线观看视频| 久久精品国产鲁丝片午夜精品| 村上凉子中文字幕在线| 欧美3d第一页| 亚洲自拍偷在线| 秋霞在线观看毛片| 亚洲av中文av极速乱| 精品一区二区三区视频在线| 国产女主播在线喷水免费视频网站 | 国产一区二区在线观看日韩| 成人亚洲精品av一区二区| 国产av在哪里看| 少妇人妻精品综合一区二区 | 毛片一级片免费看久久久久| 国产精品一区二区三区四区免费观看| 日韩成人av中文字幕在线观看| 国产av麻豆久久久久久久| 亚洲国产欧美在线一区| 成人一区二区视频在线观看| 中文精品一卡2卡3卡4更新| 十八禁国产超污无遮挡网站| 国产成年人精品一区二区| 久久久久久久久久久丰满| 成人国产麻豆网| 99久国产av精品| 久久99蜜桃精品久久| 日本撒尿小便嘘嘘汇集6| 国内少妇人妻偷人精品xxx网站| 亚洲精品国产av成人精品| 精品不卡国产一区二区三区| 18禁在线无遮挡免费观看视频| 欧美+日韩+精品| 偷拍熟女少妇极品色| 91精品一卡2卡3卡4卡| 综合色av麻豆| 国产精品不卡视频一区二区| 免费观看人在逋| 伦精品一区二区三区| 有码 亚洲区| 高清毛片免费看| 在现免费观看毛片| 成人综合一区亚洲| 五月伊人婷婷丁香| 久久午夜亚洲精品久久| 日韩欧美在线乱码| 国产免费一级a男人的天堂| 天堂√8在线中文| 99热这里只有是精品在线观看| a级一级毛片免费在线观看| 日韩中字成人| 桃色一区二区三区在线观看| 国产一区二区三区av在线 | 又粗又爽又猛毛片免费看| 成人av在线播放网站| 久久久久久九九精品二区国产| 毛片女人毛片| 欧美一区二区精品小视频在线| 国产精华一区二区三区| 国产免费一级a男人的天堂| 好男人在线观看高清免费视频| 成人亚洲欧美一区二区av| 能在线免费看毛片的网站| 三级男女做爰猛烈吃奶摸视频| 日日摸夜夜添夜夜爱| 成人高潮视频无遮挡免费网站| 99九九线精品视频在线观看视频| 国产精品人妻久久久久久| 最近视频中文字幕2019在线8| 国产av麻豆久久久久久久| 久久久久久久午夜电影| kizo精华| 免费观看人在逋| 少妇熟女aⅴ在线视频| 亚洲欧美日韩卡通动漫| 亚洲成人av在线免费| 高清毛片免费观看视频网站| 亚洲成人久久性| 简卡轻食公司| 国产三级在线视频| 欧美极品一区二区三区四区| 久久精品国产清高在天天线| 青春草视频在线免费观看| 少妇的逼好多水| 麻豆一二三区av精品| 丰满乱子伦码专区| 国产午夜精品论理片| 综合色丁香网| 老熟妇乱子伦视频在线观看| 成人无遮挡网站| 亚洲国产高清在线一区二区三| 高清毛片免费观看视频网站| 99久久久亚洲精品蜜臀av| 最近的中文字幕免费完整| 此物有八面人人有两片| 99热这里只有是精品在线观看| 国产亚洲5aaaaa淫片| 夜夜夜夜夜久久久久| 精华霜和精华液先用哪个| 亚洲不卡免费看| 丰满人妻一区二区三区视频av| 99久久成人亚洲精品观看| 日韩欧美三级三区| 看片在线看免费视频| 激情 狠狠 欧美| 成人永久免费在线观看视频| 六月丁香七月| 人妻夜夜爽99麻豆av| 午夜激情欧美在线| www日本黄色视频网| 高清毛片免费看| 我的女老师完整版在线观看| 国产精品野战在线观看| 欧美性感艳星| 老熟妇乱子伦视频在线观看| 亚洲四区av| ponron亚洲| 国产视频内射| 国产伦精品一区二区三区视频9| 国产久久久一区二区三区| 亚洲一级一片aⅴ在线观看| 国产亚洲av片在线观看秒播厂 | 日本-黄色视频高清免费观看| 一本久久精品| 国产三级中文精品| 亚洲av电影不卡..在线观看| 亚洲一级一片aⅴ在线观看| 亚洲丝袜综合中文字幕| 成人特级av手机在线观看| 久久久久免费精品人妻一区二区| а√天堂www在线а√下载| 99久久九九国产精品国产免费| 性插视频无遮挡在线免费观看| 三级男女做爰猛烈吃奶摸视频| 午夜老司机福利剧场| 悠悠久久av| 国产日本99.免费观看| 欧美bdsm另类| 精品久久久噜噜| 久久久精品欧美日韩精品| 久久鲁丝午夜福利片| 国产成人影院久久av| 悠悠久久av| 夜夜夜夜夜久久久久| 亚洲欧美成人综合另类久久久 | 日本-黄色视频高清免费观看| 亚洲欧美日韩高清专用| 国产av在哪里看| 午夜福利成人在线免费观看| 亚洲一级一片aⅴ在线观看| 联通29元200g的流量卡| 国产一区亚洲一区在线观看| 国产亚洲欧美98| 婷婷色av中文字幕| 亚洲中文字幕一区二区三区有码在线看| 欧美极品一区二区三区四区| 日本-黄色视频高清免费观看| 一区二区三区四区激情视频 | 男女下面进入的视频免费午夜| 最近的中文字幕免费完整| 亚洲成人久久性| 老司机影院成人| 日本-黄色视频高清免费观看| 国产成人a∨麻豆精品| 色综合色国产| 亚洲欧美日韩卡通动漫| 两性午夜刺激爽爽歪歪视频在线观看| 色哟哟·www| 青春草视频在线免费观看| 亚洲国产欧美在线一区| 亚洲av成人精品一区久久| 亚洲国产欧洲综合997久久,| 国产成人一区二区在线| 18禁裸乳无遮挡免费网站照片| 寂寞人妻少妇视频99o| 精品久久国产蜜桃| 亚洲,欧美,日韩| 性色avwww在线观看| 22中文网久久字幕| 亚洲av电影不卡..在线观看| 精品午夜福利在线看| 日韩成人伦理影院| 国产精品久久视频播放| 亚洲精品乱码久久久久久按摩| 国产单亲对白刺激| 在线免费观看的www视频| 少妇的逼水好多| 91久久精品国产一区二区成人| 插逼视频在线观看| 国产日本99.免费观看| 日韩制服骚丝袜av| 天堂av国产一区二区熟女人妻| 国产精品女同一区二区软件| 美女脱内裤让男人舔精品视频 | 在线观看66精品国产| 能在线免费观看的黄片| avwww免费| 可以在线观看毛片的网站| 国产高清有码在线观看视频| 国产午夜福利久久久久久| 亚洲一区高清亚洲精品| 亚洲av免费在线观看| 欧美高清成人免费视频www| 99久久无色码亚洲精品果冻| 日韩人妻高清精品专区| 中文字幕免费在线视频6| 天堂√8在线中文| 日韩,欧美,国产一区二区三区 | 国产在线精品亚洲第一网站| 欧美性猛交╳xxx乱大交人| 女同久久另类99精品国产91| 亚洲欧美中文字幕日韩二区| 少妇丰满av| 国产av麻豆久久久久久久| 69av精品久久久久久| 国产精品综合久久久久久久免费| 男人狂女人下面高潮的视频| 国产精品一区二区三区四区久久| 18+在线观看网站| 成人无遮挡网站| 国产亚洲5aaaaa淫片| 中文在线观看免费www的网站| 亚洲国产欧美人成| .国产精品久久| 久久精品国产亚洲网站| 99久久中文字幕三级久久日本| 欧美变态另类bdsm刘玥| 性插视频无遮挡在线免费观看| 悠悠久久av| 精品久久久久久成人av| 国产精品日韩av在线免费观看| 亚洲国产日韩欧美精品在线观看| 亚洲国产精品成人综合色| 男的添女的下面高潮视频| 高清毛片免费看| 狂野欧美激情性xxxx在线观看| 亚洲色图av天堂| 人妻少妇偷人精品九色| 一区福利在线观看| 国产精品人妻久久久久久| 国产成人影院久久av| 亚洲成人精品中文字幕电影| 黄色欧美视频在线观看| 69av精品久久久久久| 欧美+亚洲+日韩+国产| 国产在视频线在精品| 免费看美女性在线毛片视频| 国内精品美女久久久久久| 成人亚洲精品av一区二区| 一级毛片电影观看 | 悠悠久久av| 蜜桃久久精品国产亚洲av| 狂野欧美激情性xxxx在线观看| 三级经典国产精品| 久久久久性生活片| 深夜精品福利| 99久久九九国产精品国产免费| 久久久国产成人免费| 26uuu在线亚洲综合色| 久久这里有精品视频免费| 一区二区三区免费毛片| 亚洲成人精品中文字幕电影| 日韩视频在线欧美| av卡一久久| 蜜臀久久99精品久久宅男| 精品久久久久久久久久免费视频| av在线天堂中文字幕| 国产精品久久久久久久电影| 如何舔出高潮| 91久久精品国产一区二区三区| 国产白丝娇喘喷水9色精品| 美女 人体艺术 gogo| 91麻豆精品激情在线观看国产| 国产亚洲精品av在线| 99久国产av精品国产电影| www日本黄色视频网| 国产午夜精品一二区理论片| 99久久精品国产国产毛片| 少妇丰满av| 国产亚洲精品av在线| 好男人在线观看高清免费视频| 欧美日韩精品成人综合77777| 亚洲三级黄色毛片| 国产精品久久久久久av不卡| 最好的美女福利视频网| 黄色欧美视频在线观看| 精品久久久久久久末码| 少妇人妻一区二区三区视频| 国产欧美日韩精品一区二区| 一区二区三区免费毛片| 成人三级黄色视频| 小蜜桃在线观看免费完整版高清| 亚洲五月天丁香| 久久99精品国语久久久| 欧美最黄视频在线播放免费| 美女内射精品一级片tv| 亚洲乱码一区二区免费版| 啦啦啦韩国在线观看视频| 一级毛片我不卡|