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

    前混合水射流中丸粒速度的數(shù)值模擬與PIV 試驗研究*

    2023-09-07 07:54:28朱澤華劉蘭榮馬泳濤盧春生
    航空制造技術(shù) 2023年13期
    關(guān)鍵詞:水射流射流入口

    朱澤華,劉蘭榮,馬泳濤,盧春生,張 彬

    (1. 鄭州大學(xué)河南省資源與材料工業(yè)技術(shù)研究院,鄭州 450001;2. 鄭州大學(xué)機(jī)械與動力工程學(xué)院,鄭州 450001;3. School of Civil and Mechanical Engineering, Curtin University, Perth WA 6845, Australia)

    前混合水射流噴丸表層改性技術(shù)是一種新興的金屬工件表面處理技術(shù),其改性機(jī)理是將預(yù)混合的丸粒通過高壓水進(jìn)行加速,高速噴射到金屬工件表面,從而使工件產(chǎn)生塑性變形,形成殘余壓應(yīng)力層[1-2]。在前混合水射流中,射流束能量場的大小及分布直接影響改性效果,而射流束能量場與丸粒的速度密切相關(guān)[3-4],因此研究前混合水射流束中丸粒的速度分布,對揭示前混合水射流表層改性機(jī)理,提高表層改性能力及制定改性工藝參數(shù)具有重要意義。

    Basha 等[5]采用Fluent 軟件中的VOF 模型和DPM模型對后混合磨料水射流中噴嘴內(nèi)部氣-液-固三相的流動特性進(jìn)行了數(shù)值模擬,基于VOF 模型模擬了水-氣兩相流,基于DPM 模型模擬了顆粒軌跡。Long 等[6]采用了相同的耦合策略并研究了顆粒尺寸、密度、形狀因子對后混合射流噴嘴出口顆粒速度的影響。Liu 等[7]也采用VOF + DPM 的策略,并求解了超高速磨料水射流外流場的速度分布。Qiang 等[8]采用Fluent 軟件中歐拉模型和DPM 模型模擬了后混合射流系統(tǒng)中切割頭內(nèi)部的射流流動和顆?;旌线^程,分析了顆粒入口角度、入口位置、聚焦管收斂角對流體的加速過程的影響。

    在PIV測速方面,Thongkaew等[9]結(jié)合PIV和LIF (激光誘導(dǎo)熒光)技術(shù),對后混合磨料射流中的粒子軸向速度進(jìn)行試驗研究。Zeleňák 等[10]也采用相同技術(shù)對后混合磨料水射流的顆粒速度進(jìn)行試驗研究,并通過PTV算法分析單個粒子的運動。章文峰等[11]基于PIV 技術(shù)和圖像處理技術(shù)對前混合水射流中顆粒速度進(jìn)行試驗研究,同時得到了磨料的速度和位置信息,發(fā)現(xiàn)在14 mm 噴距范圍內(nèi),磨料速度先增大后減小,PIV 試驗采用錐直形噴嘴和球形陶瓷丸粒,噴嘴出口直徑為6 mm,丸粒密度為2.7 g/cm3。

    綜上所述,上述的研究多為單一的通過仿真模擬或者PIV 測速試驗研究液固混合射流中丸粒的速度,研究對象多為后混合水射流。后混合水射流是指在水射流形成之后,利用混合腔內(nèi)負(fù)壓引入磨料,最終形成后混合水射流,因此,后混合水射流中磨料與水混合不均勻,能量交換效率低。數(shù)值模擬時對于顆粒的建模多采用DPM 模型,并且只采用了單相耦合,未能考慮水丸之間的相互作用。本文通過耦合計算流體力學(xué)軟件Fluent和離散元軟件EDEM,對氣液固三相射流束進(jìn)行了數(shù)值模擬研究。結(jié)合PIV 測試技術(shù)和圖像處理技術(shù)對前混合水射流束的鑄鋼丸粒速度進(jìn)行了試驗測試,驗證了仿真模型的可信度。基于CFD-DEM 耦合方法探究噴嘴入口壓力、混合比 (丸粒體積分?jǐn)?shù))、噴嘴結(jié)構(gòu)參數(shù)對丸粒速度的影響,為制定工藝參數(shù)和優(yōu)化噴嘴結(jié)構(gòu)提供技術(shù)基礎(chǔ)。

    1 計算模型

    1.1 CFD 模型的構(gòu)建

    1.1.1 數(shù)學(xué)模型

    高壓射流通常處于湍流狀態(tài),所以系統(tǒng)除了求解連續(xù)性方程、N-S 方程等基本控制方程外,還需要遵循額外的湍流輸運方程[12-13]。采用雷諾平均法模擬湍流流動時,會使雷諾時均N-S 方程 (式(1))產(chǎn)生新的雷諾應(yīng)力未知項,即,為了封閉方程組,本文采用RealizableK-e湍流模型求解該項。

    式中,ρ為流體密度;μ為動力黏性系數(shù);p為壓力;ui和uj為流體速度;t為時間;x、z為位置;i和j的取值范圍為 1,2,3;Si為動量守恒方程的廣義源項。

    此外,前混合水射流包含氣液固3 相,屬于多相流動問題。在該多相流系統(tǒng)中,可以將氣液兩相看成是連續(xù)相,將丸??闯墒请x散相[5-6,14]。離散相的求解過程在1.2 節(jié)中介紹。關(guān)于連續(xù)相的求解,對氣液兩相利用連續(xù)介質(zhì)假設(shè)理論,在Euler 坐標(biāo)系下,采用有限體積法對其進(jìn)行模擬?;谠摾碚?,本文采用Fluent 中的VOF 模型對前混合水射流的連續(xù)相進(jìn)行模擬,需要注意的是,在VOF 模型中,連續(xù)相中的氣液兩相共享一組動量方程,但是考慮到離散相體積分?jǐn)?shù)的影響,需要對VOF 的連續(xù)性方程和N-S 方程進(jìn)行修改,修改后的連續(xù)性方程和N-S 方程[15]分別為

    式中,εf為網(wǎng)格孔隙率;Fs為表面張力;u為網(wǎng)格內(nèi)流體速度;ρf為網(wǎng)格內(nèi)流體密度;τ-為黏性應(yīng)力張量;g為重力加速度;fpf為流體-顆粒相互作用項的反作用力。

    1.1.2 物理模型及邊界條件

    選取水射流技術(shù)常用的錐直形噴嘴,計算區(qū)域包括噴嘴內(nèi)部和噴嘴外部直徑20 mm、長度60 mm 的圓柱區(qū)域。結(jié)構(gòu)化網(wǎng)格模型及邊界條件如圖1 所示,考慮到噴嘴壁面,在靠近噴嘴壁面處生成邊界層網(wǎng)格,并在靠近氣液兩相分界處加密網(wǎng)格,用于改善VOF 兩相界面追蹤效果。將Inlet 設(shè)置為壓力入口邊界條件,Wall 和Wall1 設(shè)置為無滑移速度壁面邊界條件,Wall2 和Outlet 設(shè)置為壓力出口邊界條件,壓力值設(shè)為 0.10 MPa(大氣壓)。

    圖1 網(wǎng)格和邊界條件示意圖Fig.1 Schematic diagram of grid and boundary conditions

    1.2 DEM 模型的構(gòu)建

    在前混合水射流多相流系統(tǒng)中,對于離散相的求解,可以在拉格朗日坐標(biāo)系下基于離散元 (DEM)方法進(jìn)行,該方法對所有顆粒進(jìn)行追蹤[16],考慮顆粒與顆粒、顆粒與邊界的相互作用?;谠摲椒?,本文采用EDEM對前混合水射流的離散相進(jìn)行模擬,并采用Hertz-Mindlin 無滑動接觸模型求解顆粒單元之間的接觸力。

    在EDEM 設(shè)置過程中,瑞利時間步取17%,總仿真時間為0.01 s。試驗中采用鑄鋼丸和硬質(zhì)合金噴嘴,材料參數(shù)如表1 所示。

    表1 材料參數(shù)設(shè)置Table 1 Parameter setting of materials

    1.3 相間作用力

    由前面的分析可知,對于前混合水射流多相流系統(tǒng),通過CFD-DEM 耦合方法在歐拉坐標(biāo)系中描述連續(xù)相的運動,在拉格朗日坐標(biāo)系下描述離散相的運動,通過相間作用力來實現(xiàn)連續(xù)相和離散相的雙向耦合。

    在前混合水射流系統(tǒng)中,所考慮的顆粒-流體相互作用力如式(4)所示[15]。

    式中,F(xiàn)pf是顆粒k受到的所有流體-顆粒之間的相互作用力;Fd是曳力,通過自由流曳力模型求解[17];是壓力梯度力;Fvm是虛擬質(zhì)量力;FB是Basset 力;FSaff是Saffman 升力;FMag是Magnus 升力。

    2 PIV 測試試驗

    2.1 PIV 試驗裝置介紹

    前混合水射流PIV 測速專用試驗臺如圖2 所示,包括兩部分,一是前混合水射流噴射系統(tǒng),用于提供前混合水射流;二是PIV 測速系統(tǒng) (美國TSI 公司)。PIV測速的工作過程為: (1)對相機(jī)進(jìn)行標(biāo)定,計算出校準(zhǔn)因子,得到圖像坐標(biāo)與真實坐標(biāo)之間的轉(zhuǎn)換關(guān)系; (2)噴射射流,相機(jī)與脈沖激光同步,在兩個連續(xù)激光脈沖期間捕獲射流圖像,完成圖像采集; (3)通過計算機(jī)系統(tǒng)對拍攝的一系列圖片進(jìn)行分析處理,進(jìn)而得到速度矢量及速度相關(guān)量。

    相機(jī)的分辨率為 2360×1776 像素,幀率為16 幀/s,最小跨幀時間為195 ns;激光器的脈沖頻率為15 Hz,脈沖能量為200 mJ。

    2.2 PIV 圖像處理

    PIV 試驗的拍攝范圍為從噴嘴出口到距離噴嘴出口約65 mm 的范圍,通過PIV 測速試驗采集的射流原始圖像的放大圖如圖3(a)和 (b)所示。原始圖像可能存在少量由于激光亮度過大導(dǎo)致的過度曝光點或丸粒邊界與射流之間的噪音,本文通過Matlab 圖像處理技術(shù)對原始圖像進(jìn)行預(yù)處理,從而削弱上述因素。最后通過基于PTV[18]算法的Matlab程序?qū)︻A(yù)處理后的A幀、B 幀圖像進(jìn)行分析處理,得到丸粒的速度及位置信息,速度矢量圖如圖3(c)所示。

    圖3 PIV 原始圖像及速度矢量圖Fig.3 PIV original images and velocity vector diagram

    經(jīng)過上述處理后,得到一對圖像中丸粒的速度相關(guān)信息。需要注意的是,丸粒的識別主要是通過給定的丸粒與背景之間的灰度值閾值和來判定的,但是在同一幅圖像中,每個丸粒與背景的灰度值差值有所不同,所以可能會誤識別丸粒,從而出現(xiàn)錯誤的速度值,因此需要對異常值進(jìn)行剔除,本文采用基于拉依達(dá)準(zhǔn)則的Matlab程序?qū)λ俣犬惓V颠M(jìn)行剔除[19]。

    2.3 試驗內(nèi)容

    通過對比三因素 (噴嘴入口壓力、混合比、噴嘴長度)、三水平下的仿真和PIV 測試結(jié)果,從多角度驗證仿真模型的準(zhǔn)確性??紤]到試驗安全性及測試條件,噴嘴入口壓力、混合比、噴嘴長度均限制在一定范圍內(nèi)。評價指標(biāo)包括兩個,一是外流場中丸粒的速度變化趨勢,用于定性地說明二者的符合程度;二是距離噴嘴出口10 mm、30 mm、50 mm 截面處丸粒的平均速度,分析仿真值與試驗值的誤差大小,用于定量說明仿真模型的準(zhǔn)確性。試驗時在磨料罐中添加800 mL 丸粒,每次拍攝100 組圖片。PIV 驗證試驗的試驗參數(shù)如表2 所示。

    表2 PIV 驗證試驗的參數(shù)設(shè)置Table 2 Parameter setting in PIV verification experiments

    2.4 試驗結(jié)果分析

    不同噴嘴入口壓力、混合比和噴嘴長度下,丸粒速度在外流場中的分布情況分別如圖4 ~ 6 所示。可以看出,在任一噴嘴入口壓力、混合比和噴嘴長度下,丸粒速度變化趨勢的仿真結(jié)果和試驗結(jié)果基本相同,在60 mm的靶距范圍內(nèi),兩種方法得到的丸粒速度基本保持不變。需要注意的是,本文的結(jié)論和文獻(xiàn)[11]結(jié)論不同,這可能是由于噴嘴出口直徑和丸粒密度不同引起的。

    圖4 不同噴嘴入口壓力下試驗與仿真的結(jié)果對比Fig.4 Comparison of experimental and simulation results under diffe ent nozzle inlet pressures

    圖5 不同混合比下試驗與仿真的結(jié)果對比Fig.5 Comparison of experimental and simulation results under diffe ent mixing ratios

    進(jìn)一步地,通過式 (5)對丸粒仿真速度的誤差率進(jìn)行分析,分別計算不同參數(shù)下距離噴嘴出口10 mm、30 mm、50 mm 截面處丸粒的仿真速度誤差率。通過分析可知,在4 MPa、6 MPa、8 MPa 3 組壓力下,丸粒仿真速度的誤差率絕對值最大分別為3.53%、2.10%、2.97%;在2.41%、4.47%、6.10% 3 組混合比下,丸粒仿真速度的誤差率絕對值最大分別為4.63%、4.84%、4.24%;在26 mm、31 mm、41 mm 3 組噴嘴長度下,丸粒仿真速度的誤差率絕對值最大分別為4.87%、5.20%、4.72%。3 組結(jié)果均表明在距離噴嘴出口距離為10 mm、30 mm、50 mm處,丸粒的仿真結(jié)果與試驗結(jié)果基本相同。

    式中,w為丸粒仿真速度的誤差率;vave,s為丸粒仿真平均速度;vave,e為丸粒試驗平均速度。在計算平均速度時,由于在某一時刻落在指定截面上的丸粒數(shù)量過少甚至為0,所以將該截面前后1 mm 范圍區(qū)域的丸粒均視為該截面上的丸粒。

    綜上所述,通過對不同噴嘴入口壓力、混合比和噴嘴長度下的仿真速度及試驗速度大小、誤差進(jìn)行對比分析,均驗證了仿真模型的正確性。

    3 前混合水射流中丸粒速度影響因素研究

    丸粒速度是影響前混合水射流改性效果的重要參數(shù),在進(jìn)行影響因素分析時,以經(jīng)過PIV 試驗驗證的仿真模型為基礎(chǔ),采用單因素分析法研究了噴嘴入口壓力、混合比、噴嘴結(jié)構(gòu)參數(shù)對前混合水射流中丸粒速度的影響規(guī)律,其中各參數(shù)的取值范圍以涵蓋實際前混合水射流噴丸表層改性技術(shù)為基礎(chǔ)。

    3.1 噴嘴入口壓力

    由于前混合水射流能量利用率高,工作壓力通常在40 MPa 以內(nèi),噴嘴入口壓力對丸粒最大速度的影響規(guī)律如圖7 所示,丸粒最大速度隨著噴嘴入口壓力的增大而增大。由液體孔口出流理論可知,噴嘴入口壓力增大使水的速度增大,高速水對丸粒起到加速作用,進(jìn)而使丸粒速度增大。壓力是影響丸粒速度的重要因素,通過調(diào)節(jié)壓力來調(diào)節(jié)丸粒速度是較為簡單的方式,因此對所測數(shù)據(jù)進(jìn)行擬合,得到了在上述約束條件下丸粒最大速度vmax,p與噴嘴入口壓力p的關(guān)系式為

    圖7 噴嘴入口壓力對丸粒最大速度的影響規(guī)律Fig.7 Effect of nozzle inlet p essure on the maximum velocity of shot

    3.2 混合比

    在改性工藝中,混合比取值通常在10%以內(nèi),混合比對丸粒最大速度的影響規(guī)律如圖8 所示,丸粒最大速度隨著混合比的增大而減小?;旌媳葟?%增大到10%,丸粒最大速度從226.96 m/s 降到201.46 m/s,減小了11.24%?;旌媳仍酱?,單位時間內(nèi)高速水需要加速的丸粒數(shù)量就越多,水的負(fù)載變大,這會使水的速度降低,進(jìn)而導(dǎo)致丸粒速度降低。在制定改性工藝時,若選取較小的混合比,可得到較高的丸粒速度,但是單位時間內(nèi)打擊至工件表面的丸粒數(shù)量較少,進(jìn)而影響改性效果。

    圖8 混合比對丸粒最大速度的影響規(guī)律Fig.8 Effect of mixing ratio on the maximum velocity of sho

    3.3 噴嘴長度

    噴嘴長度對丸粒最大速度的影響規(guī)律如圖9 所示,丸粒最大速度隨著圓柱段長度的增大而減小。圓柱段長度從18 mm 增大到66 mm,丸粒最大速度從229.83 m/s 降到211.75 m/s,減小了7.87%。圓柱段長度越長,水和噴嘴壁面的摩擦距離越長,并且丸粒與壁面的碰撞概率增大,這都會導(dǎo)致丸粒速度減小。

    圖9 噴嘴長度對丸粒最大速度的影響規(guī)律Fig.9 Effect of nozzle lengths on the maximum velocity of sho

    3.4 噴嘴出口直徑

    噴嘴出口直徑通常根據(jù)改性零部件的具體尺寸來確定,噴嘴出口直徑對丸粒最大速度的影響規(guī)律如圖10 所示,丸粒最大速度隨著噴嘴出口直徑的增大而增大。噴嘴出口直徑從1 mm 增大到2.5 mm,丸粒最大速度從219 m/s 增到240.33 m/s,增加了9.74%。噴嘴出口直徑越大,在相同混合比的情況下,丸粒與壁面碰撞的概率就越小,速度損失就越小。

    圖10 噴嘴出口直徑對丸粒最大速度的影響規(guī)律Fig.10 Effect of nozzle outlet diamete on the maximum velocity of shot

    3.5 噴嘴錐角

    在改性工藝中,為提高丸粒速度的同時減小噴嘴磨損,錐角不宜過大或過小,通常在30°左右,噴嘴錐角對丸粒最大速度的影響規(guī)律如圖11 所示,丸粒最大速度隨著錐角的增大而減小。錐角從10°增大到63°,丸粒最大速度從228.85 m/s 降到218.38 m/s,減小了4.58%。丸粒在噴嘴中主要是沿著噴嘴向前流動,錐角越大,丸粒與壁面碰撞的概率就越大,速度損失就越大。

    圖11 噴嘴錐角對丸粒最大速度的影響規(guī)律Fig.11 Effect of nozzle cone angle on the maximum velocity of sho

    4 結(jié)論

    (1)基于CFD-DEM 耦合方法模擬了前混合水射流噴嘴內(nèi)外流場的三相流動,在考慮顆粒碰撞作用的前提下,實現(xiàn)了顆粒與流體之間相互交換動量、能量的雙向耦合。

    (2)PIV 測速試驗對仿真模型的驗證結(jié)果表明,通過改變壓力、混合比和長徑比,在60 mm 靶距范圍內(nèi),仿真結(jié)果與試驗結(jié)果吻合度較好,驗證了該仿真模型的正確性。

    (3)前混合水射流中丸粒最大速度隨著噴嘴入口壓力、噴嘴出口直徑的增大而增大,隨著混合比、噴嘴長度、噴嘴錐角的增大而減小。改性效果不僅和丸粒速度相關(guān),還和丸粒性質(zhì)、分布相關(guān),應(yīng)綜合考慮多方因素,制定合理的工藝路線,得到最優(yōu)的改性效果。

    猜你喜歡
    水射流射流入口
    深海逃逸艙射流注水均壓過程仿真分析
    低壓天然氣泄漏射流擴(kuò)散特性研究
    煤氣與熱力(2022年4期)2022-05-23 12:45:00
    基于新一代稱重設(shè)備的入口治超勸返系統(tǒng)分析
    基于超高壓水射流的壓縮機(jī)智能拆解設(shè)備設(shè)計
    秘密入口
    高壓水射流切割混凝土試驗研究
    建筑科技(2018年6期)2018-08-30 03:41:04
    作品三
    第九道 靈化閣入口保衛(wèi)戰(zhàn)
    磨料水射流技術(shù)及其在水下結(jié)構(gòu)物切割中的應(yīng)用
    焊接(2015年5期)2015-07-18 11:03:40
    射流齒形噴嘴射流流場與氣動聲學(xué)分析
    av天堂中文字幕网| 久久久国产一区二区| 亚洲精品一区蜜桃| 老司机影院成人| 亚洲精品成人av观看孕妇| 午夜影院在线不卡| 在线天堂最新版资源| 精品国产乱码久久久久久小说| 色哟哟·www| 一级二级三级毛片免费看| 国产高清有码在线观看视频| 99热国产这里只有精品6| 人人妻人人澡人人看| 亚洲国产精品专区欧美| 亚洲精品国产色婷婷电影| 亚洲国产精品成人久久小说| 91精品国产九色| 亚洲欧美成人精品一区二区| 91午夜精品亚洲一区二区三区| 亚洲av成人精品一二三区| 少妇人妻 视频| 69精品国产乱码久久久| 纵有疾风起免费观看全集完整版| 全区人妻精品视频| 夫妻性生交免费视频一级片| 久久热精品热| 久久人人爽人人爽人人片va| 亚洲av免费高清在线观看| 纯流量卡能插随身wifi吗| 女的被弄到高潮叫床怎么办| 毛片一级片免费看久久久久| 亚洲精品日本国产第一区| 91成人精品电影| 亚洲av二区三区四区| 免费观看a级毛片全部| 精品99又大又爽又粗少妇毛片| 三上悠亚av全集在线观看 | av视频免费观看在线观看| 日韩精品免费视频一区二区三区 | 精品亚洲成a人片在线观看| 五月开心婷婷网| 少妇人妻一区二区三区视频| 黄色怎么调成土黄色| 国产视频内射| 又黄又爽又刺激的免费视频.| 一级毛片aaaaaa免费看小| 丁香六月天网| 久久这里有精品视频免费| 天天操日日干夜夜撸| 中文精品一卡2卡3卡4更新| 我要看黄色一级片免费的| 亚洲国产精品国产精品| 亚洲成人av在线免费| 日本vs欧美在线观看视频 | 免费观看av网站的网址| 高清黄色对白视频在线免费看 | 国产在视频线精品| 国产亚洲av片在线观看秒播厂| av福利片在线观看| 亚洲av男天堂| 一本大道久久a久久精品| 又爽又黄a免费视频| 夫妻性生交免费视频一级片| 我的女老师完整版在线观看| 亚洲高清免费不卡视频| 亚洲va在线va天堂va国产| 午夜老司机福利剧场| a级毛片免费高清观看在线播放| 久久久久网色| 久久这里有精品视频免费| 男女啪啪激烈高潮av片| 青春草亚洲视频在线观看| 国产视频内射| 国产精品99久久久久久久久| 69精品国产乱码久久久| 国产高清有码在线观看视频| 欧美日韩av久久| 欧美日韩av久久| 人妻少妇偷人精品九色| 免费在线观看成人毛片| 国产成人一区二区在线| 午夜av观看不卡| 久久97久久精品| 欧美成人精品欧美一级黄| 亚洲在久久综合| 婷婷色综合大香蕉| 欧美日韩视频高清一区二区三区二| 亚洲av免费高清在线观看| 亚洲欧美中文字幕日韩二区| 亚洲欧美一区二区三区国产| 亚洲av不卡在线观看| 岛国毛片在线播放| 亚洲av.av天堂| 久久国内精品自在自线图片| 精品一品国产午夜福利视频| 日日啪夜夜撸| 色视频在线一区二区三区| 日本黄色片子视频| 久久婷婷青草| 99re6热这里在线精品视频| 欧美激情国产日韩精品一区| 三级经典国产精品| 日韩人妻高清精品专区| 在线观看国产h片| 久久精品国产亚洲网站| 久久鲁丝午夜福利片| 精品国产乱码久久久久久小说| 亚洲高清免费不卡视频| 久久国产精品男人的天堂亚洲 | 国产精品欧美亚洲77777| 久久精品久久精品一区二区三区| 亚洲不卡免费看| 99精国产麻豆久久婷婷| 九色成人免费人妻av| 免费黄色在线免费观看| 欧美+日韩+精品| 少妇丰满av| 美女内射精品一级片tv| 一级爰片在线观看| 亚洲中文av在线| 女人久久www免费人成看片| 国产av码专区亚洲av| 国产精品人妻久久久影院| 熟女人妻精品中文字幕| 97超碰精品成人国产| 国产 一区精品| 国产免费视频播放在线视频| 亚洲国产毛片av蜜桃av| 亚洲精品乱码久久久久久按摩| 人人澡人人妻人| 久久99热6这里只有精品| 日日啪夜夜爽| 免费人成在线观看视频色| 免费大片黄手机在线观看| 亚洲av男天堂| 亚洲欧美一区二区三区国产| 国产一区二区三区综合在线观看 | 亚洲精品亚洲一区二区| 校园人妻丝袜中文字幕| av在线app专区| 色94色欧美一区二区| 亚洲激情五月婷婷啪啪| 国产一区二区在线观看日韩| 五月伊人婷婷丁香| 一级片'在线观看视频| 久久人人爽人人片av| 在线精品无人区一区二区三| 亚洲国产av新网站| 国产淫片久久久久久久久| 91精品一卡2卡3卡4卡| 久久国产乱子免费精品| 国产成人a∨麻豆精品| 日韩一区二区三区影片| 午夜影院在线不卡| 日日摸夜夜添夜夜添av毛片| 各种免费的搞黄视频| 好男人视频免费观看在线| 免费少妇av软件| 精品国产露脸久久av麻豆| 婷婷色麻豆天堂久久| 亚洲精品日韩在线中文字幕| 99精国产麻豆久久婷婷| freevideosex欧美| 午夜激情久久久久久久| 狂野欧美白嫩少妇大欣赏| 内射极品少妇av片p| 亚洲国产av新网站| 丝袜脚勾引网站| 欧美人与善性xxx| 99久久精品热视频| 国产精品99久久99久久久不卡 | 久久久精品94久久精品| 亚洲婷婷狠狠爱综合网| 两个人免费观看高清视频 | 观看免费一级毛片| 日日啪夜夜撸| 午夜老司机福利剧场| 国产在线男女| 久久久久精品久久久久真实原创| 久久久久网色| 免费黄频网站在线观看国产| 亚洲无线观看免费| 国产av精品麻豆| 精华霜和精华液先用哪个| 亚洲自偷自拍三级| 欧美3d第一页| 在线天堂最新版资源| 国产亚洲最大av| 国产爽快片一区二区三区| 制服丝袜香蕉在线| 校园人妻丝袜中文字幕| 精品一区二区三卡| 三级国产精品欧美在线观看| 草草在线视频免费看| 欧美国产精品一级二级三级 | 亚洲,欧美,日韩| 日韩制服骚丝袜av| a级毛色黄片| 美女脱内裤让男人舔精品视频| 亚洲欧洲精品一区二区精品久久久 | 亚洲成人手机| 一级av片app| 伊人久久国产一区二区| 亚洲欧美日韩东京热| 3wmmmm亚洲av在线观看| av女优亚洲男人天堂| 人妻制服诱惑在线中文字幕| 国产精品欧美亚洲77777| 女的被弄到高潮叫床怎么办| 高清黄色对白视频在线免费看 | 国产黄色视频一区二区在线观看| 国产日韩欧美在线精品| 日本黄大片高清| 亚洲内射少妇av| 三级国产精品片| 精品人妻熟女毛片av久久网站| 亚洲精品视频女| 国产亚洲最大av| 免费在线观看成人毛片| 午夜免费观看性视频| 亚洲av二区三区四区| 少妇人妻 视频| 毛片一级片免费看久久久久| 中文字幕亚洲精品专区| 高清黄色对白视频在线免费看 | 美女脱内裤让男人舔精品视频| 国产日韩一区二区三区精品不卡 | 97超碰精品成人国产| 精品久久久久久久久亚洲| 中文字幕av电影在线播放| 色视频www国产| 亚洲美女搞黄在线观看| 成人午夜精彩视频在线观看| 国产黄色免费在线视频| 在线观看国产h片| av国产精品久久久久影院| 亚洲国产精品专区欧美| 少妇裸体淫交视频免费看高清| 免费看日本二区| 欧美日本中文国产一区发布| 一区在线观看完整版| 日本爱情动作片www.在线观看| 青春草国产在线视频| 国产在线男女| 国产精品国产三级专区第一集| 青春草亚洲视频在线观看| 在线亚洲精品国产二区图片欧美 | 国产白丝娇喘喷水9色精品| av免费在线看不卡| 亚洲欧美清纯卡通| 91成人精品电影| 亚洲经典国产精华液单| 一区在线观看完整版| 成人漫画全彩无遮挡| 欧美xxxx性猛交bbbb| 欧美激情国产日韩精品一区| 国产伦精品一区二区三区视频9| 人妻人人澡人人爽人人| 青春草视频在线免费观看| 纵有疾风起免费观看全集完整版| 精品久久国产蜜桃| 国产69精品久久久久777片| 亚洲国产日韩一区二区| 亚洲av成人精品一二三区| 国产黄色视频一区二区在线观看| 免费观看av网站的网址| 亚洲婷婷狠狠爱综合网| 高清视频免费观看一区二区| 多毛熟女@视频| 日本91视频免费播放| 亚洲天堂av无毛| 狂野欧美白嫩少妇大欣赏| 肉色欧美久久久久久久蜜桃| av卡一久久| 日日摸夜夜添夜夜添av毛片| 多毛熟女@视频| 老女人水多毛片| av免费观看日本| 免费黄色在线免费观看| 交换朋友夫妻互换小说| 高清欧美精品videossex| 精品酒店卫生间| 丰满迷人的少妇在线观看| 免费久久久久久久精品成人欧美视频 | 欧美精品一区二区免费开放| 欧美xxxx性猛交bbbb| 免费观看在线日韩| 99久久综合免费| 欧美人与善性xxx| 欧美少妇被猛烈插入视频| 一级爰片在线观看| 最近最新中文字幕免费大全7| 美女视频免费永久观看网站| 国产精品成人在线| 国产乱来视频区| 高清黄色对白视频在线免费看 | 亚洲激情五月婷婷啪啪| 曰老女人黄片| 日本-黄色视频高清免费观看| 热re99久久国产66热| 内射极品少妇av片p| 国产亚洲最大av| 欧美97在线视频| 少妇人妻久久综合中文| 欧美成人午夜免费资源| 久久久久精品久久久久真实原创| 晚上一个人看的免费电影| 亚洲一区二区三区欧美精品| 国产亚洲精品久久久com| h日本视频在线播放| 在线观看美女被高潮喷水网站| 午夜视频国产福利| 一边亲一边摸免费视频| 在线观看人妻少妇| 日韩大片免费观看网站| 十八禁高潮呻吟视频 | 日本欧美视频一区| 各种免费的搞黄视频| 日日啪夜夜爽| 久久久久久久久大av| 亚洲精品国产av蜜桃| 欧美精品一区二区免费开放| 免费看av在线观看网站| 一区二区三区精品91| 国产 精品1| 天堂中文最新版在线下载| 99久久综合免费| 黄色一级大片看看| 国产成人免费无遮挡视频| 久久精品国产亚洲av涩爱| 亚洲国产精品一区三区| 人人妻人人看人人澡| 男女边吃奶边做爰视频| 久久久精品免费免费高清| 中文精品一卡2卡3卡4更新| 在线观看人妻少妇| 一本色道久久久久久精品综合| 丝瓜视频免费看黄片| 日本色播在线视频| 大香蕉97超碰在线| 亚洲,欧美,日韩| 国产成人免费观看mmmm| 麻豆成人午夜福利视频| 国产成人一区二区在线| 日本与韩国留学比较| 精品国产国语对白av| 精品一品国产午夜福利视频| 国产精品秋霞免费鲁丝片| 国产精品伦人一区二区| 日本av免费视频播放| 一级黄片播放器| 日韩av不卡免费在线播放| 久久久久精品性色| av网站免费在线观看视频| 精品一区二区三区视频在线| 我的老师免费观看完整版| 欧美 亚洲 国产 日韩一| 国产精品久久久久久久电影| 大香蕉97超碰在线| 交换朋友夫妻互换小说| 51国产日韩欧美| 精品人妻偷拍中文字幕| 亚洲欧美日韩另类电影网站| 国产精品福利在线免费观看| 妹子高潮喷水视频| 看非洲黑人一级黄片| 亚洲精品456在线播放app| 国产 一区精品| 亚洲精品中文字幕在线视频 | 欧美精品亚洲一区二区| 亚洲精品乱码久久久久久按摩| 国产黄色视频一区二区在线观看| 日产精品乱码卡一卡2卡三| 亚洲精品成人av观看孕妇| 亚洲精品久久久久久婷婷小说| 最近的中文字幕免费完整| 亚洲精品日韩在线中文字幕| 免费在线观看成人毛片| 人妻夜夜爽99麻豆av| 日日爽夜夜爽网站| a级一级毛片免费在线观看| 极品人妻少妇av视频| 性色av一级| 少妇人妻一区二区三区视频| 久久久久久人妻| 亚洲欧洲国产日韩| 观看免费一级毛片| 国产男女内射视频| 免费看不卡的av| 一级二级三级毛片免费看| 欧美+日韩+精品| 观看美女的网站| 亚洲国产日韩一区二区| 精品国产国语对白av| 精品国产露脸久久av麻豆| 高清欧美精品videossex| 寂寞人妻少妇视频99o| 老司机亚洲免费影院| 亚洲av欧美aⅴ国产| 黄色视频在线播放观看不卡| 只有这里有精品99| 中文字幕制服av| 一区二区av电影网| 18禁动态无遮挡网站| 在线观看av片永久免费下载| 最近手机中文字幕大全| 高清黄色对白视频在线免费看 | 亚洲精品日韩av片在线观看| 久久精品久久久久久久性| 成人毛片60女人毛片免费| 欧美日韩国产mv在线观看视频| 久久久久久久久久成人| 色视频在线一区二区三区| 黄色配什么色好看| 最后的刺客免费高清国语| 五月开心婷婷网| 日韩成人av中文字幕在线观看| 亚洲精华国产精华液的使用体验| 天美传媒精品一区二区| 亚洲欧美日韩东京热| 男男h啪啪无遮挡| 人体艺术视频欧美日本| 91精品国产国语对白视频| 在线精品无人区一区二区三| 精品少妇内射三级| 亚洲av日韩在线播放| 欧美xxxx性猛交bbbb| 自拍偷自拍亚洲精品老妇| 如何舔出高潮| 国产精品无大码| 久久人人爽人人片av| 国产午夜精品久久久久久一区二区三区| 国产精品久久久久久久久免| 久久久久久伊人网av| 国产色婷婷99| 亚洲国产精品一区三区| 看十八女毛片水多多多| 高清欧美精品videossex| 成人毛片a级毛片在线播放| 免费黄色在线免费观看| 3wmmmm亚洲av在线观看| 老女人水多毛片| 色哟哟·www| 久久av网站| av国产精品久久久久影院| 国产男女超爽视频在线观看| 一本一本综合久久| 欧美亚洲 丝袜 人妻 在线| 91精品国产九色| 极品人妻少妇av视频| 黄色欧美视频在线观看| 欧美性感艳星| 欧美最新免费一区二区三区| 国产在线一区二区三区精| 色婷婷av一区二区三区视频| h视频一区二区三区| 亚洲国产精品一区三区| 黄色日韩在线| 亚洲欧美中文字幕日韩二区| 国产精品无大码| av线在线观看网站| 自拍欧美九色日韩亚洲蝌蚪91 | 亚洲成人手机| 边亲边吃奶的免费视频| 又大又黄又爽视频免费| 久久久久久久大尺度免费视频| 国产熟女午夜一区二区三区 | 少妇精品久久久久久久| 少妇人妻久久综合中文| 老熟女久久久| 亚洲怡红院男人天堂| 亚洲精品乱码久久久v下载方式| av福利片在线| 香蕉精品网在线| 在线观看免费日韩欧美大片 | 曰老女人黄片| 在线天堂最新版资源| 美女cb高潮喷水在线观看| 伦精品一区二区三区| 丰满人妻一区二区三区视频av| 乱码一卡2卡4卡精品| 亚洲综合色惰| 精品亚洲成国产av| 一本—道久久a久久精品蜜桃钙片| 极品少妇高潮喷水抽搐| 黑人猛操日本美女一级片| 尾随美女入室| 综合色丁香网| 精品人妻熟女毛片av久久网站| 天美传媒精品一区二区| 久久精品国产亚洲av涩爱| 日韩成人av中文字幕在线观看| 97超碰精品成人国产| 一级毛片电影观看| a级一级毛片免费在线观看| 成人特级av手机在线观看| a级一级毛片免费在线观看| 夫妻午夜视频| 80岁老熟妇乱子伦牲交| √禁漫天堂资源中文www| 免费人成在线观看视频色| 欧美少妇被猛烈插入视频| 黄色一级大片看看| 曰老女人黄片| 少妇被粗大猛烈的视频| 成人综合一区亚洲| 观看美女的网站| 你懂的网址亚洲精品在线观看| 亚洲精品第二区| 亚洲欧洲国产日韩| 成人午夜精彩视频在线观看| 毛片一级片免费看久久久久| a级毛片在线看网站| 亚洲人成网站在线播| 一边亲一边摸免费视频| 欧美激情极品国产一区二区三区 | 日韩中字成人| 亚洲人成网站在线播| 国产熟女午夜一区二区三区 | 极品少妇高潮喷水抽搐| 亚洲av.av天堂| 国产一区有黄有色的免费视频| xxx大片免费视频| 欧美激情国产日韩精品一区| 大香蕉97超碰在线| 欧美日韩视频精品一区| 黄色欧美视频在线观看| 一区二区三区四区激情视频| 国产免费一级a男人的天堂| 香蕉精品网在线| 七月丁香在线播放| 精品久久久噜噜| 国产精品欧美亚洲77777| 18+在线观看网站| av天堂久久9| 亚洲,欧美,日韩| av国产精品久久久久影院| 伦理电影大哥的女人| 亚洲婷婷狠狠爱综合网| 国产av码专区亚洲av| 国内精品宾馆在线| 永久网站在线| 国产精品国产三级国产专区5o| 日本vs欧美在线观看视频 | 青青草视频在线视频观看| 看免费成人av毛片| 丰满少妇做爰视频| 夜夜看夜夜爽夜夜摸| 我要看黄色一级片免费的| 亚洲第一区二区三区不卡| av女优亚洲男人天堂| 久久久午夜欧美精品| 成人二区视频| 久久久久精品久久久久真实原创| 99久久人妻综合| 精品午夜福利在线看| 欧美老熟妇乱子伦牲交| 午夜福利影视在线免费观看| 国产一区二区三区av在线| 日本欧美视频一区| 免费观看无遮挡的男女| 国内少妇人妻偷人精品xxx网站| 亚洲不卡免费看| 久久ye,这里只有精品| 插逼视频在线观看| 超碰97精品在线观看| 国产一区亚洲一区在线观看| 久久鲁丝午夜福利片| 精品人妻熟女av久视频| 精品国产一区二区久久| 欧美高清成人免费视频www| 一区二区三区乱码不卡18| 深夜a级毛片| 美女中出高潮动态图| 国产成人a∨麻豆精品| 亚洲精品久久久久久婷婷小说| av在线老鸭窝| 日日爽夜夜爽网站| 国产黄色免费在线视频| 亚洲国产毛片av蜜桃av| 亚洲,一卡二卡三卡| 男女啪啪激烈高潮av片| 中文字幕人妻熟人妻熟丝袜美| 国产综合精华液| 欧美精品一区二区大全| 精品亚洲乱码少妇综合久久| 99九九线精品视频在线观看视频| 亚洲av日韩在线播放| 高清毛片免费看| 成人影院久久| 大话2 男鬼变身卡| 免费黄网站久久成人精品| 久久午夜综合久久蜜桃| 大话2 男鬼变身卡| 丝瓜视频免费看黄片| 亚洲情色 制服丝袜| 亚洲精品aⅴ在线观看| 亚洲精品一区蜜桃| 黄色怎么调成土黄色| 日韩人妻高清精品专区| 日本91视频免费播放| 亚洲色图综合在线观看| 亚洲精品乱码久久久久久按摩| 国产高清三级在线| 日日摸夜夜添夜夜爱| 国产精品国产三级专区第一集| 精品久久久噜噜| 大片免费播放器 马上看| 91aial.com中文字幕在线观看| 中文字幕久久专区| 精品国产国语对白av| av免费观看日本| 蜜桃在线观看..| 亚洲精品日本国产第一区| 日本vs欧美在线观看视频 |