• <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
    欧美精品一区二区大全| 一级毛片aaaaaa免费看小| 99热6这里只有精品| 亚洲av免费高清在线观看| 亚洲18禁久久av| 中文字幕人妻熟人妻熟丝袜美| 亚洲av电影在线观看一区二区三区 | 国产高清有码在线观看视频| 淫秽高清视频在线观看| 亚洲国产欧美在线一区| 久久人人爽人人爽人人片va| 日韩强制内射视频| av免费在线看不卡| 免费av不卡在线播放| 最近的中文字幕免费完整| 最新中文字幕久久久久| 日韩精品青青久久久久久| 九草在线视频观看| 乱码一卡2卡4卡精品| 中文字幕亚洲精品专区| 亚洲最大成人手机在线| 一本一本综合久久| 国产黄片视频在线免费观看| 91久久精品国产一区二区三区| 久久99精品国语久久久| 天堂av国产一区二区熟女人妻| 久久精品国产鲁丝片午夜精品| videossex国产| 久久97久久精品| 麻豆国产97在线/欧美| 97超视频在线观看视频| 久久这里只有精品中国| 神马国产精品三级电影在线观看| 欧美日韩综合久久久久久| 肉色欧美久久久久久久蜜桃 | 精品99又大又爽又粗少妇毛片| 激情 狠狠 欧美| 国产爱豆传媒在线观看| 国产精品一区二区三区四区免费观看| 直男gayav资源| 国产精品国产三级国产专区5o| 国产亚洲av嫩草精品影院| 日本一二三区视频观看| 久久久久九九精品影院| 国产av国产精品国产| 久久综合国产亚洲精品| 成人漫画全彩无遮挡| 最近视频中文字幕2019在线8| 精品久久国产蜜桃| 久久综合国产亚洲精品| 久久精品久久精品一区二区三区| 久久午夜福利片| 久久久精品欧美日韩精品| 亚洲熟女精品中文字幕| 免费观看在线日韩| 欧美激情久久久久久爽电影| 少妇的逼水好多| 日本一本二区三区精品| 亚洲精品456在线播放app| 中文字幕制服av| 亚洲国产成人一精品久久久| 国产v大片淫在线免费观看| 麻豆精品久久久久久蜜桃| 午夜激情福利司机影院| 熟女人妻精品中文字幕| 成年人午夜在线观看视频 | 亚洲精品国产av蜜桃| 欧美极品一区二区三区四区| 亚洲人成网站高清观看| 日韩一区二区视频免费看| 免费黄频网站在线观看国产| 免费观看的影片在线观看| 一区二区三区四区激情视频| 色哟哟·www| 久久久国产一区二区| av又黄又爽大尺度在线免费看| 久久久国产一区二区| 日韩欧美精品免费久久| 老女人水多毛片| 久久久久国产网址| 国产女主播在线喷水免费视频网站 | 亚洲精品第二区| 综合色丁香网| a级毛片免费高清观看在线播放| 爱豆传媒免费全集在线观看| 美女国产视频在线观看| 久久韩国三级中文字幕| 国产综合懂色| 中文天堂在线官网| 视频中文字幕在线观看| av线在线观看网站| 欧美xxⅹ黑人| 亚洲综合色惰| 日韩强制内射视频| 一级毛片电影观看| 亚洲精品日韩av片在线观看| 久久久久国产网址| 国产精品一区二区三区四区免费观看| 欧美人与善性xxx| 色综合色国产| 国产成人精品久久久久久| 全区人妻精品视频| 免费观看在线日韩| 午夜福利网站1000一区二区三区| 国产老妇女一区| 老司机影院毛片| 久久这里有精品视频免费| 亚洲av不卡在线观看| 国产av国产精品国产| 国产成年人精品一区二区| 国产真实伦视频高清在线观看| 国产男女超爽视频在线观看| 国产精品嫩草影院av在线观看| 国产女主播在线喷水免费视频网站 | 麻豆国产97在线/欧美| 国内精品美女久久久久久| 老女人水多毛片| 欧美人与善性xxx| 三级毛片av免费| av在线播放精品| 99热网站在线观看| 国产一区二区三区av在线| 精品熟女少妇av免费看| 亚洲欧美精品自产自拍| 在线天堂最新版资源| 伦理电影大哥的女人| 亚洲精华国产精华液的使用体验| 日本一本二区三区精品| 国产精品久久久久久久久免| 午夜福利成人在线免费观看| 99热这里只有是精品50| 日韩成人伦理影院| 国产成人freesex在线| 国内精品一区二区在线观看| 人人妻人人澡人人爽人人夜夜 | 亚洲国产av新网站| 人妻少妇偷人精品九色| 免费看光身美女| av国产免费在线观看| 嫩草影院新地址| 亚洲国产成人一精品久久久| 欧美高清成人免费视频www| 成人亚洲精品一区在线观看 | 欧美+日韩+精品| 国产综合懂色| 99热这里只有精品一区| 日日啪夜夜爽| 精品国产一区二区三区久久久樱花 | eeuss影院久久| 亚洲第一区二区三区不卡| 国产精品久久久久久久电影| h日本视频在线播放| 国产av码专区亚洲av| av.在线天堂| 亚洲欧美一区二区三区黑人 | 国产精品不卡视频一区二区| 日韩在线高清观看一区二区三区| 老师上课跳d突然被开到最大视频| 国产一级毛片在线| 午夜福利成人在线免费观看| av在线老鸭窝| 特级一级黄色大片| 成年版毛片免费区| 亚洲av日韩在线播放| 日韩精品青青久久久久久| 九九久久精品国产亚洲av麻豆| 国产男人的电影天堂91| 欧美xxxx黑人xx丫x性爽| 真实男女啪啪啪动态图| 两个人的视频大全免费| 日日干狠狠操夜夜爽| 中文字幕亚洲精品专区| kizo精华| 在线天堂最新版资源| 99视频精品全部免费 在线| 十八禁网站网址无遮挡 | 老司机影院成人| 国产老妇伦熟女老妇高清| 久久99热6这里只有精品| 中文字幕制服av| 在线播放无遮挡| 亚洲欧美一区二区三区黑人 | 国产男人的电影天堂91| 最近最新中文字幕免费大全7| 国产精品国产三级国产av玫瑰| 免费看美女性在线毛片视频| 一二三四中文在线观看免费高清| 麻豆成人午夜福利视频| 又爽又黄无遮挡网站| 国产免费又黄又爽又色| 亚洲熟女精品中文字幕| 91久久精品国产一区二区三区| 午夜老司机福利剧场| 久久久色成人| 干丝袜人妻中文字幕| 国产精品蜜桃在线观看| 亚洲久久久久久中文字幕| 建设人人有责人人尽责人人享有的 | 一个人看视频在线观看www免费| 国产伦精品一区二区三区四那| 99久久九九国产精品国产免费| 黄片无遮挡物在线观看| 国产亚洲午夜精品一区二区久久 | 免费高清在线观看视频在线观看| 日本午夜av视频| 搡老妇女老女人老熟妇| 国产精品久久久久久久久免| 国产免费一级a男人的天堂| 最近视频中文字幕2019在线8| 精品久久久久久成人av| 午夜日本视频在线| 高清欧美精品videossex| 免费av观看视频| 日韩不卡一区二区三区视频在线| 少妇熟女欧美另类| 日本免费在线观看一区| 久久97久久精品| 成人美女网站在线观看视频| 亚洲乱码一区二区免费版| 毛片一级片免费看久久久久| kizo精华| 国产亚洲91精品色在线| 丝瓜视频免费看黄片| 国产爱豆传媒在线观看| 一区二区三区乱码不卡18| 亚洲av中文字字幕乱码综合| 色吧在线观看| 亚洲精品456在线播放app| av.在线天堂| 青青草视频在线视频观看| 精品久久久久久久久久久久久| 久久精品人妻少妇| 久久国内精品自在自线图片| 白带黄色成豆腐渣| a级毛片免费高清观看在线播放| 亚洲美女搞黄在线观看| av播播在线观看一区| 插阴视频在线观看视频| 日韩,欧美,国产一区二区三区| 美女黄网站色视频| 久久久久久久大尺度免费视频| 麻豆久久精品国产亚洲av| 天美传媒精品一区二区| 美女被艹到高潮喷水动态| 亚洲精品成人久久久久久| 99久国产av精品国产电影| 特级一级黄色大片| 男女国产视频网站| 精品久久久噜噜| 欧美日韩精品成人综合77777| 99久久精品国产国产毛片| 日日摸夜夜添夜夜添av毛片| 2018国产大陆天天弄谢| 色尼玛亚洲综合影院| 久久久午夜欧美精品| av国产久精品久网站免费入址| 国产精品人妻久久久久久| 久久久久久伊人网av| 国产精品一二三区在线看| 亚洲国产高清在线一区二区三| 久久综合国产亚洲精品| 成人午夜精彩视频在线观看| 国产精品久久久久久精品电影| 午夜激情欧美在线| 插逼视频在线观看| 久久6这里有精品| 久久久精品免费免费高清| 99久久精品热视频| 久久久久久久大尺度免费视频| 成人毛片a级毛片在线播放| 一二三四中文在线观看免费高清| 人妻一区二区av| 亚洲精品色激情综合| 极品教师在线视频| 日本与韩国留学比较| 亚洲国产最新在线播放| 国产 一区 欧美 日韩| 高清午夜精品一区二区三区| 日本猛色少妇xxxxx猛交久久| av女优亚洲男人天堂| 中文字幕久久专区| 久久99精品国语久久久| 91在线精品国自产拍蜜月| 国产v大片淫在线免费观看| 日本熟妇午夜| 亚洲av男天堂| 国产精品美女特级片免费视频播放器| 久久这里有精品视频免费| 在线观看av片永久免费下载| 一级毛片久久久久久久久女| 色尼玛亚洲综合影院| 成人亚洲欧美一区二区av| 亚洲国产av新网站| 成人性生交大片免费视频hd| 久久精品熟女亚洲av麻豆精品 | 国产视频首页在线观看| 又爽又黄无遮挡网站| 内射极品少妇av片p| 亚洲欧美日韩无卡精品| 国产一区亚洲一区在线观看| 欧美日韩精品成人综合77777| 少妇熟女aⅴ在线视频| 日日啪夜夜撸| 国产成人精品婷婷| 高清欧美精品videossex| 又爽又黄a免费视频| 十八禁网站网址无遮挡 | 国产一区二区三区av在线| 在线免费十八禁| 最近中文字幕2019免费版| 69人妻影院| 午夜免费男女啪啪视频观看| 国产精品日韩av在线免费观看| 免费黄频网站在线观看国产| 99热网站在线观看| 精品亚洲乱码少妇综合久久| 看非洲黑人一级黄片| 日本黄色片子视频| 只有这里有精品99| 午夜激情欧美在线| 1000部很黄的大片| 毛片女人毛片| 午夜福利在线观看免费完整高清在| 欧美激情在线99| 国产精品不卡视频一区二区| 亚洲最大成人av| 欧美一级a爱片免费观看看| 九九久久精品国产亚洲av麻豆| 久久草成人影院| 亚洲精品久久久久久婷婷小说| 成年人午夜在线观看视频 | 一区二区三区四区激情视频| 十八禁网站网址无遮挡 | 精品酒店卫生间| 我要看日韩黄色一级片| 欧美变态另类bdsm刘玥| 国国产精品蜜臀av免费| 久久精品熟女亚洲av麻豆精品 | 成人亚洲精品av一区二区| 国产黄色免费在线视频| 简卡轻食公司| 国产乱人视频| 午夜激情欧美在线| 午夜福利在线在线| 久久人人爽人人片av| 国产永久视频网站| 国产人妻一区二区三区在| 欧美xxⅹ黑人| 日韩亚洲欧美综合| 寂寞人妻少妇视频99o| 欧美成人一区二区免费高清观看| 久久久久久久久久人人人人人人| 性插视频无遮挡在线免费观看| 人妻制服诱惑在线中文字幕| 欧美zozozo另类| 亚洲欧美日韩东京热| 在线观看美女被高潮喷水网站| kizo精华| 免费大片18禁| 亚洲三级黄色毛片| 国产欧美另类精品又又久久亚洲欧美| 观看免费一级毛片| 亚洲自偷自拍三级| 亚洲在线自拍视频| av又黄又爽大尺度在线免费看| 国产在视频线在精品| 精品久久久久久电影网| 亚洲精品乱码久久久v下载方式| 97人妻精品一区二区三区麻豆| 色视频www国产| 亚洲精品视频女| 狂野欧美白嫩少妇大欣赏| 亚洲欧美精品自产自拍| av女优亚洲男人天堂| 麻豆久久精品国产亚洲av| 国内精品宾馆在线| 一区二区三区高清视频在线| 在线免费十八禁| 成人高潮视频无遮挡免费网站| 中文精品一卡2卡3卡4更新| 国产综合懂色| 亚洲一级一片aⅴ在线观看| 日韩欧美国产在线观看| 婷婷六月久久综合丁香| 人妻系列 视频| 最近中文字幕2019免费版| 欧美xxxx黑人xx丫x性爽| 午夜免费男女啪啪视频观看| 国产精品国产三级国产av玫瑰| 亚洲精品乱码久久久v下载方式| 日韩不卡一区二区三区视频在线| 伊人久久国产一区二区| 99久久中文字幕三级久久日本| 午夜久久久久精精品| 久久久久久久久大av| 国产在视频线在精品| 久久久久久久久大av| 欧美成人午夜免费资源| 色综合亚洲欧美另类图片| 熟妇人妻不卡中文字幕| 久久久久久伊人网av| 中文字幕人妻熟人妻熟丝袜美| 精华霜和精华液先用哪个| 干丝袜人妻中文字幕| 国产真实伦视频高清在线观看| 又爽又黄无遮挡网站| 少妇人妻一区二区三区视频| 亚洲人成网站在线播| 又大又黄又爽视频免费| 高清在线视频一区二区三区| 我的女老师完整版在线观看| 国产免费又黄又爽又色| 久久久精品94久久精品| 亚洲av中文字字幕乱码综合| 欧美性感艳星| 国产精品一区二区三区四区免费观看| 好男人视频免费观看在线| 人体艺术视频欧美日本| 在现免费观看毛片| 国产精品综合久久久久久久免费| 久久久久性生活片| 婷婷色综合大香蕉| 九九爱精品视频在线观看| 亚洲怡红院男人天堂| 日韩欧美国产在线观看| 国产欧美另类精品又又久久亚洲欧美| 成人欧美大片| 丰满乱子伦码专区| 人妻一区二区av| 69人妻影院| 三级经典国产精品| 男插女下体视频免费在线播放| 亚洲欧美一区二区三区黑人 | 纵有疾风起免费观看全集完整版 | 九草在线视频观看| 丰满少妇做爰视频| 永久免费av网站大全| 国产麻豆成人av免费视频| 久久精品人妻少妇| 亚洲精品456在线播放app| 最新中文字幕久久久久| 在线观看人妻少妇| 一级片'在线观看视频| 精品久久久久久久人妻蜜臀av| 淫秽高清视频在线观看| 国产一区二区三区综合在线观看 | 国产国拍精品亚洲av在线观看| 天堂av国产一区二区熟女人妻| 嫩草影院精品99| 国产成人freesex在线| 91精品伊人久久大香线蕉| 九九在线视频观看精品| 简卡轻食公司| 亚洲自拍偷在线| 成人亚洲精品av一区二区| 我的老师免费观看完整版| 久久久久久伊人网av| 久久综合国产亚洲精品| 亚洲四区av| 国内精品一区二区在线观看| 狂野欧美激情性xxxx在线观看| 日韩在线高清观看一区二区三区| 国产高清有码在线观看视频| 亚洲成人av在线免费| 国产精品嫩草影院av在线观看| 久久精品人妻少妇| 91精品国产九色| 九九久久精品国产亚洲av麻豆| 国产高清三级在线| 淫秽高清视频在线观看| 免费大片黄手机在线观看| 美女内射精品一级片tv| 亚洲av日韩在线播放| 亚洲精品第二区| 日本色播在线视频| 天堂影院成人在线观看| av又黄又爽大尺度在线免费看| 国产在线男女| 少妇裸体淫交视频免费看高清| 日韩国内少妇激情av| 欧美xxxx黑人xx丫x性爽| 三级国产精品欧美在线观看| 国产精品一二三区在线看| 少妇丰满av| 国产高清三级在线| or卡值多少钱| 白带黄色成豆腐渣| 中文乱码字字幕精品一区二区三区 | 一个人免费在线观看电影| 伊人久久精品亚洲午夜| 婷婷色综合www| 久久精品人妻少妇| 午夜精品国产一区二区电影 | 日本猛色少妇xxxxx猛交久久| 国产伦理片在线播放av一区| 18禁在线无遮挡免费观看视频| 精品不卡国产一区二区三区| 国产免费视频播放在线视频 | 综合色丁香网| 可以在线观看毛片的网站| 777米奇影视久久| 乱码一卡2卡4卡精品| 久久精品夜夜夜夜夜久久蜜豆| 久热久热在线精品观看| 免费高清在线观看视频在线观看| 纵有疾风起免费观看全集完整版 | 一级黄片播放器| 亚洲精品乱码久久久久久按摩| 插阴视频在线观看视频| 亚洲色图av天堂| 日韩成人av中文字幕在线观看| 又爽又黄a免费视频| 伦精品一区二区三区| 免费大片18禁| 亚洲av中文字字幕乱码综合| 亚洲精品乱久久久久久| 成年免费大片在线观看| 大香蕉久久网| 久久综合国产亚洲精品| 青春草视频在线免费观看| 日韩三级伦理在线观看| 国产久久久一区二区三区| 婷婷色综合大香蕉| 午夜福利在线观看免费完整高清在| 99久久精品国产国产毛片| 国产有黄有色有爽视频| 久久久久九九精品影院| av.在线天堂| 99久国产av精品国产电影| 寂寞人妻少妇视频99o| 亚洲国产色片| 国产永久视频网站| 亚洲人与动物交配视频| 亚洲av国产av综合av卡| 日韩 亚洲 欧美在线| 久久精品人妻少妇| 99热6这里只有精品| 午夜日本视频在线| 国产中年淑女户外野战色| 亚洲精品,欧美精品| 尤物成人国产欧美一区二区三区| 国产精品爽爽va在线观看网站| 99久久九九国产精品国产免费| 亚洲国产色片| 久久久久精品久久久久真实原创| 美女被艹到高潮喷水动态| 69av精品久久久久久| 国产精品一区二区性色av| 如何舔出高潮| 国产黄a三级三级三级人| www.av在线官网国产| 亚洲国产色片| 亚洲欧洲国产日韩| 日韩一区二区视频免费看| 久久6这里有精品| av网站免费在线观看视频 | 国产白丝娇喘喷水9色精品| 久久久久精品性色| 亚洲欧美成人综合另类久久久| 一级毛片aaaaaa免费看小| 一二三四中文在线观看免费高清| 日韩人妻高清精品专区| 亚洲精品乱码久久久v下载方式| 久久久精品94久久精品| h日本视频在线播放| 非洲黑人性xxxx精品又粗又长| 亚洲伊人久久精品综合| 午夜免费男女啪啪视频观看| 一级毛片我不卡| 免费看a级黄色片| 亚洲人成网站高清观看| 久久精品久久精品一区二区三区| 国产高潮美女av| 欧美高清性xxxxhd video| av一本久久久久| 九九在线视频观看精品| 婷婷色av中文字幕| 2018国产大陆天天弄谢| 久久久久久久久久久免费av| 97精品久久久久久久久久精品| 亚洲国产精品专区欧美| 亚洲无线观看免费| 亚洲四区av| 国产不卡一卡二| 成人午夜精彩视频在线观看| videossex国产| 亚洲欧美精品专区久久| 久久久久久久亚洲中文字幕| 国精品久久久久久国模美| 男女啪啪激烈高潮av片| 日韩欧美国产在线观看| 午夜亚洲福利在线播放| 成人毛片a级毛片在线播放| 丝瓜视频免费看黄片| 国产精品精品国产色婷婷| 亚洲欧美一区二区三区黑人 | 在线免费观看的www视频| 天天一区二区日本电影三级| 午夜爱爱视频在线播放| 免费av毛片视频| 伦精品一区二区三区| 麻豆成人午夜福利视频| 精品久久久久久电影网| 高清视频免费观看一区二区 | 亚洲自拍偷在线| 久久精品熟女亚洲av麻豆精品 | 日韩国内少妇激情av| 水蜜桃什么品种好| 一级毛片 在线播放| 亚洲欧美日韩东京热| 在线免费观看不下载黄p国产| 日日撸夜夜添| 久久鲁丝午夜福利片| 日韩 亚洲 欧美在线| 丰满人妻一区二区三区视频av|