• <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空中探測器的設計
    溝槽對湍流邊界層中展向渦影響的實驗研究
    開縫圓柱縫隙傾斜角對脫落渦的影響
    久久久久国内视频| 成年人午夜在线观看视频| 天堂俺去俺来也www色官网| 一级a爱片免费观看的视频| 女性生殖器流出的白浆| 在线观看午夜福利视频| 中文字幕色久视频| 欧美日韩国产mv在线观看视频| 动漫黄色视频在线观看| 一a级毛片在线观看| 女人久久www免费人成看片| 色精品久久人妻99蜜桃| 欧美中文综合在线视频| 又紧又爽又黄一区二区| 亚洲国产看品久久| 国产精品久久久久久人妻精品电影| 精品一区二区三区四区五区乱码| 亚洲中文日韩欧美视频| 好看av亚洲va欧美ⅴa在| 亚洲一区中文字幕在线| 精品亚洲成a人片在线观看| 99国产极品粉嫩在线观看| 99在线人妻在线中文字幕 | 午夜91福利影院| 12—13女人毛片做爰片一| 在线观看免费视频网站a站| 日韩熟女老妇一区二区性免费视频| 后天国语完整版免费观看| av一本久久久久| 成年人午夜在线观看视频| 精品久久久久久,| 人人妻人人澡人人爽人人夜夜| 亚洲色图 男人天堂 中文字幕| 国产片内射在线| 国产精品一区二区精品视频观看| 精品国产美女av久久久久小说| 黄色视频,在线免费观看| 精品熟女少妇八av免费久了| 亚洲午夜理论影院| 无遮挡黄片免费观看| 黄色视频,在线免费观看| 国产成人av教育| 精品国产乱子伦一区二区三区| 一边摸一边做爽爽视频免费| 亚洲成人国产一区在线观看| 欧美成狂野欧美在线观看| 久久影院123| 欧美日韩av久久| 搡老岳熟女国产| 国产有黄有色有爽视频| 少妇的丰满在线观看| 香蕉久久夜色| 99re在线观看精品视频| 国产精品 国内视频| 欧美一级毛片孕妇| 下体分泌物呈黄色| 精品人妻在线不人妻| 男女免费视频国产| 色综合欧美亚洲国产小说| 精品久久久精品久久久| 亚洲第一av免费看| av网站免费在线观看视频| 婷婷丁香在线五月| 欧美日韩av久久| av在线播放免费不卡| 亚洲av第一区精品v没综合| 欧美丝袜亚洲另类 | 精品第一国产精品| 欧美日韩亚洲高清精品| 人人妻人人澡人人爽人人夜夜| 18禁黄网站禁片午夜丰满| 啦啦啦视频在线资源免费观看| 午夜精品国产一区二区电影| 美女高潮到喷水免费观看| 69精品国产乱码久久久| xxxhd国产人妻xxx| 韩国精品一区二区三区| 国产亚洲精品久久久久5区| 夜夜夜夜夜久久久久| 黑人操中国人逼视频| 老熟妇仑乱视频hdxx| 国产一区二区三区综合在线观看| 国产成人av教育| 黄色a级毛片大全视频| 精品一区二区三区视频在线观看免费 | 不卡一级毛片| 视频区欧美日本亚洲| 无人区码免费观看不卡| 搡老熟女国产l中国老女人| 日本五十路高清| 日本黄色日本黄色录像| 亚洲人成伊人成综合网2020| 亚洲va日本ⅴa欧美va伊人久久| 国产又色又爽无遮挡免费看| 亚洲人成电影观看| 人人澡人人妻人| 国产99白浆流出| 他把我摸到了高潮在线观看| 欧美成狂野欧美在线观看| 午夜免费鲁丝| 午夜免费成人在线视频| 久久中文字幕人妻熟女| 法律面前人人平等表现在哪些方面| 欧美激情久久久久久爽电影 | 女人高潮潮喷娇喘18禁视频| 免费高清在线观看日韩| 18禁裸乳无遮挡动漫免费视频| 亚洲自偷自拍图片 自拍| 亚洲国产精品sss在线观看 | 中文字幕人妻熟女乱码| 女人爽到高潮嗷嗷叫在线视频| 久久久国产欧美日韩av| 亚洲一码二码三码区别大吗| 亚洲欧洲精品一区二区精品久久久| 国产精品1区2区在线观看. | 亚洲欧洲精品一区二区精品久久久| 亚洲成人免费电影在线观看| av有码第一页| 18禁裸乳无遮挡动漫免费视频| 欧美国产精品va在线观看不卡| 99re在线观看精品视频| 青草久久国产| 丝袜美足系列| 免费黄频网站在线观看国产| 国产成人影院久久av| 黄片播放在线免费| 女人被躁到高潮嗷嗷叫费观| 国产精品av久久久久免费| 女人被狂操c到高潮| 亚洲va日本ⅴa欧美va伊人久久| 国产又色又爽无遮挡免费看| 亚洲熟女精品中文字幕| 精品久久久精品久久久| 一进一出好大好爽视频| 无人区码免费观看不卡| 老司机在亚洲福利影院| 欧美日本中文国产一区发布| 亚洲一区高清亚洲精品| 在线观看66精品国产| 宅男免费午夜| 久久这里只有精品19| 欧美乱码精品一区二区三区| 久久狼人影院| 午夜精品久久久久久毛片777| 王馨瑶露胸无遮挡在线观看| 人妻 亚洲 视频| 亚洲人成电影观看| 黑人操中国人逼视频| 国产又爽黄色视频| 国产欧美日韩一区二区精品| 中文字幕精品免费在线观看视频| 久久精品国产清高在天天线| aaaaa片日本免费| 免费在线观看视频国产中文字幕亚洲| 日韩欧美三级三区| 国产亚洲欧美98| 亚洲国产精品合色在线| 中文字幕人妻丝袜一区二区| bbb黄色大片| 精品久久久精品久久久| 亚洲第一青青草原| 日韩有码中文字幕| 一本大道久久a久久精品| 国产不卡一卡二| 最新的欧美精品一区二区| av国产精品久久久久影院| 午夜免费成人在线视频| 超色免费av| 亚洲色图综合在线观看| videos熟女内射| 欧美在线一区亚洲| 亚洲色图综合在线观看| av欧美777| 成年人黄色毛片网站| 丝袜美足系列| 1024香蕉在线观看| 国产精品亚洲av一区麻豆| 曰老女人黄片| 露出奶头的视频| 亚洲avbb在线观看| 久久久久久久国产电影| 久久亚洲真实| 午夜日韩欧美国产| 一进一出好大好爽视频| 女人被狂操c到高潮| 少妇 在线观看| 日本黄色日本黄色录像| 成人国语在线视频| 妹子高潮喷水视频| 69av精品久久久久久| 黄片播放在线免费| av在线播放免费不卡| 如日韩欧美国产精品一区二区三区| 在线观看免费高清a一片| 丁香欧美五月| 日本一区二区免费在线视频| 亚洲中文字幕日韩| 少妇裸体淫交视频免费看高清 | 日日爽夜夜爽网站| 丁香欧美五月| 在线观看免费午夜福利视频| av免费在线观看网站| 亚洲av成人av| 一级片'在线观看视频| 国产99白浆流出| 精品国产乱码久久久久久男人| 正在播放国产对白刺激| 亚洲国产毛片av蜜桃av| 日韩有码中文字幕| 一边摸一边做爽爽视频免费| 国内毛片毛片毛片毛片毛片| 亚洲在线自拍视频| 老熟女久久久| 夜夜爽天天搞| 日日摸夜夜添夜夜添小说| 欧美av亚洲av综合av国产av| 69精品国产乱码久久久| 变态另类成人亚洲欧美熟女 | 亚洲精品一二三| 美女 人体艺术 gogo| 久久久久国产精品人妻aⅴ院 | 国产免费男女视频| 亚洲av成人av| 国产精品亚洲av一区麻豆| 国产亚洲精品久久久久久毛片 | 国产精品99久久99久久久不卡| 啦啦啦 在线观看视频| 久久草成人影院| 日韩中文字幕欧美一区二区| 51午夜福利影视在线观看| 精品久久久久久久久久免费视频 | 亚洲精品久久成人aⅴ小说| 欧美日韩福利视频一区二区| 黑人巨大精品欧美一区二区mp4| 国产三级黄色录像| 久久ye,这里只有精品| 国产深夜福利视频在线观看| 国产av又大| 757午夜福利合集在线观看| 十八禁网站免费在线| 自拍欧美九色日韩亚洲蝌蚪91| 妹子高潮喷水视频| а√天堂www在线а√下载 | 亚洲五月天丁香| 一级毛片精品| 另类亚洲欧美激情| 亚洲avbb在线观看| 午夜福利视频在线观看免费| 亚洲色图 男人天堂 中文字幕| 亚洲va日本ⅴa欧美va伊人久久| 天堂动漫精品| 国产精品欧美亚洲77777| 自线自在国产av| 在线观看午夜福利视频| 亚洲国产精品合色在线| 俄罗斯特黄特色一大片| 欧美久久黑人一区二区| 在线观看免费视频网站a站| 三级毛片av免费| 国产精品亚洲一级av第二区| 国产伦人伦偷精品视频| 精品高清国产在线一区| 精品乱码久久久久久99久播| 欧美日韩亚洲国产一区二区在线观看 | 亚洲av美国av| 多毛熟女@视频| 亚洲伊人色综图| 黄片大片在线免费观看| 久久精品亚洲熟妇少妇任你| 午夜福利在线免费观看网站| 精品国产乱子伦一区二区三区| 一a级毛片在线观看| 三级毛片av免费| 黄网站色视频无遮挡免费观看| 国产有黄有色有爽视频| 色播在线永久视频| 老汉色av国产亚洲站长工具| 成人18禁在线播放| 91成人精品电影| 欧美成狂野欧美在线观看| 曰老女人黄片| 亚洲七黄色美女视频| 美女福利国产在线| 亚洲av美国av| 亚洲欧洲精品一区二区精品久久久| 国产野战对白在线观看| 国产精品美女特级片免费视频播放器 | 免费av中文字幕在线| 免费看a级黄色片| 久久九九热精品免费| www日本在线高清视频| 国产一区二区激情短视频| 男女高潮啪啪啪动态图| 美女午夜性视频免费| 精品国产一区二区久久| 黄色 视频免费看| 母亲3免费完整高清在线观看| 天天躁日日躁夜夜躁夜夜| 久久久久精品人妻al黑| 亚洲在线自拍视频| 日韩欧美三级三区| 国产精品一区二区在线不卡| avwww免费| 亚洲成a人片在线一区二区| 身体一侧抽搐| 中文欧美无线码| 欧美日韩精品网址| 亚洲av电影在线进入| 亚洲成a人片在线一区二区| 一边摸一边做爽爽视频免费| 久久99一区二区三区| 法律面前人人平等表现在哪些方面| 国产精品综合久久久久久久免费 | 国产91精品成人一区二区三区| 老司机午夜福利在线观看视频| 18禁观看日本| 在线观看免费午夜福利视频| 12—13女人毛片做爰片一| 99riav亚洲国产免费| 亚洲一区中文字幕在线| 最近最新中文字幕大全免费视频| 一本大道久久a久久精品| 日韩熟女老妇一区二区性免费视频| 国产在线精品亚洲第一网站| 在线免费观看的www视频| 午夜久久久在线观看| 无遮挡黄片免费观看| 午夜老司机福利片| 老熟妇仑乱视频hdxx| 在线观看免费日韩欧美大片| 亚洲欧美日韩另类电影网站| 亚洲av欧美aⅴ国产| 欧美 亚洲 国产 日韩一| 成人免费观看视频高清| 国产99久久九九免费精品| 午夜免费观看网址| 伦理电影免费视频| 日韩欧美一区二区三区在线观看 | 好男人电影高清在线观看| 欧美成狂野欧美在线观看| 精品少妇一区二区三区视频日本电影| 午夜福利在线观看吧| 久久久久国产精品人妻aⅴ院 | 亚洲熟女毛片儿| а√天堂www在线а√下载 | 香蕉国产在线看| 欧美成人午夜精品| 桃红色精品国产亚洲av| 亚洲精品成人av观看孕妇| 亚洲aⅴ乱码一区二区在线播放 | 在线观看免费视频网站a站| 精品久久蜜臀av无| 一级a爱片免费观看的视频| 女人精品久久久久毛片| 午夜影院日韩av| 18禁黄网站禁片午夜丰满| 免费在线观看影片大全网站| 制服诱惑二区| 黄色 视频免费看| 中文欧美无线码| 中文亚洲av片在线观看爽 | 亚洲专区中文字幕在线| 在线观看日韩欧美| 少妇粗大呻吟视频| 免费观看a级毛片全部| 国产99白浆流出| 十八禁人妻一区二区| 老司机在亚洲福利影院| 电影成人av| 亚洲欧美精品综合一区二区三区| 母亲3免费完整高清在线观看| 91字幕亚洲| 久久香蕉激情| 三上悠亚av全集在线观看| 欧美成狂野欧美在线观看| 亚洲中文字幕日韩| 如日韩欧美国产精品一区二区三区| 一区二区三区激情视频| 国产av又大| 国产精品久久久久久精品古装| 嫁个100分男人电影在线观看| av免费在线观看网站| 母亲3免费完整高清在线观看| 美女国产高潮福利片在线看| 国产精品 国内视频| 777久久人妻少妇嫩草av网站| 亚洲黑人精品在线| 咕卡用的链子| 亚洲av第一区精品v没综合| av不卡在线播放| 国内毛片毛片毛片毛片毛片| 99国产精品一区二区三区| 欧美激情高清一区二区三区| 亚洲自偷自拍图片 自拍| 黄频高清免费视频| 欧美不卡视频在线免费观看 | 国内久久婷婷六月综合欲色啪| 一区二区日韩欧美中文字幕| 99久久综合精品五月天人人| 精品一区二区三卡| 亚洲国产中文字幕在线视频| 在线永久观看黄色视频| 国产aⅴ精品一区二区三区波| 国产蜜桃级精品一区二区三区 | 久久久久久久午夜电影 | 亚洲av成人不卡在线观看播放网| 日韩一卡2卡3卡4卡2021年| 伊人久久大香线蕉亚洲五| 国产精品美女特级片免费视频播放器 | 精品午夜福利视频在线观看一区| 黄色女人牲交| 人人妻人人添人人爽欧美一区卜| 日本黄色日本黄色录像| 久久精品国产综合久久久| 免费观看人在逋| 欧美精品啪啪一区二区三区| 一边摸一边抽搐一进一小说 | 亚洲人成伊人成综合网2020| 天天影视国产精品| 正在播放国产对白刺激| 亚洲熟妇熟女久久| 亚洲国产欧美网| 少妇裸体淫交视频免费看高清 | 每晚都被弄得嗷嗷叫到高潮| 日韩一卡2卡3卡4卡2021年| 无遮挡黄片免费观看| 高清av免费在线| 久久久精品免费免费高清| 伦理电影免费视频| 成年人免费黄色播放视频| 国产成人精品久久二区二区免费| 51午夜福利影视在线观看| 777久久人妻少妇嫩草av网站| 欧美亚洲 丝袜 人妻 在线| 亚洲精品国产区一区二| 很黄的视频免费| 波多野结衣av一区二区av| 美女扒开内裤让男人捅视频| 18禁国产床啪视频网站| 国产成人免费无遮挡视频| 午夜福利在线免费观看网站| 久久久久国产精品人妻aⅴ院 | 丝瓜视频免费看黄片| 久久久久久亚洲精品国产蜜桃av| 一区福利在线观看| xxx96com| 国产一区二区三区在线臀色熟女 | 天天影视国产精品| 搡老乐熟女国产| 啦啦啦 在线观看视频| 窝窝影院91人妻| 岛国毛片在线播放| 老司机深夜福利视频在线观看| 夜夜爽天天搞| 美女视频免费永久观看网站| 夜夜夜夜夜久久久久| avwww免费| 久久久久久久久免费视频了| 黄色视频不卡| 不卡av一区二区三区| 视频在线观看一区二区三区| 纯流量卡能插随身wifi吗| e午夜精品久久久久久久| 久久精品熟女亚洲av麻豆精品| 高清在线国产一区| 国产一区二区三区综合在线观看| 黑人操中国人逼视频| 欧美在线一区亚洲| 美女视频免费永久观看网站| aaaaa片日本免费| 精品视频人人做人人爽| 精品国产乱子伦一区二区三区| 天天添夜夜摸| 又黄又粗又硬又大视频| 啪啪无遮挡十八禁网站| 国产精品国产高清国产av | 亚洲精品成人av观看孕妇| 大香蕉久久成人网| 国产亚洲精品久久久久久毛片 | a在线观看视频网站| 亚洲第一欧美日韩一区二区三区| 在线观看舔阴道视频| 欧美黑人欧美精品刺激| 亚洲欧美激情综合另类| 丝袜在线中文字幕| 亚洲精品中文字幕一二三四区| 日日爽夜夜爽网站| 精品一区二区三卡| 天天添夜夜摸| 亚洲免费av在线视频| ponron亚洲| 热99久久久久精品小说推荐| 免费av中文字幕在线| 成年动漫av网址| 如日韩欧美国产精品一区二区三区| 97人妻天天添夜夜摸| 日韩欧美一区视频在线观看| 成熟少妇高潮喷水视频| 一进一出抽搐gif免费好疼 | 国产免费av片在线观看野外av| 男女之事视频高清在线观看| 91九色精品人成在线观看| 手机成人av网站| 制服诱惑二区| 老司机福利观看| 国产高清激情床上av| 国产精品电影一区二区三区 | 欧美人与性动交α欧美精品济南到| 黑人巨大精品欧美一区二区蜜桃| 99热国产这里只有精品6| 国产一区在线观看成人免费| 91大片在线观看| 午夜福利一区二区在线看| 黄色毛片三级朝国网站| 深夜精品福利| 一本综合久久免费| 免费看十八禁软件| av免费在线观看网站| 国产极品粉嫩免费观看在线| av在线播放免费不卡| 高清av免费在线| 黄色视频不卡| 欧美日韩精品网址| 国产精品久久视频播放| av线在线观看网站| 欧美中文综合在线视频| 国产av一区二区精品久久| 一级毛片高清免费大全| 在线看a的网站| 国产精品国产高清国产av | 国产深夜福利视频在线观看| 国产精品秋霞免费鲁丝片| 亚洲片人在线观看| 国产成人欧美在线观看 | 男人舔女人的私密视频| 亚洲全国av大片| 人妻久久中文字幕网| 亚洲成av片中文字幕在线观看| 一本大道久久a久久精品| 一夜夜www| 久久久久久免费高清国产稀缺| 亚洲精品乱久久久久久| 亚洲伊人色综图| av天堂久久9| 亚洲九九香蕉| 欧美精品亚洲一区二区| 91国产中文字幕| 精品亚洲成国产av| 欧美久久黑人一区二区| 91精品三级在线观看| 亚洲国产精品合色在线| 亚洲一卡2卡3卡4卡5卡精品中文| 午夜视频精品福利| 欧美另类亚洲清纯唯美| 18禁观看日本| 国产高清videossex| 亚洲av电影在线进入| 日日摸夜夜添夜夜添小说| 大陆偷拍与自拍| 国产精品亚洲一级av第二区| 欧美日韩一级在线毛片| 免费看a级黄色片| 久久久久久久久免费视频了| 日本a在线网址| 老司机在亚洲福利影院| 老司机午夜十八禁免费视频| 亚洲第一青青草原| 搡老乐熟女国产| 欧美日韩乱码在线| 国产精品久久久久久精品古装| 电影成人av| 亚洲欧美激情在线| 天堂中文最新版在线下载| 免费日韩欧美在线观看| 一进一出抽搐动态| 国产91精品成人一区二区三区| 国产精品一区二区在线不卡| 性色av乱码一区二区三区2| 欧美日韩亚洲高清精品| av一本久久久久| 日本精品一区二区三区蜜桃| 国产精品乱码一区二三区的特点 | 日韩三级视频一区二区三区| 成年版毛片免费区| 久久亚洲真实| 精品亚洲成a人片在线观看| 黄色怎么调成土黄色| 捣出白浆h1v1| 久热爱精品视频在线9| 在线av久久热| 久久精品国产清高在天天线| 亚洲精品av麻豆狂野| 在线观看午夜福利视频| 国产一区有黄有色的免费视频| www.精华液| 俄罗斯特黄特色一大片| 飞空精品影院首页| 欧美成狂野欧美在线观看| 在线观看免费视频网站a站| 韩国av一区二区三区四区| 欧美日韩视频精品一区| 国产深夜福利视频在线观看| 91老司机精品| av超薄肉色丝袜交足视频| 亚洲成人免费av在线播放| 亚洲综合色网址| 亚洲av日韩精品久久久久久密| 日韩中文字幕欧美一区二区| 国产精华一区二区三区| 国产精品香港三级国产av潘金莲| 91在线观看av| 中文字幕另类日韩欧美亚洲嫩草| 岛国在线观看网站| 脱女人内裤的视频| 露出奶头的视频|