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

    螺桿間隙及入口流量對(duì)三螺桿擠出機(jī)組合螺桿流場(chǎng)的影響

    2017-03-07 11:36:06王天書王宇航朱向哲

    王天書 , 王宇航, 蘇 海, 朱向哲

    (1.遼寧石油化工大學(xué) 機(jī)械工程學(xué)院,遼寧 撫順 113001; 2.中國(guó)石化管道儲(chǔ)運(yùn)有限公司,江蘇 徐州 221008;.錦州美聯(lián)橋汽車部件有限公司,遼寧 錦州 121013)

    螺桿間隙及入口流量對(duì)三螺桿擠出機(jī)組合螺桿流場(chǎng)的影響

    王天書1,2, 王宇航1, 蘇 海3, 朱向哲1

    (1.遼寧石油化工大學(xué) 機(jī)械工程學(xué)院,遼寧 撫順 113001; 2.中國(guó)石化管道儲(chǔ)運(yùn)有限公司,江蘇 徐州 221008;.錦州美聯(lián)橋汽車部件有限公司,遼寧 錦州 121013)

    利用Polyflow軟件,采用三維擬瞬態(tài)有限元模型,分析了螺桿間隙及入口流量對(duì)三螺桿擠出機(jī)組合螺桿聚合物流體的復(fù)雜流動(dòng)和混合特性的影響。通過(guò)可視化云圖分析了幾何中心截面軸向速度、壓力、剪切速率分布;采用粒子示蹤技術(shù),通過(guò)統(tǒng)計(jì)后處理曲線,分析了組合螺桿在不同條件下的軸向壓力、軸向面拉伸率、混合指數(shù)概率、分離尺度、累積停留時(shí)間概率函數(shù)、停留時(shí)間分布函數(shù)。結(jié)果表明,適當(dāng)增大螺桿間隙能夠增加停留時(shí)間,提高產(chǎn)品混合質(zhì)量;入口流量較小時(shí)混合效果較好。

    三螺桿擠出機(jī); 組合螺桿; 流場(chǎng)分析; 螺桿間隙; Polyflow軟件

    三螺桿擠出機(jī)是聚合物工業(yè)近年來(lái)逐步發(fā)展起來(lái)的新型設(shè)備,能夠更好地滿足高效節(jié)能的制品需求,應(yīng)用于石油化工產(chǎn)品的后續(xù)加工。關(guān)于三螺桿擠出機(jī)的數(shù)值模擬已有較多報(bào)導(dǎo)[1-3]。三螺桿擠出機(jī)獨(dú)特的中心區(qū)及多個(gè)嚙合區(qū),使其幾何結(jié)構(gòu)更為復(fù)雜,而組合螺桿的不同組合形式更加深了這種復(fù)雜程度。螺桿是物料混合的核心元件,在混合過(guò)程中,幾何中心區(qū)的面積在螺桿旋轉(zhuǎn)過(guò)程中呈現(xiàn)“大—小—大”的動(dòng)態(tài)周期性變化,而螺桿之間的間隙對(duì)中心區(qū)面積有一定影響,尤其在中心區(qū)面積最小時(shí)刻。同時(shí),多個(gè)嚙合區(qū)的嚙合間隙、機(jī)筒壁與螺桿頂徑之間的間隙,也是物料混合中剪切作用激烈的區(qū)域,因此螺桿間隙是影響三螺桿擠出機(jī)混合的重要因素之一[4-6]。

    不同的邊界條件也會(huì)導(dǎo)致三螺桿擠出機(jī)組合螺桿熔融混合效果的差異[7-8]。本文邊界條件給定入口流量,定量研究其對(duì)物料混合的影響。通過(guò)對(duì)軸截面的速度分布、剪切速率分布等云圖以及軸向平均面拉伸率分布、軸向壓力分布、面拉伸率概率、分離尺度等混合評(píng)價(jià)參數(shù)進(jìn)行對(duì)比分析,重點(diǎn)研究了三螺桿擠出機(jī)組合螺桿間隙及入口流量對(duì)流動(dòng)混合的影響,以期為實(shí)際設(shè)計(jì)生產(chǎn)提供理論參考。

    1 有限元模型和模擬條件

    本文研究的三螺桿擠出機(jī)組合螺桿的主要幾何參數(shù)見(jiàn)表1。有限元模型的建立,采用Polyflow前處理模塊GAMBIT軟件包,其本身所具有的強(qiáng)大功能及快速更新,在目前CFD前處理軟件中穩(wěn)居上游。為更精確地觀察螺桿與機(jī)筒之間的間隙、螺桿之間的間隙及螺桿中心區(qū)導(dǎo)致的流場(chǎng)變化,采用兩層邊界層劃分,對(duì)流道中心區(qū)網(wǎng)格細(xì)劃,邊界層厚度為螺桿間隙厚度。流道有限元模型網(wǎng)格總數(shù)為56 160個(gè),節(jié)點(diǎn)總數(shù)為69 322個(gè);螺桿有限元模型的網(wǎng)格總數(shù)為32 125個(gè),節(jié)點(diǎn)總數(shù)為9 637個(gè)。采用網(wǎng)格重疊技術(shù)將流道及螺桿組合[9-11]。三螺桿擠出機(jī)流道及螺桿有限元模型如圖1所示。

    表1 三螺桿擠出機(jī)組合螺桿的主要幾何參數(shù)

    1 流道入口; 2 流道外壁; 3 流道內(nèi)壁; 4 流道出口

    (b) 螺桿

    為了研究螺桿間隙對(duì)組合螺桿流體流場(chǎng)的影響,螺桿間隙選取0.3、0.5、0.8 mm三組模型。對(duì)不同入口流量進(jìn)行對(duì)比分析時(shí),采用螺桿間隙為0.5 mm的組合螺桿,入口流量分別設(shè)定為2.5×10-6、6.0×10-6m3/s及10.0×10-6m3/s。邊界條件:機(jī)筒內(nèi)表面設(shè)定法向速度vn與切向速度vs為0,螺桿轉(zhuǎn)速為50 r/min,流動(dòng)出口設(shè)定法向力Fn與切向力Fs為0。采用黏度對(duì)溫度不敏感的聚乙烯材料,其參數(shù)如下:密度為740 kg/m3,極限黏度η∞為0 Pa·s,零剪切黏度η0為5 200 Pa·s,冪率指數(shù)n為0.5。

    2 螺桿間隙對(duì)流場(chǎng)的影響

    2.1 截面軸向速度分布

    不同螺桿間隙下組合螺桿中心截面軸向速度分布云圖如圖2所示。

    (a) 0.3 mm

    (b) 0.5 mm

    (c) 0.8 mm

    由圖2可知,在反向嚙合塊位置,均產(chǎn)生藍(lán)色回流區(qū)域;螺紋元件嚙合位置有一定的正向輸送能力;螺桿間隙不同時(shí)軸向速度分布云圖變化并不明顯,但是隨著螺桿間隙逐漸增大,軸向速度略有降低。以上現(xiàn)象說(shuō)明,隨著螺桿間隙的增加,軸向速度仍然保持原有的分布趨勢(shì),但能夠在一定程度上降低軸向速度,這對(duì)于物料混合是有利的。

    2.2 截面壓力分布

    不同螺桿間隙下組合螺桿中心軸截面的壓力分布云圖如圖3所示。

    (a) 0.3 mm

    (b) 0.5 mm

    (c) 0.8 mm

    由圖3可知,螺桿間隙不同時(shí),組合螺桿中心軸截面的壓力分布趨勢(shì)基本相同,聚合物粒子從右端入口到左端出口,經(jīng)歷壓力逐漸降低的過(guò)程。在嚙合塊位置,壓力降低的趨勢(shì)表現(xiàn)更為明顯。這是因?yàn)椋簢Ш蠅K元件屬于混合元件,有降壓的作用,而螺紋元件屬于輸送元件,有建壓作用;隨著每組螺桿間隙的減小,壓差變化幅度也依次遞減,但壓力分布仍保持原有趨勢(shì)。一般情況下,負(fù)壓差越大,越有利于聚合物粒子的充分混合,從圖3可以看出,螺桿間隙為0.3 mm的組合螺桿壓差最大。

    2.3 截面剪切速率分布

    不同螺桿間隙下組合螺桿中心軸截面的剪切速率分布云圖如圖4所示。

    (a) 0.3 mm

    (b) 0.5 mm

    (c) 0.8 mm

    從圖4可以看出,高剪切區(qū)域集中在嚙合塊嚙合的位置,嚙合塊的嚙合面積相對(duì)螺紋元件更大,對(duì)物料的剪切作用更強(qiáng),不同螺桿間隙的每組速度分布趨勢(shì)基本相同,隨著螺桿間隙的增大,嚙合位置高剪切速率區(qū)域云圖色彩逐漸變淡,剪切速率呈現(xiàn)減弱態(tài)勢(shì)。這是由于流道狹窄使螺棱間剪切作用更為強(qiáng)烈。因此,小螺桿間隙使聚合物粒子在嚙合區(qū)受到的剪切作用更為激烈。

    2.4 軸向壓力分布

    為進(jìn)一步深入研究,分析了組合螺桿軸向壓力,壓力分布曲線如圖5所示。從圖5可以看出,從入口到出口各組壓力整體呈現(xiàn)下降趨勢(shì),與圖3(截面壓力分布云圖)相吻合。在通常情況下,負(fù)壓差越大,一定程度上越能夠降低軸向速度,增加停留時(shí)間,對(duì)聚合物粒子充分混合起積極作用。螺桿間隙越大,軸向平均壓力越小,也就是說(shuō)壓力隨著螺桿間隙的減小而增大,適當(dāng)減小螺桿間隙有利于提高分散混合效果。

    圖5 不同螺桿間隙下組合螺桿軸向壓力分布曲線

    2.5 混合指數(shù)概率

    組合螺桿在不同螺桿間隙下的混合指數(shù)概率分布曲線如圖6所示。如果混合指數(shù)為0~0.5,則流動(dòng)形式為剪切流動(dòng);如果混合指數(shù)大于0.5,則流動(dòng)形式為拉伸流動(dòng),而拉伸流動(dòng)在混合中的貢獻(xiàn)大于剪切流動(dòng)。從圖6可以看出,隨著螺桿間隙的增大,混合指數(shù)大的概率也變大,表明隨著螺桿間隙增大,分散混合效果越來(lái)越好。

    圖6 組合螺桿在不同螺桿間隙下的混合指數(shù)概率分布曲線

    2.6 停留時(shí)間分布

    采用粒子示蹤技術(shù)[12],研究粒子在軌跡上的流動(dòng)規(guī)律,計(jì)算結(jié)束后在Polystat中進(jìn)行后處理設(shè)置并得到相應(yīng)的處理結(jié)果。組合螺桿累積停留時(shí)間和停留時(shí)間分布(RTD)曲線如圖7所示,為便于觀察將曲線的局部放大圖置于圖7中。從圖7可以看出,每組曲線較為接近。為了進(jìn)行定量分析,通常用ΔT表示T75-T0。式中,T0表示首個(gè)粒子流出的時(shí)間,s,T0越大對(duì)混合越有利;T75表示75%的粒子流出流場(chǎng)的時(shí)間,s,T75越小對(duì)混合越有利;ΔT越小表明混合效果越好。不同螺桿間隙各組RTD分布參數(shù)見(jiàn)表2。統(tǒng)計(jì)發(fā)現(xiàn),T0隨著螺桿間隙的增加而增大,T75隨著螺桿間隙的增加而減小,因此ΔT也隨著螺桿間隙的增加而減小,這對(duì)物料的均勻混合產(chǎn)生積極的作用。因此,適當(dāng)增大螺桿間隙能夠增加停留時(shí)間,提高產(chǎn)品混合質(zhì)量。

    (a) 累積停留時(shí)間

    (b) 停留時(shí)間

    螺桿間隙/mmT0/sT75/sΔT/s0.31.782811.32179.53890.52.928310.51667.58830.84.191510.50766.3161

    3 入口流量對(duì)流場(chǎng)的影響

    3.1 截面軸向速度分布

    不同入口流量下的組合螺桿中心軸截面的速度分布云圖如圖8所示。由圖8可知,每組在不同流量下的速度分布具有相同趨勢(shì);速度隨著入口流量的增加而增加;各組螺紋元件與嚙合塊元件相接的嚙合塊、中心區(qū)位置的速度云圖顏色也明顯地加深。這是因?yàn)椋郝菁y元件是輸送元件,有建壓能力,而嚙合塊是混合元件,能降低壓力和速度,兩者的交界位置有較大的壓差變化。從總體來(lái)看,速度隨著入口流量的增大而增大,停留時(shí)間逐漸減小,這對(duì)于物料充分混合是不利的。因而,入口流量相對(duì)較小的組合,能夠有效降低軸向速度,增加停留時(shí)間,混合能力更強(qiáng)。

    (a) 2.5×10-6 m3/s

    (b) 6.0×10-6 m3/s

    (c) 10.0×10-6 m3/s

    3.2 截面壓力分布

    不同入口流量下組合螺桿中心軸截面的壓力分布云圖如圖9所示。由圖9可知,入口流量為2.5×10-6m3/s時(shí),在入口端嚙合塊位置和出口端嚙合塊位置對(duì)稱分布局部最大壓力和最小壓力,在螺紋元件與嚙合塊元件交界位置,對(duì)稱分布局部低壓區(qū)和高壓區(qū);其余兩組壓力云圖基本呈現(xiàn)遞減趨勢(shì);隨著入口流量的增加,負(fù)壓差逐漸增大,這對(duì)于物料混合是有利的。

    (a) 2.5×10-6 m3/s

    (b) 6.0×10-6 m3/s

    (c) 10.0×10-6 m3/s

    3.3 截面剪切速率分布

    不同入口流量下組合螺桿中心軸截面剪切速率分布云圖如圖10所示。

    (a) 2.5×10-6 m3/s

    (b) 6.0×10-6 m3/s

    (c) 10.0×10-6 m3/s

    由圖10可以看出,入口流量不同時(shí),截面的剪切速率分布趨勢(shì)基本相同,這說(shuō)明改變?nèi)肟诹髁坎](méi)有改變剪切速率原有的分布趨勢(shì);隨著入口流量的增大,剪切速率略有增大。這對(duì)于聚合物粒子的分散混合有一定的增強(qiáng)作用。

    3.4 軸向平均壓力分布

    為了進(jìn)一步深入研究,分析了組合螺桿軸向平均壓力,結(jié)果如圖11所示。從圖11可以看出,當(dāng)流量為2.5×10-6m3/s時(shí),壓力有上下波動(dòng)趨勢(shì),這種波動(dòng)更直觀地反映了嚙合塊與螺紋元件的不同功能;當(dāng)流量為6.0×10-6m3/s和10.0×10-6m3/s時(shí),壓力隨軸向距離的增加而下降,隨著入口流量的增加,負(fù)壓差也隨之增大,與圖9的分析結(jié)果一致。

    圖11 不同入口流量下組合螺桿軸向壓力分布曲線

    3.5 面拉伸率

    不同入口流量下組合螺桿軸向面拉伸率分布曲線如圖12所示。

    圖12 不同入口流量下組合螺桿軸向面拉伸率分布曲線

    從圖12可以看出,沿軸向呈現(xiàn)指數(shù)增長(zhǎng)趨勢(shì),這是有效層流混合的必要條件之一,表明三螺桿擠出機(jī)組合螺桿的分散混合能力較好。隨著入口流量的依次減小,面拉伸率曲線明顯升高,并且入口流量為2.5×10-6m3/s的面拉伸率遠(yuǎn)遠(yuǎn)高于入口流量為6.0×10-6m3/s和10.0×10-6m3/s時(shí)的面拉伸率。入口流量為6.0×10-6m3/s時(shí)的面拉伸率比入口流量為10.0×10-6m3/s時(shí)高出41.4%,入口流量為2.5×10-6m3/s時(shí)的面拉伸率比入口流量為6.0×10-6m3/s時(shí)高出118.0%,這種巨大的差距表明,隨著入口流量的降低,拉伸作用增強(qiáng),并且這種差距隨著流量的減小而增大,粒子在較小入口流量中混合所受到的拉伸作用更強(qiáng),而面拉伸率是衡量粒子高效混合的重要參數(shù)之一,拉伸作用比剪切作用對(duì)擠出機(jī)的分散混合更為有利。

    3.6 分離尺度

    組合螺桿在不同入口流量下的分離尺度對(duì)比曲線如圖 13所示,為便于觀察將局部放大圖置于圖13中。

    圖13 組合螺桿在不同入口流量下的分離尺度對(duì)比曲線

    由圖13可知,各組分離尺度在開始時(shí)刻都有一段急劇下降,說(shuō)明在此階段分布混合能力較強(qiáng),隨后分離尺度曲線較為平穩(wěn),各自在一個(gè)小范圍內(nèi)波動(dòng);分離尺度隨著入口流量的增加而增加。分離尺度是對(duì)擠出機(jī)整體分布混合能力進(jìn)行評(píng)價(jià)的參數(shù),并不能從分布曲線中得到混合能力較弱的具體位置,分離尺度減小表明粒子分布的濃度差減小,分布混合能力較好。

    3.7 停留時(shí)間分布

    不同流量下累積停留時(shí)間及停留時(shí)間分布曲線如圖14所示。

    (a) 累積停留時(shí)間

    (b) 停留時(shí)間

    累積停留時(shí)間分布曲線是衡量粒子分布混合效率高低的重要參數(shù)。由圖14(a)可知,入口流量小的曲線位于流量大的曲線的右側(cè),這表明入口流量小則停留時(shí)間長(zhǎng),減小入口流量能夠增加物料的停留時(shí)間,使分布混合的效果更為理想。

    從物料混合的均一性角度考慮,最小停留時(shí)間起點(diǎn)越大越好,也就是首個(gè)粒子流出的時(shí)間越長(zhǎng)越好。從圖14(b)可以看出,入口流量越小,停留時(shí)間分布起點(diǎn)越大;同時(shí),隨著入口流量的減小,曲線的橫向?qū)挾纫惨来卧龃?,停留時(shí)間分布曲線的橫向分布越寬說(shuō)明縱向混合能力越強(qiáng),能夠獲得混合均一性更好的產(chǎn)品??傊?,入口流量越小,停留時(shí)間分布起點(diǎn)越小,橫向?qū)挾仍酱?,停留時(shí)間分布更為合理,混合效果更佳。

    4 結(jié) 論

    (1)小螺桿間隙使聚合物粒子在嚙合區(qū)的剪切作用更為激烈,但適當(dāng)增大螺桿間隙,能夠降低軸向速度,使高剪切作用頻率降低,并且螺桿間隙增大能夠增強(qiáng)拉伸流動(dòng),有利于粒子充分混合;小螺桿間隙能夠降低分離尺度,但通過(guò)停留時(shí)間分布曲線發(fā)現(xiàn),螺桿間隙大的組合螺桿ΔT值較小。因此,從總體來(lái)看適當(dāng)增大螺桿間隙能夠增加停留時(shí)間,提高產(chǎn)品的混合質(zhì)量。

    (2)增大入口流量,剪切速率略有增大,但入口流量相對(duì)較小的組合螺桿,能夠有效降低軸向速度,即縮短高剪切作用的時(shí)間;隨著入口流量的降低,拉伸作用更為頻繁,粒子在較小入口流量中混合所受到的拉伸作用更強(qiáng);盡管入口流量增大能夠提高負(fù)壓差,但隨著入口流量的降低,分離尺度降低,分布混合能力增強(qiáng)。同時(shí),入口流量越小,停留時(shí)間分布更為合理,混合效果更佳。因此,從總體來(lái)看較小的入口流量混合效果更為理想。

    [1] 朱向哲, 謝禹鈞, 苗一.雙螺桿反應(yīng)器同向嚙合螺紋元件的三維流場(chǎng)分析[J].石油化工高等學(xué)校學(xué)報(bào),2004,17(3):66-69.

    [2] Zhu X Z,Xie Y J,Yuan H Q.Numerical simulation of extrusion characteristics for co-rotating tri-screw extruder[J]. Polym Plastics Technol.Eng.,2007,46:401-407.

    [3] Zhu X Z,Yuan H Q,Wang W Q.Numerical simulation of flow characteristics in new co-rotating triangle arrayed triple screw extruders[J].Journal of Materials Processing Technology,2009,209(7):3289-3299.

    [4] 龐軍艦,何亞?wèn)|,信春玲,等.三螺桿擠出機(jī)螺桿元件混合特性的數(shù)值研究[J].中國(guó)塑料,2012,26(1):103-108.

    [5] Wang G,Zhu X Z,He Y D,et al.Effects on screw clearances and blends ratios on the mixing characteristics of tri-screw extruders in cross section with FEM[J].Engineering Applications of Computational Fluid Mechanics,2013,7(1):74-89.

    [6] 郝永輝, 陳立, 何延?xùn)|,等. 螺桿構(gòu)型對(duì)三螺桿擠出機(jī)混合效率的影響[J]. 遼寧石油化工大學(xué)學(xué)報(bào), 2015,35(3):47-50.

    [7] 苗立榮,潘龍,薛平,等.單螺桿擠出機(jī)螺旋溝槽固體段壓力的實(shí)驗(yàn)研究[J].當(dāng)代化工,2011,40(4):347-350.

    [8] 張迅,朱向哲,王閣.滑移邊界條件對(duì)同向嚙合三螺桿擠出機(jī)混合性能的影響[J].石油化工高等學(xué)校學(xué)報(bào),2012, 25(4):68-72.

    [9] Vergnes B,Valle G D,Delamare L.A global computer software for polymer flows in co-rotating twin screw extruders[J].Polymer Engineering & Science,1998,38(11):1781-1792.

    [10] Zhang X M,F(xiàn)eng L F,Hoppe S,et al.Local residence time, residence revolution,and residence volume distributions in twin-screw extruders[J].Polymer Engineering & Science,2008,48(1):19-28.

    [11] Zhang X M,F(xiàn)eng L F,Chen W X, et al.Numerical simulation and experimental validation of mixing performance of kneading discs in a twin screw extruder[J].Polymer Engineering & Science,2009,49(9):1772-1783.

    [12] Connelly R K,Kokini J L.The effect of shear thinning and differential viscoelasticity on mixing in a model 2D mixer as determined using FEM with particle tracking[J].Journal of Non-Newtonian Fluid Mechanics,2004,123(1):1-17.

    (編輯 宋錦玉)

    Effects of Different Clearance and Inlet Flow Rate on the Mixing Characteristics of the Fluid Flow Field of Combined Screws with Tri-Screw Extruder

    Wang Tianshu1,2, Wang Yuhang1, Su Hai3, Zhu Xiangzhe1

    (1.SchoolofMechanicalEngineering,LiaoningShihuaUniversity,F(xiàn)ushunLiaoning113001,China; 2.PipelineTransportationandStorageCompany,Sinopec,XuzhouJiangsu221008,China; 3.JinzhouUBridgeAutoPartsCo.,Ltd.,JinzhouLiaoning121013,China)

    The three-dimensional (3D) transient finite element model was established by using Polyflow software. The effects of screw clearance and inlet flow rate of combined screws on the complex flow and mixing characteristics of polymer melt were discussed. The axial velocity, pressure and shear rate distribution of the geometric center region were analyzed through the visual nephograms. Based on the particle tracing method, the post-processing software of statistic analysis was used to analyze the axial pressure, axial area stretching rate, the probability of mixing index, segregation scale, cumulative residence time distribution and residence time distribution in different conditions. The results showed that the increase of the screw clearance could increase the residence time distribution and improve the mixing quality of the product, and the smaller inlet flow mixing effect was more ideal.

    Tri-screw extruders; Combined screws; Flow field analysis; Screw clearance; Polyflow

    1672-6952(2017)01-0046-07

    投稿網(wǎng)址:http://journal.lnpu.edu.cn

    2016-05-05

    2016-06-30

    國(guó)家自然科學(xué)基金項(xiàng)目(50903042,51473073)。

    王天書(1988-),女,碩士研究生,從事高效節(jié)能石化裝備方面的研究;E-mail:812357296@qq.com。

    朱向哲(1974-),男,博士,教授,從事高效節(jié)能石化裝備方面的研究;E-mail:xzzhu@126.com。

    TQ320.5+1

    A

    10.3969/j.issn.1672-6952.2017.01.010

    久久性视频一级片| 非洲黑人性xxxx精品又粗又长| 日本一区二区免费在线视频| 成人三级做爰电影| 久久精品亚洲熟妇少妇任你| 亚洲成人精品中文字幕电影| 极品教师在线免费播放| 国产区一区二久久| 久久青草综合色| 天天一区二区日本电影三级 | 亚洲精品av麻豆狂野| 欧美乱码精品一区二区三区| 人人妻,人人澡人人爽秒播| 国产黄a三级三级三级人| 色精品久久人妻99蜜桃| 亚洲国产中文字幕在线视频| 亚洲va日本ⅴa欧美va伊人久久| 亚洲情色 制服丝袜| 中文字幕精品免费在线观看视频| 两性夫妻黄色片| 这个男人来自地球电影免费观看| 久久久精品欧美日韩精品| 天天添夜夜摸| 精品久久久久久久毛片微露脸| 免费av毛片视频| 一二三四社区在线视频社区8| 免费高清在线观看日韩| www.999成人在线观看| 一进一出好大好爽视频| www.999成人在线观看| 宅男免费午夜| 午夜a级毛片| 国产在线观看jvid| 91成人精品电影| 国产亚洲精品第一综合不卡| 日韩精品青青久久久久久| 日韩精品青青久久久久久| 久久亚洲精品不卡| 纯流量卡能插随身wifi吗| 欧美一区二区精品小视频在线| 亚洲欧美激情在线| 女性被躁到高潮视频| 国产亚洲精品综合一区在线观看 | 精品国产一区二区三区四区第35| 一级a爱片免费观看的视频| 精品乱码久久久久久99久播| 日韩欧美三级三区| 欧美激情极品国产一区二区三区| 久久国产精品男人的天堂亚洲| 黄色毛片三级朝国网站| 久久久久精品国产欧美久久久| 久久久久久久久免费视频了| 亚洲国产日韩欧美精品在线观看 | 日韩 欧美 亚洲 中文字幕| 免费在线观看影片大全网站| 侵犯人妻中文字幕一二三四区| 亚洲一区二区三区色噜噜| 精品人妻1区二区| 欧美激情 高清一区二区三区| 国产伦人伦偷精品视频| 日韩av在线大香蕉| 三级毛片av免费| 色综合欧美亚洲国产小说| 国产日韩一区二区三区精品不卡| 欧美成人免费av一区二区三区| 国内精品久久久久久久电影| 久久久久久久精品吃奶| 少妇裸体淫交视频免费看高清 | 91av网站免费观看| 在线免费观看的www视频| 国产精品日韩av在线免费观看 | 亚洲人成伊人成综合网2020| 丰满人妻熟妇乱又伦精品不卡| 久久青草综合色| 两性午夜刺激爽爽歪歪视频在线观看 | 国内毛片毛片毛片毛片毛片| 淫秽高清视频在线观看| 国产精品爽爽va在线观看网站 | 波多野结衣高清无吗| 欧美激情 高清一区二区三区| 日韩av在线大香蕉| 亚洲精品国产一区二区精华液| 亚洲自偷自拍图片 自拍| 天天躁夜夜躁狠狠躁躁| 国产精品免费一区二区三区在线| 国产精品二区激情视频| 欧美黄色片欧美黄色片| 美女高潮到喷水免费观看| 久久天躁狠狠躁夜夜2o2o| 国产高清videossex| 黑人巨大精品欧美一区二区mp4| 1024香蕉在线观看| 一进一出抽搐动态| 欧美乱色亚洲激情| 国产精品久久久人人做人人爽| 欧美成人一区二区免费高清观看 | 成在线人永久免费视频| 午夜久久久在线观看| 少妇 在线观看| 国产精品,欧美在线| 久久九九热精品免费| 脱女人内裤的视频| 亚洲国产精品999在线| 最近最新中文字幕大全免费视频| 精品久久久久久久久久免费视频| 神马国产精品三级电影在线观看 | 亚洲av日韩精品久久久久久密| 极品人妻少妇av视频| 精品久久久精品久久久| 精品国产美女av久久久久小说| 亚洲人成电影观看| 很黄的视频免费| 亚洲国产欧美网| 亚洲伊人色综图| 叶爱在线成人免费视频播放| 色尼玛亚洲综合影院| 一二三四在线观看免费中文在| 1024香蕉在线观看| 日本vs欧美在线观看视频| 一本大道久久a久久精品| 一区福利在线观看| 69精品国产乱码久久久| 两性午夜刺激爽爽歪歪视频在线观看 | 又黄又粗又硬又大视频| 国产不卡一卡二| 丝袜人妻中文字幕| 99久久99久久久精品蜜桃| 国产精品久久久av美女十八| 黄片大片在线免费观看| 亚洲在线自拍视频| 欧美av亚洲av综合av国产av| 黄片大片在线免费观看| 女人爽到高潮嗷嗷叫在线视频| 99re在线观看精品视频| 久9热在线精品视频| 91在线观看av| 宅男免费午夜| 在线av久久热| 久久久久久免费高清国产稀缺| 女生性感内裤真人,穿戴方法视频| 国内精品久久久久久久电影| 久热爱精品视频在线9| 国产成年人精品一区二区| 操出白浆在线播放| 国语自产精品视频在线第100页| 亚洲第一青青草原| 十分钟在线观看高清视频www| 757午夜福利合集在线观看| 亚洲视频免费观看视频| 国产欧美日韩一区二区精品| 不卡av一区二区三区| 色综合亚洲欧美另类图片| 制服诱惑二区| 9191精品国产免费久久| 日韩 欧美 亚洲 中文字幕| 午夜免费激情av| 丁香六月欧美| 亚洲片人在线观看| 国产极品粉嫩免费观看在线| 亚洲最大成人中文| 变态另类丝袜制服| 亚洲男人的天堂狠狠| 亚洲电影在线观看av| 91成年电影在线观看| 成熟少妇高潮喷水视频| 18禁观看日本| 啦啦啦 在线观看视频| 成年人黄色毛片网站| 一级毛片女人18水好多| 亚洲电影在线观看av| 亚洲中文av在线| 少妇的丰满在线观看| 亚洲五月色婷婷综合| 啦啦啦 在线观看视频| 中出人妻视频一区二区| 女人被躁到高潮嗷嗷叫费观| 欧美在线一区亚洲| 日日摸夜夜添夜夜添小说| 欧美在线黄色| 一进一出抽搐动态| 精品福利观看| 国产精品精品国产色婷婷| 国产精品永久免费网站| 午夜免费成人在线视频| 香蕉国产在线看| 麻豆一二三区av精品| 美女午夜性视频免费| 波多野结衣一区麻豆| 国产伦一二天堂av在线观看| 国产精品,欧美在线| 极品人妻少妇av视频| av免费在线观看网站| av有码第一页| av天堂久久9| 日韩欧美国产一区二区入口| 亚洲五月婷婷丁香| 少妇熟女aⅴ在线视频| 人人妻人人爽人人添夜夜欢视频| 最新在线观看一区二区三区| 免费在线观看影片大全网站| 久久精品91无色码中文字幕| 女人爽到高潮嗷嗷叫在线视频| 一进一出好大好爽视频| 99久久综合精品五月天人人| 久久久久久免费高清国产稀缺| 欧美激情高清一区二区三区| 黑人欧美特级aaaaaa片| 一a级毛片在线观看| 欧美日韩一级在线毛片| 在线观看免费视频日本深夜| 高清黄色对白视频在线免费看| 我的亚洲天堂| 精品人妻1区二区| 国产伦人伦偷精品视频| 国产欧美日韩一区二区精品| 黄频高清免费视频| 日本撒尿小便嘘嘘汇集6| 久久国产亚洲av麻豆专区| 身体一侧抽搐| 99热只有精品国产| 嫁个100分男人电影在线观看| 久久人人爽av亚洲精品天堂| 精品不卡国产一区二区三区| 国产精品久久久人人做人人爽| 又大又爽又粗| 国产aⅴ精品一区二区三区波| 久久欧美精品欧美久久欧美| 99在线视频只有这里精品首页| 亚洲一区中文字幕在线| 极品教师在线免费播放| 黄色a级毛片大全视频| 国产在线观看jvid| 亚洲成a人片在线一区二区| 在线观看www视频免费| 免费观看人在逋| 高潮久久久久久久久久久不卡| 丁香欧美五月| e午夜精品久久久久久久| 妹子高潮喷水视频| 免费看十八禁软件| 欧美日韩瑟瑟在线播放| 日韩视频一区二区在线观看| 脱女人内裤的视频| 亚洲,欧美精品.| 精品久久蜜臀av无| 亚洲伊人色综图| 青草久久国产| 久热爱精品视频在线9| 免费看美女性在线毛片视频| 中文字幕久久专区| 日韩欧美三级三区| 伦理电影免费视频| 国产主播在线观看一区二区| 欧美成人免费av一区二区三区| 91字幕亚洲| 视频在线观看一区二区三区| 女性被躁到高潮视频| 亚洲无线在线观看| 国产精品电影一区二区三区| 丰满的人妻完整版| 一卡2卡三卡四卡精品乱码亚洲| 一进一出抽搐gif免费好疼| 母亲3免费完整高清在线观看| 亚洲色图综合在线观看| 岛国视频午夜一区免费看| 黄片大片在线免费观看| 9热在线视频观看99| 国产一区二区三区在线臀色熟女| 亚洲色图综合在线观看| 国产欧美日韩精品亚洲av| 法律面前人人平等表现在哪些方面| 韩国精品一区二区三区| 99精品久久久久人妻精品| 国产亚洲av高清不卡| 18禁美女被吸乳视频| 99国产精品99久久久久| 啦啦啦 在线观看视频| av天堂久久9| 精品国产一区二区三区四区第35| 久久久国产成人精品二区| 天天躁夜夜躁狠狠躁躁| 久9热在线精品视频| 亚洲欧美日韩另类电影网站| 日本免费a在线| 亚洲精品av麻豆狂野| 亚洲第一青青草原| 国产97色在线日韩免费| 在线观看66精品国产| 色综合站精品国产| 国产高清激情床上av| 久久精品成人免费网站| 看黄色毛片网站| 91字幕亚洲| 亚洲精品久久国产高清桃花| 亚洲av成人一区二区三| 精品福利观看| www.999成人在线观看| 中文字幕色久视频| 欧美色欧美亚洲另类二区 | 天天躁夜夜躁狠狠躁躁| 99国产精品一区二区蜜桃av| 在线av久久热| 在线十欧美十亚洲十日本专区| 黄片小视频在线播放| 免费不卡黄色视频| 老熟妇乱子伦视频在线观看| 国产精品二区激情视频| 国产区一区二久久| 成人国产综合亚洲| 在线永久观看黄色视频| 黄片大片在线免费观看| 两个人看的免费小视频| 老熟妇仑乱视频hdxx| 最近最新中文字幕大全免费视频| 精品一区二区三区av网在线观看| 国产野战对白在线观看| 亚洲国产精品成人综合色| 怎么达到女性高潮| 亚洲伊人色综图| 最近最新免费中文字幕在线| 日日夜夜操网爽| 91成人精品电影| 亚洲av成人av| 成人三级做爰电影| 丝袜在线中文字幕| 露出奶头的视频| 精品日产1卡2卡| 欧美日韩黄片免| 日本免费a在线| 久久国产精品人妻蜜桃| 亚洲伊人色综图| 中国美女看黄片| 国产亚洲精品综合一区在线观看 | 制服人妻中文乱码| 黄色 视频免费看| 90打野战视频偷拍视频| 色av中文字幕| 少妇 在线观看| 精品久久久久久久毛片微露脸| 日本在线视频免费播放| 97碰自拍视频| 麻豆av在线久日| 久久 成人 亚洲| 99国产极品粉嫩在线观看| 母亲3免费完整高清在线观看| 国产99久久九九免费精品| 亚洲第一av免费看| 别揉我奶头~嗯~啊~动态视频| 美女国产高潮福利片在线看| 日本精品一区二区三区蜜桃| 午夜免费鲁丝| АⅤ资源中文在线天堂| 亚洲精品av麻豆狂野| 国产黄a三级三级三级人| 免费在线观看日本一区| 国产亚洲欧美在线一区二区| 成人精品一区二区免费| av在线天堂中文字幕| 在线观看舔阴道视频| 国产精品自产拍在线观看55亚洲| 无限看片的www在线观看| 91av网站免费观看| 欧美日韩瑟瑟在线播放| 亚洲精华国产精华精| 老熟妇乱子伦视频在线观看| 欧美成人一区二区免费高清观看 | 香蕉国产在线看| 999久久久精品免费观看国产| a级毛片在线看网站| 国产1区2区3区精品| 在线十欧美十亚洲十日本专区| 欧美不卡视频在线免费观看 | 久久久久精品国产欧美久久久| 人人妻,人人澡人人爽秒播| www.www免费av| 熟妇人妻久久中文字幕3abv| 亚洲,欧美精品.| 免费久久久久久久精品成人欧美视频| www.自偷自拍.com| 中国美女看黄片| 亚洲欧美日韩另类电影网站| 欧美成狂野欧美在线观看| 老汉色∧v一级毛片| 久久香蕉激情| www.自偷自拍.com| 性欧美人与动物交配| 丁香欧美五月| 欧美国产日韩亚洲一区| 欧美黄色片欧美黄色片| 黄色毛片三级朝国网站| 亚洲无线在线观看| 又黄又粗又硬又大视频| 亚洲国产欧美一区二区综合| 亚洲人成网站在线播放欧美日韩| 99在线视频只有这里精品首页| 欧美一区二区精品小视频在线| 国产国语露脸激情在线看| 伊人久久大香线蕉亚洲五| 欧美色欧美亚洲另类二区 | 色尼玛亚洲综合影院| 他把我摸到了高潮在线观看| 日韩大尺度精品在线看网址 | 757午夜福利合集在线观看| 欧美激情高清一区二区三区| 老司机在亚洲福利影院| 老鸭窝网址在线观看| 国产视频一区二区在线看| 久久性视频一级片| ponron亚洲| 非洲黑人性xxxx精品又粗又长| 国产三级在线视频| 成熟少妇高潮喷水视频| 欧美老熟妇乱子伦牲交| 脱女人内裤的视频| av在线天堂中文字幕| 99国产精品一区二区三区| 亚洲黑人精品在线| 在线观看免费视频日本深夜| 在线观看日韩欧美| 亚洲成人精品中文字幕电影| 亚洲国产中文字幕在线视频| 国产成年人精品一区二区| 日日爽夜夜爽网站| 欧美日韩乱码在线| 国产精品免费视频内射| 香蕉国产在线看| 亚洲全国av大片| 久久久久久亚洲精品国产蜜桃av| 国产又色又爽无遮挡免费看| 国产日韩一区二区三区精品不卡| 国产精品日韩av在线免费观看 | 日韩精品免费视频一区二区三区| 国产激情久久老熟女| 丰满的人妻完整版| 亚洲色图综合在线观看| 免费观看人在逋| 午夜激情av网站| 午夜免费激情av| 国产精品综合久久久久久久免费 | 成人av一区二区三区在线看| 女人精品久久久久毛片| 精品人妻在线不人妻| 黄片播放在线免费| 日韩视频一区二区在线观看| 女性被躁到高潮视频| 19禁男女啪啪无遮挡网站| 国产1区2区3区精品| 国内毛片毛片毛片毛片毛片| 国产精品免费视频内射| 99精品久久久久人妻精品| 亚洲av成人不卡在线观看播放网| 成人三级做爰电影| 99精品欧美一区二区三区四区| 一级毛片高清免费大全| 亚洲一码二码三码区别大吗| 午夜两性在线视频| 久久这里只有精品19| 国产一区在线观看成人免费| 两个人看的免费小视频| 亚洲狠狠婷婷综合久久图片| 老司机在亚洲福利影院| 欧美人与性动交α欧美精品济南到| 神马国产精品三级电影在线观看 | 女人精品久久久久毛片| av视频免费观看在线观看| 夜夜躁狠狠躁天天躁| 亚洲狠狠婷婷综合久久图片| 精品福利观看| 中文字幕人妻熟女乱码| 免费少妇av软件| 黄色丝袜av网址大全| 啦啦啦 在线观看视频| 一本大道久久a久久精品| 久久天躁狠狠躁夜夜2o2o| 每晚都被弄得嗷嗷叫到高潮| 国产97色在线日韩免费| 1024香蕉在线观看| 亚洲av成人一区二区三| 最近最新免费中文字幕在线| 999久久久国产精品视频| 老司机午夜福利在线观看视频| 在线观看一区二区三区| 免费观看精品视频网站| 国产伦人伦偷精品视频| 亚洲电影在线观看av| 国产精品秋霞免费鲁丝片| 久久国产精品男人的天堂亚洲| 成人亚洲精品av一区二区| 日韩 欧美 亚洲 中文字幕| 99久久综合精品五月天人人| 免费在线观看视频国产中文字幕亚洲| 国产亚洲精品第一综合不卡| 亚洲人成电影免费在线| 亚洲欧美激情在线| 久久精品成人免费网站| 精品福利观看| 一级片免费观看大全| 日韩精品中文字幕看吧| 日日干狠狠操夜夜爽| 搡老熟女国产l中国老女人| 国产精品精品国产色婷婷| 精品国产国语对白av| 美女高潮到喷水免费观看| 国产精品电影一区二区三区| 亚洲熟妇熟女久久| 淫妇啪啪啪对白视频| 国产精品一区二区三区四区久久 | 男男h啪啪无遮挡| av在线天堂中文字幕| 亚洲va日本ⅴa欧美va伊人久久| 欧美日韩乱码在线| 久久国产精品人妻蜜桃| 一二三四社区在线视频社区8| 99精品在免费线老司机午夜| 国产精品日韩av在线免费观看 | 亚洲九九香蕉| 国产乱人伦免费视频| 一级作爱视频免费观看| 啦啦啦免费观看视频1| 黄片播放在线免费| 亚洲av第一区精品v没综合| 不卡一级毛片| 亚洲第一欧美日韩一区二区三区| 欧美老熟妇乱子伦牲交| 国产伦一二天堂av在线观看| 成人手机av| 亚洲国产精品sss在线观看| 成人18禁高潮啪啪吃奶动态图| 国产精品精品国产色婷婷| 亚洲欧美激情在线| 欧美黄色淫秽网站| 91精品三级在线观看| 久久精品亚洲熟妇少妇任你| 波多野结衣一区麻豆| 美女午夜性视频免费| 亚洲色图综合在线观看| 久久亚洲真实| 九色亚洲精品在线播放| 国产精品久久视频播放| 久久影院123| 天天一区二区日本电影三级 | 美女扒开内裤让男人捅视频| 欧美在线黄色| 国产精品久久电影中文字幕| 亚洲片人在线观看| 法律面前人人平等表现在哪些方面| 久热这里只有精品99| 亚洲国产中文字幕在线视频| 精品国产美女av久久久久小说| 精品久久久久久久久久免费视频| 搡老妇女老女人老熟妇| 黄色 视频免费看| 欧美日韩亚洲国产一区二区在线观看| 亚洲成人国产一区在线观看| 欧美一级a爱片免费观看看 | 麻豆av在线久日| 精品国产亚洲在线| 欧美日韩福利视频一区二区| 亚洲伊人色综图| 欧美一级毛片孕妇| 欧美在线黄色| www.熟女人妻精品国产| 亚洲av美国av| 可以在线观看毛片的网站| 亚洲人成77777在线视频| 色婷婷久久久亚洲欧美| 高清毛片免费观看视频网站| 午夜福利高清视频| 狂野欧美激情性xxxx| 巨乳人妻的诱惑在线观看| 国产av一区在线观看免费| 亚洲av成人av| 日日摸夜夜添夜夜添小说| videosex国产| 亚洲国产欧美网| 欧美成人性av电影在线观看| 变态另类成人亚洲欧美熟女 | 欧美精品亚洲一区二区| 中国美女看黄片| 亚洲中文字幕一区二区三区有码在线看 | 一级片免费观看大全| 精品一区二区三区视频在线观看免费| 欧美人与性动交α欧美精品济南到| 亚洲欧美日韩另类电影网站| 免费看十八禁软件| 日韩成人在线观看一区二区三区| 日韩大尺度精品在线看网址 | 午夜福利高清视频| 日韩欧美国产在线观看| 成在线人永久免费视频| 欧美亚洲日本最大视频资源| 欧美中文综合在线视频| 欧美日韩亚洲综合一区二区三区_| 黄色 视频免费看| 丰满人妻熟妇乱又伦精品不卡| 中文字幕另类日韩欧美亚洲嫩草| 亚洲 欧美 日韩 在线 免费| 性欧美人与动物交配| 国产一区二区三区综合在线观看| 一级a爱视频在线免费观看| 97超级碰碰碰精品色视频在线观看| 少妇裸体淫交视频免费看高清 | xxx96com| 亚洲国产精品成人综合色| 国产精品影院久久| 亚洲欧美日韩高清在线视频| 一级片免费观看大全| 亚洲 国产 在线| 亚洲第一电影网av| 久久久久久久午夜电影| 搞女人的毛片| 色精品久久人妻99蜜桃| 国产成人欧美在线观看| 日韩 欧美 亚洲 中文字幕| 久久人人爽av亚洲精品天堂|