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

    散斑能見度光譜法測量筒倉內顆粒流的顆粒溫度

    2022-07-06 08:16:06陳泉鄭澤希李然孫其誠楊暉
    化工學報 2022年6期
    關鍵詞:散斑流率筒倉

    陳泉,鄭澤希,李然,3,孫其誠,楊暉

    (1 上海理工大學光電信息與計算機工程學院,上海 200093; 2 上海理工大學機械工程學院,上海 200093; 3 上海理工大學醫(yī)療器械與食品學院,上海 200093; 4 清華大學水沙科學與水利水電工程國家重點實驗室,北京 100084)

    引 言

    筒倉作為常見的顆粒物料存儲設備廣泛應用于糧食存儲、醫(yī)藥加工和化工生產等領域[1-4]。理想的筒倉能夠在最小的占地面積上存儲最多的顆粒,并且以所需的質量流率正常流出[5]。然而,在運行時,由于離散顆粒的無序運動導致了筒倉內固體顆粒無法按照質量流(mass flow)形式正常流出,引起了筒壁與顆粒之間的黏滑運動[6],間歇性堵塞[7]以及非局域流變[8]等問題。這些問題嚴重影響了產品的生產質量,甚至危害筒倉的安全運行。因此,測量筒倉內顆粒的無序運動,掌握其變化規(guī)律,在筒倉設計和化工調控工程中具有重要意義。

    顆粒溫度(granular temperature)反映了介尺度條件下顆粒運動的無序程度[9]。Ogawa 等[10]根據(jù)微觀粒子的熱運動提出了顆粒溫度的概念來表征顆粒系統(tǒng)的無序運動,用δv2~<δvδv>表示,其中δv表示無規(guī)則運動顆粒的速度波動,<>表示系綜平均。顆粒溫度與熱力學溫度相似,二者都產生壓力,并且控制著顆粒系統(tǒng)的質量、動量和能量傳播。近20年來,人們采用離散單元法(discrete element method,DEM)對筒倉內顆粒流的顆粒溫度進行了大量研究,建立了介觀顆粒溫度與宏觀顆粒流型、質量流率之間的關系模型。2007 年,Ahn[11]采用基于硬球接觸模型的DEM 方法模擬不同質量流率條件下圓柱形筒倉內球形顆粒的運動,根據(jù)顆粒的矢量速度計算了筒倉內顆粒溫度的軸向分布,揭示了顆粒動能的耗散率與顆粒溫度值滿足指數(shù)為1.5 的冪函數(shù)關系。2013 年,Tewari 等[12]采用基于顆粒動理學的DEM 方法模擬二維筒倉中球形顆粒的運動,探究堵塞前顆粒流的動力學特征。根據(jù)單位時間步長內顆粒的位移計算出筒倉內顆粒流的顆粒溫度,發(fā)現(xiàn)了當質量流率較高時,筒倉中心區(qū)域的顆粒溫度值較低;而當顆粒流接近堵塞時,顆粒溫度的分布發(fā)生了明顯逆轉,在筒壁處的顆粒溫度存在極小值。2015 年,Regele 等[13]采用基于靜摩擦接觸模型的DEM 方法模擬球形顆粒通過具有狹窄通道的三維筒倉,統(tǒng)計了單元顆粒的平移速度和旋轉速度,計算了不同流態(tài)區(qū)域的顆粒溫度分布。研究表明:孔口附近快速顆粒流區(qū)域內存在較低的顆粒溫度值,但該位置處顆粒的旋轉速度較高。擴散波光譜技術(diffusing wave spectroscopy,DWS)和散斑能見度光譜技術(speckle visibility spectroscopy,SVS)是兩種能夠直接測量顆粒溫度的方法[14]。DWS 技術要求被測量顆粒系統(tǒng)的運動狀態(tài)具有穩(wěn)定性,以滿足光子在自相關時間內的擴散傳播過程;SVS 技術是在DWS技術基礎上發(fā)展起來的,該測量技術的優(yōu)勢在于具有良好的魯棒性,適用于筒倉內復雜顆粒流的顆粒溫度測量,并且SVS 技術具有微秒級時間分辨率(由線陣CCD相機的采集幀率決定,最快可以到10 μs)和微米級空間分辨率(激光波長的1/4,約0.1 μm)。2005 年,Bandyopadhyay 等[15]根據(jù)擴散波光譜法首次提出了散斑能見度光譜技術,并應用于振動流化床的顆粒溫度測量,觀測到顆粒溫度的測量值與系統(tǒng)振動頻率具有相同的周期性,進而驗證了該技術能夠準確測量密集顆粒系統(tǒng)的無序運動。2010 年,Katsuragi 等[16]利用散斑能見度光譜技術測量了斜坡流中不同位置的顆粒溫度分布,通過對比顆粒溫度與顆粒堆積密度,發(fā)現(xiàn)了顆粒溫度隨顆粒堆積深度的變化是非線性的。2019 年,本課題組[17-19]應用該技術測量了轉筒內雪崩顆粒流的顆粒溫度,觀察到雪崩前顆粒堆積的重排現(xiàn)象,發(fā)現(xiàn)了不規(guī)則顆粒的重排運動可以減小雪崩的持續(xù)時間。

    因此,本文采用散斑能見度光譜技術,對卸料過程中筒倉內兩種粒徑的球形顆粒分別開展顆粒溫度測量。探究了穩(wěn)態(tài)流動過程中顆粒溫度的分布特征以及顆粒粒徑對顆粒溫度的影響。最后,分析了堵塞過程中顆粒溫度的變化規(guī)律。

    1 實驗材料和方法

    為了掌握筒倉內顆粒流的復雜動力學,本文搭建了基于線陣CCD 相機的散斑能見度光譜實驗裝置,測量筒倉內玻璃珠在卸料過程中的顆粒溫度。

    1.1 顆粒材料

    圖1(a)是顯微鏡下均值粒徑為0.94 mm 的玻璃珠顆粒。由于透明材料可以幫助觀察顆粒床中一定堆積深度(4~6 mm)的顆粒運動,因此采用透明玻璃珠進行筒倉內顆粒流的測量實驗。首先,采用標準篩網(孔徑分別為0.88 和1.00 mm;孔徑分別為1.40 和1.70 mm)對玻璃珠進行了兩次篩分處理,獲得了兩種不同粒徑的實驗樣品顆粒。然后,隨機取出部分樣品顆粒放置在顯微鏡(Olympus SZX16)的工作臺上,放大25 倍后進行圖像采集。最后,根據(jù)顆粒的投影面積計算顆粒的等效圓面積直徑[20]。連續(xù)統(tǒng)計300 顆玻璃珠樣品,獲得實驗顆粒粒徑的累積分布,如圖1(b)所示。表1為兩種粒徑玻璃珠的特征參數(shù)。表中β是顆粒堆積的下休止角[21]。

    表1 實驗顆粒材料的特征參數(shù)Table 1 Characteristic of granular materials used in this work

    圖1 實驗玻璃珠的圖片和顆粒粒徑的累積分布(fCDF代表顆粒粒徑的累積分布函數(shù);d代表顆粒粒徑)Fig.1 The image of glass beads and the cumulative particle size distributions of the glass beads

    1.2 實驗裝置系統(tǒng)

    圖2為實驗裝置實物圖。實驗裝置系統(tǒng)包括激光光源,凸透鏡,線陣CCD 相機,濾光片,電子天平秤和矩形筒倉。光源是波長為532 nm 的半導體激光器,輸出功率為100 mW,光束直徑為2 mm。出射激光經凸透鏡擴束得到直徑為8 mm 的入射光束。當入射光與堆積顆粒發(fā)生相互作用時,產生的前向散射光被線陣CCD 相機(加拿大DALSA 的Spyder mod 3)捕獲。線陣CCD 相機的像素尺寸為14 μm×14 μm,像元個數(shù)為1024 個,采樣幀率設置為20 kHz。為了減弱其他光源信號對實驗測量結果的影響,在線陣CCD相機的表面粘貼了532 nm的窄帶濾光片。平底矩形筒倉由長方體容器和金屬底座兩部分組成。長方體容器采用透明有機玻璃材料加工制成,壁厚為10 mm,高度和寬度分別為240 和150 mm,前后壁面的間隙為8 mm。筒倉底座由金屬鋁加工制成,表面采用電鍍工藝進行氧化黑處理??卓诰匦尾鄣膶挾葹? mm,利用可調插銷控制孔口矩形槽的長度D1,使得D1在[4 mm,18 mm]范圍內變化。電子天平秤(SF-400)通過RS232 接口與計算機相連接,用于記錄筒倉內填充顆粒的質量流率。電子秤的最大采樣頻率為5 Hz,重復性誤差為1 g。以O點為坐標原點建立笛卡兒坐標系,其中,O點位于孔口中心的內壁上,x軸與y軸形成的平面與筒倉的前壁面平行,如圖2 中矩形筒倉的正視圖所示。平底矩形筒倉是關于中心垂線對稱的,因此在x?[0,75 mm],y?[0,240 mm]范圍內測量顆粒溫度。每個測量點的大小為直徑8 mm 的圓形區(qū)域(約40 個顆粒),兩個相鄰測量點之間的中心距離為10 mm。測量點A位于筒倉孔口附近的快速顆粒流區(qū)域。測量點A為直徑8 mm 的圓形區(qū)域,中心坐標為A(0,10 mm,0)。實驗前,先利用軟橡膠塞堵住孔口矩形槽,再從筒倉頂部采用均勻分散加載方式進行玻璃珠的加載。

    圖2 實驗裝置實物圖(D1代表筒倉孔口的長度)Fig.2 The annotated photo of the experimental setup

    1.3 散斑能見度光譜法

    圖3為散斑能見度光譜法測量顆粒溫度的原理圖。假設單個光子通過離散顆粒樣品時發(fā)生了N次散射,則根據(jù)擴散波光譜原理,在t時刻光子散射N次的總路徑長度S為[22]

    式中,ki(t)為t時刻光子經過i次散射后的波矢量;ri(t)為t時刻顆粒的空間位置坐標。同時,光子的相位可以表示為[23]

    式中,φ(t)為光子相位隨時間的變化函數(shù)。根據(jù)式(1)和式(2),得到散射光強E(t)的表達式為

    式中,ES為光程為S光子的光強幅值。

    因此,散射光場的歸一化自相關函數(shù)可以表示為

    式中,λ為入射激光的波長;δv= <δv2>,δv2表示顆粒溫度;Γ= 4πδv λ,是函數(shù)g1(t) 的衰減線寬。

    根據(jù)統(tǒng)計光學原理,散斑圖像的對比度可以用單位曝光時間內散射光強的方差σ與均值的比值來表示[24],即

    式中,T為采集相機的曝光時間。根據(jù)高斯光束的Siegert 公式,建立了散斑圖像對比度與散射光場自相關函數(shù)之間的關系[25]

    式中,β為測量系統(tǒng)的相干因子。將式(5)代入式(7),得

    顆粒溫度可以根據(jù)式(8)反演獲得。注意,由于受到CCD 相機的電路噪聲和測量環(huán)境的干擾,相干因子β無法被準確測量[15]。為了解決該問題,一般將2 倍和1 倍曝光時間下的散斑對比度值相除,消去系統(tǒng)相干因子,來獲得顆粒溫度δv2。另外,一般要求散斑面積的等效直徑約等于10 倍顆粒粒徑。

    由于光子的波粒二象性,在測量區(qū)域內不同光子的傳播會產生干涉現(xiàn)象,從而在CCD 相機上顯示出“顆粒狀”的散斑圖像,如圖3 中灰度圖所示。在相同的曝光時間條件下,顆粒的運動速度越快,散斑圖像的模糊程度越高,而散斑圖像的對比度就越低。圖3 中紅色曲線是根據(jù)圖3 中的散斑圖像計算獲得的顆粒溫度變化曲線。散斑能見度光譜測量技術的重復性誤差小于3%[26]。

    圖3 散斑能見度光譜法測量顆粒溫度的原理示意圖[ES(t)代表光程為S的散射光子在t時刻的光強幅值;ki代表入射激光經過i次散射后的波矢量]Fig.3 Schematic diagram of measuring granular temperature by speckle visibility spectroscopy

    2 實驗結果與討論

    設置孔口矩形槽的長度D1=12 mm,將均值粒徑為0.94 和1.55 mm 的玻璃珠加載到平底矩形筒倉內進行實驗。采用散斑能見度光譜法測量卸料過程中顆粒流的顆粒溫度,每個測量點重復測量10 次。由于筒倉內顆粒流的運動是關于筒倉中心軸線對稱的,因此在顆粒溫度測量過程中,僅測量筒倉內顆粒流的右半部分。

    2.1 顆粒溫度的特征

    圖4 為卸料過程中筒倉內測量點A的顆粒溫度變化曲線。可以看出,顆粒溫度曲線包含了初始瞬態(tài)流,穩(wěn)態(tài)顆粒流和排空三種不同的卸料階段。在穩(wěn)態(tài)顆粒流中,顆粒溫度值保持相對穩(wěn)定。該結果說明了卸料過程中離散顆粒的運動在介尺度條件下具有穩(wěn)定性;而在非穩(wěn)態(tài)顆粒流中,由于顆粒介質之間的摩擦和非彈性碰撞具有高度的非線性[27],導致了顆粒溫度曲線存在較大的梯度變化。另外,對比卸料過程中顆粒流的三種不同階段,穩(wěn)態(tài)顆粒流的持續(xù)時間最長。因此,在接下來的研究中重點關注穩(wěn)態(tài)顆粒流的顆粒溫度。

    圖4 筒倉內顆粒流測量點A(0,10 mm,0)的顆粒溫度變化曲線(代表穩(wěn)態(tài)流中粒徑為0.94 mm顆粒的平均顆粒溫度;代表穩(wěn)態(tài)流中粒徑為1.55 mm顆粒的平均顆粒溫度)Fig.4 The curve of granular temperature at measuring point A(0,10 mm,0)in silo

    表2 為穩(wěn)態(tài)顆粒流過程中筒倉內A位置顆粒溫度與質量流率的均值和方差??梢钥闯?,10 次實驗中顆粒溫度的測量結果具有良好的一致性。值得注意的是,均值粒徑為0.94 和1.55 mm 顆粒的質量流率分別為Q1=(33.5±0.8) g/s,Q2=(26.7±0.7) g/s;而顆粒溫度為δv21=(2.02±0.01) mm2/s2,δv22=(2.35±0.02)mm2/s2。這表明了大粒徑顆粒的無序運動程度更強,即大粒徑顆粒在卸料過程中具有更高的能量耗散,這可能是導致其質量流率較低的直接原因。該現(xiàn)象揭示了介觀顆粒的無序運動是控制宏觀質量流率的主要因素之一。在實際應用過程中,通過在筒倉內部添加障礙物調節(jié)離散顆粒的無序運動,可以為調控宏觀質量流率提供新的途徑。

    2.2 顆粒溫度分布

    圖5為兩種粒徑顆粒在筒倉內穩(wěn)態(tài)顆粒流中的顆粒溫度分布??梢钥闯?,不同粒徑顆粒的顆粒溫度分布具有相似性。顆粒溫度的最大值出現(xiàn)在筒倉頂部區(qū)域,并且以“收縮”形式向下遞減傳播,如圖5(a)和(b)所示。筒倉頂部的高顆粒溫度值可能是顆粒床的雪崩運動導致的。注意,兩種粒徑顆粒在筒倉孔口附近均存在一個低顆粒溫度區(qū)域。較低的顆粒溫度值揭示了該位置處離散顆粒存在定向有序的運動,這與斑團模型中介觀顆粒運動的特征相同,并且Zhang 等[28-29]在DEM 模擬實驗中也觀察到了類似的現(xiàn)象。進一步,根據(jù)Hu 等[30]的研究結果,筒倉孔口兩側的堆積顆粒可能是造成低顆粒溫度的直接原因??卓诟浇^(qū)的顆粒堆積可以減小離散顆粒間碰撞產生的接觸力,從而降低了顆粒系統(tǒng)的顆粒溫度。這表明在實際筒倉的優(yōu)化設計中,孔口兩側添加緩沖單元,可以減弱顆粒間碰撞產生的接觸力,從而使得填充顆??梢愿菀着懦?,有利于減小筒倉的堵塞概率。

    圖5 穩(wěn)態(tài)流動中筒倉內顆粒流的顆粒溫度分布Fig.5 Granular temperature distribution of steady flow in silo

    圖5(b)和(e)分別為均值粒徑為0.94 和1.55 mm顆粒的顆粒溫度的徑向分布。可以看出,兩種粒徑顆粒的顆粒溫度的徑向分布是基本相同的。在筒倉的不同高度y處,顆粒溫度的變化趨勢是一致的,顆粒溫度值隨著x的增大而減小。并且,隨著y值的減小,顆粒溫度為0 值的測量點逐漸增多。這表明了停滯區(qū)隨著筒倉高度的降低而增多。然而,在不同的軸向位置處,顆粒溫度的值隨著y的增大而增大,如圖5(c)和(f)所示。當測量點距離筒壁約10 mm時,顆粒溫度約為0值,這表明該位置處的顆粒均處于停滯區(qū)內。

    2.3 顆粒流的堵塞測量

    堵塞現(xiàn)象在筒倉顆粒流中是隨機產生的。圖6 為堵塞過程中筒倉內A位置的顆粒溫度變化曲線。可以看出,當t>0.06 s 時,顆粒流處于穩(wěn)態(tài)流狀態(tài),顆粒溫度值相對穩(wěn)定,在2.02 mm2/s2上下波動。然而,當t=0.06 s 時,筒倉內顆粒流發(fā)生了堵塞現(xiàn)象。在0.06 s 時間內,顆粒溫度值從2.02 mm2/s2減小至0。這表明了筒倉內顆粒流的堵塞現(xiàn)象存在一個弛豫過程。進一步,在誤差小于2%的條件下,采用最小二乘法對堵塞前0.06 s 時間內的顆粒溫度進行擬合處理,結果如圖6 中插圖所示??梢钥闯?,堵塞過程中的顆粒溫度變化滿足指數(shù)函數(shù)分布。

    圖6 堵塞過程中均值粒徑為0.94 mm顆粒在筒倉內A位置的顆粒溫度變化曲線(插圖為堵塞前0.06 s時間內的顆粒溫度曲線;藍色實線代表擬合曲線;Δt表示堵塞過程的弛豫時間)Fig.6 The curve of granular temperature with an average particle size of 0.94 mm at position A in silos during blockage

    圖7是筒倉內顆粒流堵塞過程的弛豫時間Δt隨孔口矩形槽長度D1的變化曲線??梢钥闯觯谙嗤卓诔叽鐥l件下,均值粒徑為1.55 mm 顆粒的堵塞弛豫時間Δt更長。觀察到D1在[4 mm, 18 mm]范圍內時,兩種粒徑顆粒的堵塞弛豫時間Δt與孔口尺寸D1均滿足正比例關系,最短弛豫時間 Δt=(0.020 ±0.002)s。該結果表明通過測量筒倉內顆粒流的介觀顆粒溫度信號,可以提前預測顆粒流的堵塞事件,并且隨著筒倉孔口尺寸的增加,可以更早預測顆粒流的堵塞。

    圖7 筒倉內顆粒流堵塞過程的弛豫時間Δt隨孔口矩形槽長度D1的變化曲線Fig.7 The curve of relaxation time of blockage changed with orifice size in silo

    3 結 論

    基于散斑能見度光譜技術對卸料過程中筒倉內顆粒流的顆粒溫度進行了測量,對比了均值粒徑為0.94 和1.55 mm 兩種玻璃珠的時變顆粒溫度,探究宏觀質量流率與介觀顆粒溫度之間的聯(lián)系。同時,利用該測量技術高時空分辨率的特點,分析了筒倉顆粒流堵塞過程的弛豫變化特征。主要結論如下。

    (1)對比均值粒徑為0.94 和1.55 mm 顆粒在筒倉內穩(wěn)態(tài)流動中的質量流率和顆粒溫度,發(fā)現(xiàn)了質量流率較大的顆粒,其顆粒溫度值較低。該現(xiàn)象揭示了介觀顆粒的無序運動是控制宏觀質量流率的主要因素之一。

    (2)發(fā)現(xiàn)了筒倉孔口附近存在一個低顆粒溫度的特征區(qū)域,說明了該位置處的離散顆粒存在定向有序的運動。筒倉孔口兩側的堆積顆??梢杂行p弱顆粒間碰撞產生的接觸力是造成該現(xiàn)象的主要原因。

    (3)不同粒徑顆粒的顆粒溫度分布具有一致性。顆粒溫度的最大值出現(xiàn)在筒倉頂部區(qū)域,并且以“收縮”形式向下擴散傳播。在筒倉的不同高度y處,顆粒溫度的變化規(guī)律是相同的,隨著x的增大而減小。然而,在不同的軸向位置處,顆粒溫度的值隨著y的增大而增大。

    (4)利用散斑能見度光譜技術的高時間分辨率的特點,確定了筒倉內顆粒流堵塞的弛豫過程,發(fā)現(xiàn)了堵塞過程中的顆粒溫度變化滿足指數(shù)函數(shù)分布。并且弛豫時間與筒倉孔口尺寸呈正比例關系。該結果有助于揭示顆粒流隨機堵塞的預兆。

    本文實驗結果揭示了筒倉內顆粒流的運動規(guī)律,為完善顆粒材料的傳輸效率和化工過程的精準控制提供了參考數(shù)據(jù)。

    猜你喜歡
    散斑流率筒倉
    筒倉施工中滑模技術的應用實踐探討
    CDN 緩存空間調整優(yōu)化方法
    帶交換和跳躍的一維雙向自驅動系統(tǒng)的仿真研究
    激光顯示中的彩色散斑測量研究
    激光投影顯示散斑抑制方法研究
    多束光在皮膚組織中傳輸?shù)拿商乜_模擬
    用于檢驗散斑協(xié)方差矩陣估計性能的白化度評價方法
    雷達學報(2017年3期)2018-01-19 02:01:21
    聚合物氣輔共擠中熔體體積流率對擠出脹大的影響
    立筒倉通風效果的探索
    大型筒倉布料與給料設施的應用現(xiàn)狀與選型研究
    浙江電力(2016年7期)2016-04-11 12:25:46
    亚洲av欧美aⅴ国产| 人妻久久中文字幕网| 日日摸夜夜添夜夜添小说| 在线 av 中文字幕| 极品少妇高潮喷水抽搐| 1024香蕉在线观看| 成年美女黄网站色视频大全免费| 女警被强在线播放| 免费久久久久久久精品成人欧美视频| 日韩免费高清中文字幕av| 999久久久国产精品视频| 午夜福利一区二区在线看| 在线看a的网站| 免费不卡黄色视频| 精品高清国产在线一区| 国产一级毛片在线| 18禁国产床啪视频网站| 高清黄色对白视频在线免费看| 亚洲国产av影院在线观看| 国产av国产精品国产| 久久久精品区二区三区| 考比视频在线观看| 纵有疾风起免费观看全集完整版| 日韩一卡2卡3卡4卡2021年| 欧美中文综合在线视频| av一本久久久久| 午夜成年电影在线免费观看| 国产不卡av网站在线观看| 日韩欧美免费精品| 成人黄色视频免费在线看| 777米奇影视久久| 欧美精品av麻豆av| 啦啦啦视频在线资源免费观看| 国产精品1区2区在线观看. | 欧美av亚洲av综合av国产av| 久久久水蜜桃国产精品网| 久久99一区二区三区| 国产精品秋霞免费鲁丝片| 中文字幕人妻熟女乱码| 国产欧美日韩综合在线一区二区| 国产主播在线观看一区二区| 蜜桃国产av成人99| 不卡一级毛片| 老熟妇仑乱视频hdxx| 精品亚洲乱码少妇综合久久| 国产又爽黄色视频| 动漫黄色视频在线观看| 亚洲精品国产区一区二| 又紧又爽又黄一区二区| 久久亚洲精品不卡| 国产免费福利视频在线观看| 少妇的丰满在线观看| 一本久久精品| videos熟女内射| 国产男女超爽视频在线观看| 国产成人av教育| 黄色毛片三级朝国网站| 久久久久久人人人人人| 久久av网站| 成人手机av| 亚洲欧洲精品一区二区精品久久久| 日韩 亚洲 欧美在线| 久久毛片免费看一区二区三区| a级毛片在线看网站| 夫妻午夜视频| 九色亚洲精品在线播放| 久久亚洲精品不卡| 亚洲激情五月婷婷啪啪| svipshipincom国产片| 亚洲欧美精品综合一区二区三区| 久久久精品免费免费高清| 永久免费av网站大全| 午夜激情久久久久久久| 丝袜脚勾引网站| 两个人免费观看高清视频| 最近最新中文字幕大全免费视频| 婷婷丁香在线五月| 成年人午夜在线观看视频| 美女中出高潮动态图| 欧美精品一区二区大全| 一级片'在线观看视频| 国产精品久久久久久精品电影小说| 男女无遮挡免费网站观看| 纵有疾风起免费观看全集完整版| svipshipincom国产片| 亚洲欧美激情在线| 久久av网站| 久久香蕉激情| 一级毛片精品| 国产深夜福利视频在线观看| 午夜福利在线免费观看网站| 欧美日韩福利视频一区二区| 91精品伊人久久大香线蕉| 精品亚洲乱码少妇综合久久| 久久人妻熟女aⅴ| 大码成人一级视频| 成年人午夜在线观看视频| 国产免费av片在线观看野外av| 亚洲精品日韩在线中文字幕| 欧美久久黑人一区二区| 蜜桃国产av成人99| 女人久久www免费人成看片| 亚洲国产av新网站| 一本综合久久免费| 在线观看一区二区三区激情| 人人妻人人澡人人爽人人夜夜| 久久精品成人免费网站| 51午夜福利影视在线观看| 日韩有码中文字幕| 久久青草综合色| 老汉色∧v一级毛片| 午夜福利乱码中文字幕| a级毛片在线看网站| 久久国产精品人妻蜜桃| 12—13女人毛片做爰片一| 女性生殖器流出的白浆| 亚洲欧美成人综合另类久久久| 欧美久久黑人一区二区| 久久久久久人人人人人| 叶爱在线成人免费视频播放| 国产亚洲欧美精品永久| 亚洲专区国产一区二区| 精品少妇黑人巨大在线播放| 国产免费福利视频在线观看| 桃花免费在线播放| av电影中文网址| 巨乳人妻的诱惑在线观看| 国产欧美亚洲国产| 999久久久精品免费观看国产| 午夜福利乱码中文字幕| 日韩人妻精品一区2区三区| 中国国产av一级| 超碰97精品在线观看| 国产亚洲欧美在线一区二区| 成人国语在线视频| av天堂在线播放| 欧美xxⅹ黑人| 一级毛片电影观看| 热re99久久精品国产66热6| 五月开心婷婷网| 99久久精品国产亚洲精品| 免费日韩欧美在线观看| 久久久精品免费免费高清| 成年人黄色毛片网站| 91精品伊人久久大香线蕉| tube8黄色片| 伊人亚洲综合成人网| 国产一区二区 视频在线| 久久人人97超碰香蕉20202| 久久人妻福利社区极品人妻图片| 欧美日韩亚洲国产一区二区在线观看 | 久久热在线av| 999久久久国产精品视频| 中文字幕另类日韩欧美亚洲嫩草| 国产亚洲欧美精品永久| 激情视频va一区二区三区| 亚洲精品国产精品久久久不卡| 视频区图区小说| 一本大道久久a久久精品| 午夜福利乱码中文字幕| 成人黄色视频免费在线看| 天天躁日日躁夜夜躁夜夜| 一区福利在线观看| 成人手机av| 天堂8中文在线网| 久久国产精品影院| 国产又色又爽无遮挡免| 久久久久国内视频| 国产成人啪精品午夜网站| 十分钟在线观看高清视频www| 欧美另类亚洲清纯唯美| 搡老熟女国产l中国老女人| 每晚都被弄得嗷嗷叫到高潮| 国产av一区二区精品久久| 男女之事视频高清在线观看| 日本撒尿小便嘘嘘汇集6| 青春草视频在线免费观看| 欧美激情高清一区二区三区| 国产成人av教育| 中文字幕av电影在线播放| 妹子高潮喷水视频| 精品乱码久久久久久99久播| 视频在线观看一区二区三区| 熟女少妇亚洲综合色aaa.| 亚洲午夜精品一区,二区,三区| 免费黄频网站在线观看国产| 黑丝袜美女国产一区| 99热国产这里只有精品6| 窝窝影院91人妻| 国产在线观看jvid| 国产黄频视频在线观看| 99九九在线精品视频| 亚洲国产看品久久| 国产欧美日韩一区二区精品| 欧美人与性动交α欧美软件| 1024香蕉在线观看| 99热网站在线观看| 亚洲精品成人av观看孕妇| 欧美日韩中文字幕国产精品一区二区三区 | 汤姆久久久久久久影院中文字幕| 黄色毛片三级朝国网站| 日韩精品免费视频一区二区三区| 欧美成狂野欧美在线观看| 咕卡用的链子| 黄片小视频在线播放| 国产一区二区在线观看av| 国产精品一区二区精品视频观看| 老汉色av国产亚洲站长工具| 狂野欧美激情性xxxx| 国产亚洲一区二区精品| 99国产综合亚洲精品| 成人影院久久| 老熟妇仑乱视频hdxx| 免费高清在线观看日韩| 电影成人av| 精品人妻熟女毛片av久久网站| 性色av乱码一区二区三区2| 老司机午夜福利在线观看视频 | 亚洲天堂av无毛| 人人澡人人妻人| 天堂俺去俺来也www色官网| 日本精品一区二区三区蜜桃| 国产黄频视频在线观看| 看免费av毛片| 三上悠亚av全集在线观看| 国产三级黄色录像| 亚洲欧美精品自产自拍| 国产精品秋霞免费鲁丝片| 两个人看的免费小视频| 欧美成人午夜精品| a级片在线免费高清观看视频| 国产黄频视频在线观看| 午夜免费鲁丝| 女人被躁到高潮嗷嗷叫费观| 人人妻人人爽人人添夜夜欢视频| 久久九九热精品免费| 久久国产精品人妻蜜桃| 不卡一级毛片| 19禁男女啪啪无遮挡网站| 精品国产乱码久久久久久小说| 91精品伊人久久大香线蕉| 久久精品国产a三级三级三级| 亚洲性夜色夜夜综合| 国产区一区二久久| 无限看片的www在线观看| 日本欧美视频一区| 纵有疾风起免费观看全集完整版| 日日爽夜夜爽网站| 两个人免费观看高清视频| 少妇人妻久久综合中文| 国产日韩欧美视频二区| 欧美xxⅹ黑人| 少妇猛男粗大的猛烈进出视频| 国产av一区二区精品久久| 国产成人免费观看mmmm| 老司机靠b影院| 69av精品久久久久久 | 国产深夜福利视频在线观看| 日韩大码丰满熟妇| 女人精品久久久久毛片| 超碰97精品在线观看| 精品福利永久在线观看| 丰满少妇做爰视频| 久久久久精品人妻al黑| 一级毛片电影观看| av线在线观看网站| 日本a在线网址| 男女高潮啪啪啪动态图| 天堂俺去俺来也www色官网| 国产成人精品在线电影| 国产97色在线日韩免费| 国产亚洲av高清不卡| 黑人猛操日本美女一级片| 新久久久久国产一级毛片| 波多野结衣一区麻豆| 国产精品一区二区精品视频观看| 久久中文字幕一级| 50天的宝宝边吃奶边哭怎么回事| 男人添女人高潮全过程视频| 人人妻人人澡人人看| 欧美成人午夜精品| 永久免费av网站大全| 精品国内亚洲2022精品成人 | 亚洲欧美精品综合一区二区三区| 久久人妻福利社区极品人妻图片| 啦啦啦视频在线资源免费观看| 亚洲国产看品久久| 色婷婷av一区二区三区视频| 热re99久久精品国产66热6| 久久久久久久久久久久大奶| 一本大道久久a久久精品| 国产亚洲欧美在线一区二区| 9色porny在线观看| 成年动漫av网址| 操美女的视频在线观看| 老熟妇仑乱视频hdxx| 亚洲五月婷婷丁香| 狂野欧美激情性bbbbbb| 精品国产乱码久久久久久男人| 国产成人免费无遮挡视频| 久热这里只有精品99| 狂野欧美激情性xxxx| 日韩欧美国产一区二区入口| 99香蕉大伊视频| 一区二区三区激情视频| 欧美精品一区二区大全| 老司机深夜福利视频在线观看 | 9色porny在线观看| 午夜日韩欧美国产| 婷婷丁香在线五月| 欧美日韩av久久| 中文字幕精品免费在线观看视频| 亚洲av日韩在线播放| 中文字幕最新亚洲高清| 激情视频va一区二区三区| 久久久国产精品麻豆| 可以免费在线观看a视频的电影网站| 亚洲成人免费电影在线观看| 99国产综合亚洲精品| 久久影院123| 黑人猛操日本美女一级片| 亚洲色图综合在线观看| 狠狠精品人妻久久久久久综合| 国产免费现黄频在线看| 2018国产大陆天天弄谢| 亚洲三区欧美一区| 91成年电影在线观看| av网站在线播放免费| 黄色 视频免费看| 日本黄色日本黄色录像| 桃花免费在线播放| 女性被躁到高潮视频| 巨乳人妻的诱惑在线观看| 亚洲五月婷婷丁香| 搡老熟女国产l中国老女人| 考比视频在线观看| 日本一区二区免费在线视频| 国产在视频线精品| 嫩草影视91久久| 亚洲黑人精品在线| 两个人看的免费小视频| 丁香六月天网| 免费高清在线观看视频在线观看| 精品卡一卡二卡四卡免费| av视频免费观看在线观看| 国产精品秋霞免费鲁丝片| 老司机在亚洲福利影院| 岛国毛片在线播放| 国产亚洲欧美精品永久| 国产高清视频在线播放一区 | 成人国语在线视频| 一级黄色大片毛片| 欧美黑人欧美精品刺激| 黄色片一级片一级黄色片| a级毛片黄视频| 亚洲精品美女久久av网站| 欧美中文综合在线视频| 亚洲精品国产av蜜桃| 亚洲一码二码三码区别大吗| 亚洲情色 制服丝袜| 国产在线一区二区三区精| 交换朋友夫妻互换小说| 男男h啪啪无遮挡| 蜜桃在线观看..| 国产成人一区二区三区免费视频网站| 水蜜桃什么品种好| 国产精品九九99| 免费久久久久久久精品成人欧美视频| 成人av一区二区三区在线看 | 无限看片的www在线观看| 国产一卡二卡三卡精品| 久久久久久免费高清国产稀缺| 激情视频va一区二区三区| 巨乳人妻的诱惑在线观看| 久久精品熟女亚洲av麻豆精品| 亚洲精品国产av蜜桃| 性色av一级| 不卡av一区二区三区| 欧美日韩国产mv在线观看视频| 欧美乱码精品一区二区三区| 欧美人与性动交α欧美软件| 精品少妇久久久久久888优播| 国产成人免费观看mmmm| 国产欧美日韩一区二区精品| 大香蕉久久成人网| 少妇裸体淫交视频免费看高清 | 一本大道久久a久久精品| 国产日韩一区二区三区精品不卡| 日本vs欧美在线观看视频| 成人免费观看视频高清| 国产亚洲精品久久久久5区| 国产真人三级小视频在线观看| 秋霞在线观看毛片| 女人被躁到高潮嗷嗷叫费观| 日韩欧美国产一区二区入口| 丝袜喷水一区| 美女高潮喷水抽搐中文字幕| 亚洲精品中文字幕在线视频| 搡老熟女国产l中国老女人| 亚洲欧美日韩高清在线视频 | 丝瓜视频免费看黄片| 国产精品一区二区精品视频观看| netflix在线观看网站| 人人妻人人添人人爽欧美一区卜| 可以免费在线观看a视频的电影网站| 肉色欧美久久久久久久蜜桃| 亚洲九九香蕉| 69av精品久久久久久 | av线在线观看网站| 国产亚洲精品第一综合不卡| 日韩,欧美,国产一区二区三区| 亚洲av电影在线进入| 99久久人妻综合| 亚洲av成人不卡在线观看播放网 | 亚洲精品美女久久久久99蜜臀| 国产91精品成人一区二区三区 | 国产又爽黄色视频| 国产成人精品在线电影| 欧美精品一区二区免费开放| 国产精品亚洲av一区麻豆| 丰满人妻熟妇乱又伦精品不卡| 午夜久久久在线观看| 人人澡人人妻人| 国产成人精品久久二区二区免费| 一级a爱视频在线免费观看| 妹子高潮喷水视频| 午夜影院在线不卡| 国产在线免费精品| 精品国产乱码久久久久久男人| 久久精品亚洲av国产电影网| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲第一欧美日韩一区二区三区 | 日韩视频在线欧美| 欧美日韩亚洲综合一区二区三区_| 国产一区二区三区在线臀色熟女 | 亚洲中文字幕日韩| 亚洲久久久国产精品| 日本vs欧美在线观看视频| 天堂俺去俺来也www色官网| 欧美激情 高清一区二区三区| 可以免费在线观看a视频的电影网站| 免费在线观看影片大全网站| 国产一区二区三区在线臀色熟女 | 狂野欧美激情性xxxx| 99国产精品一区二区三区| 男女床上黄色一级片免费看| 国产av又大| 99九九在线精品视频| 午夜免费观看性视频| 亚洲天堂av无毛| 亚洲欧美清纯卡通| 国产免费av片在线观看野外av| 亚洲人成电影免费在线| 亚洲精品国产av成人精品| 欧美变态另类bdsm刘玥| 亚洲一卡2卡3卡4卡5卡精品中文| 免费观看a级毛片全部| 99九九在线精品视频| 老汉色av国产亚洲站长工具| 精品熟女少妇八av免费久了| 亚洲va日本ⅴa欧美va伊人久久 | 老汉色av国产亚洲站长工具| 两人在一起打扑克的视频| 国产成人a∨麻豆精品| 国产欧美日韩一区二区三 | 欧美性长视频在线观看| 90打野战视频偷拍视频| 一级毛片女人18水好多| 在线永久观看黄色视频| bbb黄色大片| 水蜜桃什么品种好| 天堂8中文在线网| 自线自在国产av| 亚洲欧美一区二区三区黑人| 中文字幕最新亚洲高清| 亚洲 国产 在线| 免费在线观看黄色视频的| 不卡av一区二区三区| 久久中文看片网| 国产男女内射视频| 亚洲精品中文字幕在线视频| 狠狠狠狠99中文字幕| 在线观看舔阴道视频| 国产男女超爽视频在线观看| 精品国内亚洲2022精品成人 | 一进一出抽搐动态| 亚洲成人手机| 亚洲成人免费av在线播放| 亚洲色图 男人天堂 中文字幕| 97精品久久久久久久久久精品| 精品人妻熟女毛片av久久网站| 成年人黄色毛片网站| 日韩 亚洲 欧美在线| 精品人妻1区二区| 少妇裸体淫交视频免费看高清 | 五月天丁香电影| 在线永久观看黄色视频| 99久久精品国产亚洲精品| 少妇 在线观看| 亚洲激情五月婷婷啪啪| 91精品伊人久久大香线蕉| 一二三四社区在线视频社区8| 亚洲精品国产一区二区精华液| 天天影视国产精品| 18禁裸乳无遮挡动漫免费视频| 纯流量卡能插随身wifi吗| 99精品久久久久人妻精品| 久久国产精品影院| 视频在线观看一区二区三区| 无限看片的www在线观看| 黄色视频,在线免费观看| 丁香六月欧美| 亚洲男人天堂网一区| 高清在线国产一区| 国产亚洲精品第一综合不卡| 人人妻人人澡人人看| 国产精品一二三区在线看| 久久精品国产亚洲av高清一级| 老司机亚洲免费影院| 叶爱在线成人免费视频播放| 啦啦啦视频在线资源免费观看| 国产精品.久久久| 99九九在线精品视频| 欧美黄色片欧美黄色片| 日本wwww免费看| 亚洲国产欧美日韩在线播放| 国产一区二区三区av在线| 蜜桃国产av成人99| 亚洲精品久久久久久婷婷小说| 色婷婷久久久亚洲欧美| 久久狼人影院| 777米奇影视久久| 韩国精品一区二区三区| 青草久久国产| 少妇 在线观看| 亚洲激情五月婷婷啪啪| 九色亚洲精品在线播放| 男女下面插进去视频免费观看| 五月天丁香电影| 欧美精品人与动牲交sv欧美| 男人爽女人下面视频在线观看| 91av网站免费观看| 免费日韩欧美在线观看| 天天添夜夜摸| 正在播放国产对白刺激| h视频一区二区三区| 国产男人的电影天堂91| av线在线观看网站| 十八禁人妻一区二区| 国产精品一区二区精品视频观看| cao死你这个sao货| 精品第一国产精品| 久久久精品免费免费高清| 中文字幕最新亚洲高清| 高清视频免费观看一区二区| 女性生殖器流出的白浆| 一区二区三区四区激情视频| 亚洲精品美女久久久久99蜜臀| 久久人人爽av亚洲精品天堂| 亚洲精品一二三| 久久午夜综合久久蜜桃| 免费在线观看完整版高清| 精品国产乱子伦一区二区三区 | 性少妇av在线| 国产精品久久久av美女十八| 亚洲精品国产一区二区精华液| 国产主播在线观看一区二区| 搡老乐熟女国产| 亚洲,欧美精品.| 日日摸夜夜添夜夜添小说| 日韩,欧美,国产一区二区三区| 脱女人内裤的视频| 亚洲人成77777在线视频| 老司机深夜福利视频在线观看 | 久久精品亚洲熟妇少妇任你| 国产伦人伦偷精品视频| 老熟妇乱子伦视频在线观看 | 久久久久视频综合| 国产淫语在线视频| 制服人妻中文乱码| 蜜桃在线观看..| 精品一区二区三区av网在线观看 | 男人舔女人的私密视频| 日韩一卡2卡3卡4卡2021年| 成人国产一区最新在线观看| 国产精品亚洲av一区麻豆| av视频免费观看在线观看| 国产成人av激情在线播放| 国产一区二区三区在线臀色熟女 | 国产精品影院久久| 精品欧美一区二区三区在线| 久久狼人影院| 丰满迷人的少妇在线观看| 久久久国产精品麻豆| 可以免费在线观看a视频的电影网站| 久久久精品94久久精品| 国产欧美日韩一区二区三 | 亚洲精品国产av成人精品| 亚洲av国产av综合av卡| 天天躁夜夜躁狠狠躁躁| 国产激情久久老熟女| 另类亚洲欧美激情| 高潮久久久久久久久久久不卡| 国产精品国产av在线观看| 久久国产亚洲av麻豆专区| 美女脱内裤让男人舔精品视频| 亚洲精品一卡2卡三卡4卡5卡 | 夜夜骑夜夜射夜夜干| 一级,二级,三级黄色视频| 美女视频免费永久观看网站| 99久久精品国产亚洲精品| 啦啦啦免费观看视频1| 精品少妇久久久久久888优播| 99精国产麻豆久久婷婷| 一边摸一边抽搐一进一出视频|