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

    基于CFD的H型真空手套箱氣固兩相流動特性分析

    2021-11-04 01:32:30蒲志新賈加亮武志龍翟孖廣
    中國粉體技術(shù) 2021年6期
    關(guān)鍵詞:出氣口進氣口箱體

    蒲志新,賈加亮,武志龍,翟孖廣

    (遼寧工程技術(shù)大學 機械工程學院,遼寧 阜新 123000)

    進入21世紀以來,我國制造業(yè)得到了快速發(fā)展,隨之對機械加工制造工藝和加工環(huán)境提出了更高的要求,真空手套箱作為一種重要的加工平臺在越來越多的領(lǐng)域得以應(yīng)用。真空手套箱的工作原理是,利用其中的惰性氣體將箱體內(nèi)的氧氣、 水蒸氣等排出或控制箱內(nèi)氣體成分和比例。真空手套箱廣泛應(yīng)用于鋰離子電池、 半導體等領(lǐng)域。除此之外在生物學領(lǐng)域也有應(yīng)用,例如厭氧菌培養(yǎng)、 細胞低氧培養(yǎng)等[1]。

    在實際生產(chǎn)過程中,手套箱的結(jié)構(gòu)會根據(jù)生產(chǎn)環(huán)境、可利用空間、箱體內(nèi)放置設(shè)備形狀等因素,做出相應(yīng)的改變。目前,國內(nèi)手套箱生產(chǎn)企業(yè)多以仿制為主,對手套箱結(jié)構(gòu)變化之后的對應(yīng)流場狀態(tài)的改變研究很少,缺乏相應(yīng)的理論支撐,導致箱體內(nèi)產(chǎn)生理想氛圍時間較長,或者不能夠達到理想氛圍要求。本文中主要對一種結(jié)構(gòu)外形類似大寫字母H,同時滿足增材制造和減材制造的手套箱內(nèi)流域進行流固耦合仿真分析,手套箱內(nèi)增材制造過程主要運用選擇性激光燒結(jié)(selective laser sintering,SLS)[2]打印工藝,該過程主要運用金屬、陶瓷等高分子材料。在打印時,材料并非是絲線狀,而是聚合物粉末,使得打印過程會產(chǎn)生“煙霧”[3]。為了優(yōu)化打印區(qū)域氣體氛圍質(zhì)量,減少煙霧量,提高產(chǎn)品質(zhì)量,解決企業(yè)在生產(chǎn)過程中遇到的問題,本文中主要研究氣固兩相流動特性,優(yōu)化箱體加工環(huán)境,降低企業(yè)人力物力消耗,為現(xiàn)有特種工況手套箱的設(shè)計提供參考。

    1 物理建模

    傳統(tǒng)手套箱系統(tǒng)是由凈化器、手套箱以及箱體內(nèi)設(shè)備共同組成,需求不同,手套箱系統(tǒng)的結(jié)構(gòu)往往也不盡相同[4]。本文中的箱體結(jié)構(gòu)是根據(jù)客戶需求設(shè)計的,該箱體左側(cè)腔室放置一臺3D打印機,右側(cè)腔室放置一臺重型壓機,工件可以在箱體中部傳送、檢驗,該H型手套箱系統(tǒng)三維模型結(jié)構(gòu)圖如圖1所示;二維主視圖和右視圖如圖2所示。

    1—激光增材制造設(shè)備; 2—手套箱孔圈; 3—大過渡艙; 4—進氣孔; 5—真空泵; 6—出氣孔; 7—減材制造設(shè)備; 8—小過渡艙; 9—凈化器。圖1 H型手套箱系統(tǒng)三維模型結(jié)構(gòu)圖Fig.1 Standard view of 3D model of H-type glove box system

    a)主視圖b)右視圖圖2 H型真空手套箱二維視圖(單位:mm)Fig.2 Two dimensional view of H-type vacuum glove box

    2 網(wǎng)格模型

    手套箱內(nèi)部結(jié)構(gòu)較為復雜,包括置物盤、 手套、 增材制造、 減材制造設(shè)備、 大小過渡倉等,在流體分析時,如果考慮這些因素則流道抽取比較困難,個別結(jié)構(gòu)對流場分析影響比較小,可以對模型進行簡化,忽略次要因素[5-6],圖3為模型簡化后的流場結(jié)構(gòu)圖。

    圖3 模型簡化后流場結(jié)構(gòu)圖Fig.3 Flow field structure after model simplification

    為了分析網(wǎng)格質(zhì)量對數(shù)值模擬結(jié)果的影響,最大網(wǎng)格尺寸參數(shù)分別設(shè)置為100、 70、 50、 30 mm,劃分好網(wǎng)格后,進氣口速度恒定為4.777 m/s,分別求出不同網(wǎng)格數(shù)量對流道對稱面直線ab上速度值的影響,見圖4。從圖中可以看出,在最大網(wǎng)格尺寸小于70 mm時,數(shù)值模擬結(jié)果較為吻合。為了降低計算機的運算量,將網(wǎng)格劃分尺寸設(shè)置為70 mm,并將局部網(wǎng)格加密,壁面增長率為1.2[7],采用的網(wǎng)格數(shù)為270 766萬左右,最小體積為111.154 2 cm3,最大體積為297.622 1 cm3,最差的網(wǎng)格質(zhì)量為0.992。網(wǎng)格劃分示意圖如圖5所示。

    圖4 不同網(wǎng)格數(shù)量對直線ab上速度值的影響Fig.4 Influence of different grid number on simulation results

    圖5 網(wǎng)格劃分示意圖Fig.5 Grid generation diagram

    3 數(shù)學模型

    3.1 控制方程的選擇

    3.1.1 基本假設(shè)

    在建立數(shù)學模型時,忽略影響較小的因素,做如下的假設(shè): 1)氣相、固相均為連續(xù)性流體且在空間上共存; 2)箱體內(nèi)兩相流動均為湍流且固相體積濃度均勻分布于工作臺表面; 3)固體顆粒為剛性球體且表面光滑; 4)僅考慮固體顆粒的阻力和重力并忽略其他力[8-9]。

    3.1.2 基本方程

    在研究過程中,忽略吸熱、放熱等能量交換,應(yīng)用的基本方程為續(xù)性控制方程和動量控制方程[10]。

    連續(xù)性方程為

    (1)

    式中:ρ為流體密度,kg/m3;t為時間,s;u、v、w分別為x、y、z方向的速度分量,m/s。

    動量方程為

    (2)

    (3)

    (4)

    式中:fx、fy、fz分別為x、y、z方向的作用力,該作用力包含體積力和表面力。

    5個模型:標準k-ε、重整化群(RNG)k-ε、可實現(xiàn)k-ε、標準k-ω和湍流模型剪切應(yīng)力輸運(SST)k-ω,是使用最廣泛的湍流模型[11]。標準k-ε模型是Launder and Spalding提出的適用范圍廣、經(jīng)濟、合理的精度模型,它是個半經(jīng)驗的公式,是從實驗現(xiàn)象中總結(jié)出來的。RNGk-ε模型來源于嚴格的統(tǒng)計技術(shù),它和標準k-ε模型很相似,但是有以下改進: 1)RNG模型在ε方程中加了一個條件,有效地改善了精度; 2)考慮到了湍流漩渦,提高了在這方面的精度; 3)RNG理論為湍流Prandtl數(shù)提供了一個解析公式,然而,標準k-ε模型使用的是用戶提供的常數(shù); 4)標準k-ε模型是一種高雷諾數(shù)的模型,RNG理論提供了一個考慮低雷諾數(shù)流動粘性的解析公式??蓪崿F(xiàn)的k-ε模型出現(xiàn)較晚,與標準k-ε模型相比有2個主要的不同點: 可實現(xiàn)的k-ε模型為湍流粘性增加了一個公式; 為耗散率增加了新的傳輸方程。標準k-ω模型是基于Wilcoxk-ω模型,為考慮低雷諾數(shù)、可壓縮性和剪切流傳播而修改的。標準的k-ε模型的一個變形就是SSTk-ω模型,這個模型由Menter發(fā)展,以便在廣泛的領(lǐng)域中可以獨立于k-ε模型,在近壁自由流中k-ω模型有廣泛的應(yīng)用范圍和精度。

    綜合考慮本項目的特點和不同湍流模型的特點,選擇標準k-ε模型公式進行求解計算[12-13]。

    湍流模型為

    (5)

    (6)

    式中:ε為湍流耗散率;Gk為因速度梯度引起的湍動能增加項;Gb為因浮力引起的湍動能k的增加項;ui為液相速度分量,m/s;μt為湍流黏性系數(shù);ρ為流體密度g/cm3;C1ε為經(jīng)驗常數(shù),一般分別取0.09、1.40、 1.92。

    3.2 邊界條件的設(shè)定

    由于我們模擬的是固體小顆粒在手套箱內(nèi)的運動,進氣口速度是已知的,因此進氣口處邊界條件選擇速度入口;出氣口處邊界條件選擇壓力出口,顆粒入口設(shè)置為wall-jet,其余的邊界設(shè)置為固體壁面邊界條件[14]。

    4 模擬結(jié)果分析

    4.1 氣固兩相速度分布

    為了更好地觀察箱體內(nèi)氣相和固相的速度分布,分別選取z為-150、 -580.5、 -101 1 mm;x為-500、 675.5、 1 840 mm截面為研究對象,分別得到當進氣口速度為4.777 m/s時,各個截面的速度離散圖和速度云圖。

    圖6為不同z值截面的速度離散圖。通過圖6 a)可以看到,進氣口速度為4.777m/s,氣體由進氣口進入流道之后,因為氣體的發(fā)散性,速度值逐漸減小,并且逐漸向周邊區(qū)域擴散。z=-150 mm右側(cè)為箱體壁面,箱體壁面有一定的粗糙度,導致速度值減小,速率較大;在z=-150 mm至z=1 840 mm區(qū)域內(nèi),氣相速度值變化較小,恒定為1 m/s左右;當處于流道出口時,速度再次增大,在流道末端速度值增加較快。由于增材制造設(shè)備產(chǎn)生固體顆粒占用一定的空間,會出現(xiàn)氣相在出氣口處速度大于進氣口處速度的現(xiàn)象。

    通過圖6 b)可以看出,在增材制造加工平臺,熱氣流上升帶動周圍固體小顆粒上升,隨后在高速氣體的帶動下,使固體顆粒速度明顯增大,從最初的0.6 m/s增大至2.5 m/s左右。之后在氣流組織的帶動下,固體顆粒伴隨氣流組織在箱體內(nèi)流動,整體速度與氣相速度正相關(guān),在箱體出口附近,氣相速度增大,則固體相顆粒速度逐漸增大,最終隨著氣相排出箱外。

    a)氣相b)固相圖6 不同z值截面的速度離散圖Fig.6 Discrete diagram of section velocity with different z values

    圖7為不同z值截面的速度云圖,圖中z為-150和-1 011 mm截面氣相在進氣口速度值最大,進入流道后,流域空間增加,氣體逐漸向四周擴散。在左側(cè)箱體壁面的作用下,氣體速度方向發(fā)生改變,整體分布表現(xiàn)為箱體壁面處速度值較大,中間區(qū)域速度較小。在出氣口位置附近,速度分布較為均勻,但是在近壁面附近的壁面粗糙,導致氣相速度較小,出氣口速度大于入氣口速度,這是因為在流場運動的過程中,不斷有固體顆粒的加入,導致流域體積空間發(fā)生變化。固相在進氣口附近速度值大于固相入口速度,其原因在于固相在流域內(nèi)游離,在進氣口附近高速氣體的帶動下,固相速度有所增加。在出氣口附近表現(xiàn)出固體顆粒粒徑較小,固相氣相伴隨性良好的特點。固相速度隨著氣相速度逐漸增加到最大,最后固體顆粒隨著高速氣體排出箱體外。在z=-580.5 mm截面,可以看到氣相和固相速度變化不大,基本趨于穩(wěn)定,在左側(cè)存在內(nèi)部設(shè)備,表現(xiàn)為氣相和固相在設(shè)備表面處速度值增大。固體顆粒釋放速度相對于進氣口氣相速度較小,所以在整個速度云圖中并不能明顯觀察到固相速度入口附近區(qū)域的速度變化。

    z=-150 mmz=-580.5 mmz=-1 011 mma)氣相z=-150 mmz=-580.5 mmz=-1 011 mmb)固相圖7 不同z值截面的速度云圖Fig.7 Velocity cloud diagram with different z-value cross section

    圖8為不同x值截面的速度離散圖,其中紅色區(qū)域為進氣口截面,黑色區(qū)域為出氣口截面,綠色區(qū)域為箱體中間對稱面。由圖中紅色區(qū)域可以看出,箱體在上、下端進氣口速度值呈喇叭形狀散開,速度減小,由于流道中設(shè)備占用空間,導致上端進氣口氣相速度衰減較為緩慢,下端進氣口速度快速衰減至一定值;黑色區(qū)域可以看出,箱體上、下端出氣口,速度分布較為均勻,表現(xiàn)出一定的規(guī)律性,由圖8 a)明顯可以看出,出氣口速度值大于進氣口速度值,進氣口速度為4.777 m/s,出氣口速度為5.7 m/s左右;從綠色區(qū)域可以看出,流道內(nèi)氣相在中心區(qū)域速度值較大,靠近箱體壁面處速度值較小,從壁面到中心區(qū)域速度逐漸增加,增加幅度較為平緩,中心區(qū)域固相速度值最大可以達到1 m/s。

    a)氣相b)固相圖8 不同x值截面的速度離散圖Fig.8 Discrete diagram of cross section velocity with different x values

    由圖8 b)中紅色區(qū)域可以看出,固相在上、下端進氣口速度值大于增材制造加工表面固相速度,箱體中部流域由于高速氣體的作用,部分固體顆粒表現(xiàn)出較高的速度;由黑色區(qū)域可看出,箱體上、下端出氣口,固相速度分布較為均勻,表現(xiàn)出一定的規(guī)律性,與氣相速度分布基本一致,呈現(xiàn)出固相隨氣相流動的現(xiàn)象;從綠色區(qū)域可以看出,流道內(nèi)中心區(qū)域速度值較大,靠近箱體壁面處速度值較小,從壁面到中心區(qū)域逐漸增加,增加幅度較為平緩,中心區(qū)域固相速度值最大可以達到1 m/s。

    圖9為氣相不同x值截面的速度云圖。由圖中可以看出,x=-500 mm截面處,由于內(nèi)部流道的布置,氣相速度主要沿箱體近壁面處分布,且速度衰減較快,流道中間流域速度較??;在進、出氣口附近區(qū)域,固相小顆粒伴隨氣相運動,固相顆粒速度小幅度增加。x=-675.5 mm截面處,氣相速度較小且變化不大,氣、固兩相速度在中部區(qū)域略大于壁面處。x=1 840 mm截面處,氣相在靠近出氣口位置速度均勻增加,到達出氣口時,速度值達到最大。在此區(qū)域固相速度隨氣相速度增加而增加,固相顆粒隨著氣體從出氣口排出箱體外部。

    x=-500 mmx=-675.5 mmx=1 840 mmx=-500 mmx=-675.5 mmx=1 840 mm圖9 氣相不同x值截面的速度云圖Fig.9 Cross section velocity nephogram of gas phase with different x values solid

    4.2 入口氣體速度對固相速度的影響

    為了探究進氣口氣體速度與箱體內(nèi)氣、 固兩相流動特性的關(guān)系,令進、 出氣口尺寸不變,多相流參數(shù)不變,通過改變?nèi)霘饪跉怏w速度,分別對氣、 固兩相流動進行數(shù)值模擬,最后對比分析在5種不同入口氣體速度下,箱體內(nèi)過點a(-500, 475.5, -580.5)、b(1 840, 475.5, -580.5)直線上的固相速度。工況1—5的氣體入口速度分別設(shè)置為1.592、 2.213、 2.654、 3.715、 4.777 m/s,氣體密度均設(shè)置為1.2 kg/m3,固體顆粒直徑設(shè)置為90 nm,固體顆粒密度設(shè)置為2 713 kg/m3。

    圖10為不同進氣口速度對直線ab上固相速度的影響,由5種工況下的固相速度對比圖可看出,隨著進氣口速度逐漸增加,固相速度在箱體內(nèi)速度變化規(guī)律類似,具體表現(xiàn)為:x=-500至x=-400 mm區(qū)域,固相第1次加速;x=-400至x=-100 mm區(qū)域,固相顆粒上升力與空氣阻力逐漸達到平衡狀態(tài),使固相速度區(qū)域平緩,并且有略微的減速現(xiàn)象;x=-100至x=300 mm區(qū)域,固相顆粒物第2次加速,速度分別增加至對應(yīng)最大值;x=300至x=1 500 mm區(qū)域,固相顆粒進入流道中間區(qū)域,速度值減小??諝庾枇εc固相顆粒速度成正比,所以進氣口氣相速度值越大,該速度衰減速度越快;x=1 500至x=1 800 mm區(qū)域的內(nèi)部設(shè)備占用流道空間,導致固相顆??焖偎p至一個較小的值,并且這個區(qū)域存在進氣口氣相速度值越大,速度衰減越快的現(xiàn)象。

    圖10 不同進氣口速度對直線ab上固相速度的影響Fig.10 Effect of different inlet velocities on solid phase velocity at straight line ab

    4.3 固相顆粒粒徑對固相速度的影響

    在箱體實際應(yīng)用過程中,顆粒物會凝并,左側(cè)3D打印區(qū)域產(chǎn)生的粉體粒徑大小會有差異,故此,為了更加深入地研究箱體內(nèi)部粉塵的運動規(guī)律,達到理想的除塵效果,在保證流場分析參數(shù)不變的情況下,改變顆粒物粒徑,對箱體內(nèi)流域中的固相顆粒速度進行分析,工況6—10的固體顆粒直徑分別設(shè)置為0.09、 0.9、 9、 5、 10 μm,氣體入口速度均設(shè)置為1.592 m/s,氣體密度設(shè)置為1.2 kg/m3。

    對5種工況進行模擬,得到5種工況下固體顆粒在ab線上的速度對比圖,見圖11。由圖可以看出,當進氣口速度v=1.592 m/s恒定時,固體顆粒在線ab上的速度總體上隨著固體顆粒粒徑的增大而減小。當顆粒直徑d小于5 μm時,表現(xiàn)為在x=-500至x=-400 mm區(qū)域,固相速度增加; 在x=-400至x=-100 mm區(qū)域,由于固相顆粒上升力與空氣阻力共同作用,固相速度出現(xiàn)小范圍減速; 在x=-100至x=300 mm區(qū)域,固相顆粒物速度再次增加; 在x=300至x=1 500 mm區(qū)域,速度值緩慢減??;x=1 500 mm至x=1 840 mm區(qū)域,由于左側(cè)箱壁的作用,固相速度快速減小。當顆粒直徑d≥10 μm時,在x=-500至x=-300 mm區(qū)域,固相顆粒物速度值急劇減?。?在x=-300 mm至x=600 mm區(qū)域,速度隨著高速氣流作用增加,固相顆粒物速度快速增加; 在x=600 mm至x=1 500 mm區(qū)域,速度值緩慢減小,由于慣性較大,速度減小較為緩慢;x=1 500 mm至x=1 840 mm區(qū)域,與顆粒直徑d≤5 μm的固體顆粒相同。綜上所述,當3D打印時產(chǎn)生的燒結(jié)顆粒直徑為5 μm以下時,設(shè)定的初始空氣流量依然可以滿足箱體內(nèi)的除塵要求。

    圖11 5種工況下固體顆粒在ab線上的速度對比圖Fig.11 Velocity comparison of solid particles on abline under five working conditions

    4.4 實驗驗證

    依據(jù)已構(gòu)建的物理模型,進行箱體強度靜力學分析,干涉性校驗之后,搭建實驗平臺,對氣固兩相流動仿真分析出的現(xiàn)象進行工程實際驗證,1∶1實物實驗平臺如圖12所示。

    a)正視圖b)側(cè)視圖圖12 1∶1實物試驗平臺Fig.12 1∶1 Physical test platform

    H型真空手套箱主要為產(chǎn)品提供低含水量、低含氧量、無塵的的加工環(huán)境,通常水、氧的體積分數(shù)均小于10-6,進而保證箱體內(nèi)產(chǎn)品質(zhì)量以及產(chǎn)品材料的原有屬性。當固相顆粒粒徑較小時,固相顆粒物對氣相有較好的伴隨性,所以本實驗通過箱體內(nèi)氧含量的變化,反映理論流體仿真分析結(jié)果的可靠性。通過實際測試,手套箱水、氧含量變化趨勢圖如圖13所示。

    a)第1次檢測圖b)第2次檢測圖圖13 氧含量和水含量變化趨勢圖Fig.13 Trend chart of oxygen content and water content

    圖13中,第1次檢測氧含量在1 h內(nèi)一直處于波動狀態(tài),其中前30 min氧含量變化的幅度較大;第2次檢測氧含量波動頻率較小,在長時間內(nèi)氧含量不能達到穩(wěn)定狀態(tài)。通過2次箱體內(nèi)部氧含量檢測試驗,驗證箱體數(shù)值模擬結(jié)果,即當前手套箱系統(tǒng)內(nèi)部流場流動不均勻、局部空氣渦流,使水、氧含量不能在較短時間內(nèi)達到生產(chǎn)制造的氛圍要求。

    5 結(jié)論

    1)對H型特殊結(jié)構(gòu)手套箱內(nèi)部氣固兩相速度進行數(shù)值模擬,分析在生產(chǎn)過程中不能快速達到加工標準氛圍原因,通過關(guān)鍵位置速度云圖、速度離散圖發(fā)現(xiàn),箱體內(nèi)速度分布不均勻,局部產(chǎn)生較嚴重渦流區(qū),進而表現(xiàn)為局部煙霧、或者局部顆粒物沉積。

    2)通過對箱體內(nèi)工況進行分析,發(fā)現(xiàn)當固相顆粒物粒徑(d=90 nm)不變時,隨著進氣口速度的增加,固相速度隨之增大,固相對氣相表現(xiàn)出良好的伴隨性。

    3)考慮到顆粒物的凝并現(xiàn)象,研究了當進氣口速度恒定,依次增加顆粒物粒徑,分析固相速度得到,當顆粒粒徑小于5 μm,固相對氣相表現(xiàn)出較好的伴隨性,當顆粒粒徑大于等于10 μm時,固相顆粒物對氣相的伴隨性較差。

    4)通過分析不同參數(shù)變化對氣相、固相速度的影響,闡釋H型手套箱內(nèi)氣固兩相流動特性,雖然不能直接適用于所有手套箱設(shè)備,但對于類似問題研究具有一定參考價值。

    猜你喜歡
    出氣口進氣口箱體
    閥門及氣相沉積設(shè)備
    低溫與特氣(2022年5期)2022-03-23 17:49:51
    基于仿真的空濾出氣口圓角尺寸對流阻的影響
    一種SBS橡膠尾氣收集處理裝置
    高牌號灰鐵前端箱體質(zhì)量提升
    飛機輔助進氣口小流量工況的性能初探
    科技視界(2017年13期)2017-09-30 16:20:01
    北極熊的運氣
    他山之石
    航空模型(2016年3期)2016-05-14 09:00:58
    超大型冷剪箱體加工難點分析
    基于ANSYS Workbench 的ATB260 減速器箱體模態(tài)分析
    一款箱體可整體收縮折疊式簾布半掛車
    專用汽車(2016年9期)2016-03-01 04:17:30
    一区二区三区高清视频在线| 好男人视频免费观看在线| 精品午夜福利在线看| 亚洲欧美精品专区久久| 91午夜精品亚洲一区二区三区| 三级国产精品欧美在线观看| 黄色一级大片看看| 特级一级黄色大片| 尾随美女入室| 一级毛片aaaaaa免费看小| 久久久久久久久久久免费av| 亚洲第一电影网av| 特大巨黑吊av在线直播| 秋霞在线观看毛片| 最后的刺客免费高清国语| 欧美zozozo另类| 狠狠狠狠99中文字幕| 观看美女的网站| 激情 狠狠 欧美| 久久久精品欧美日韩精品| 国产一区二区在线av高清观看| 国产成年人精品一区二区| 成人三级黄色视频| 久久久午夜欧美精品| 国产精品久久久久久亚洲av鲁大| 午夜福利在线观看免费完整高清在 | 精品欧美国产一区二区三| 亚洲av熟女| 99国产精品一区二区蜜桃av| 两个人视频免费观看高清| 日本爱情动作片www.在线观看| 亚洲av一区综合| 青春草亚洲视频在线观看| 国产视频首页在线观看| 免费大片18禁| 偷拍熟女少妇极品色| 免费大片18禁| 一区二区三区四区激情视频 | 国产欧美日韩精品一区二区| 国产精品人妻久久久久久| 国产伦理片在线播放av一区 | 日韩精品青青久久久久久| 婷婷亚洲欧美| 黄色配什么色好看| 国产精品.久久久| 变态另类成人亚洲欧美熟女| 亚洲最大成人中文| 亚洲国产精品合色在线| 在线a可以看的网站| 全区人妻精品视频| 精品熟女少妇av免费看| 中文亚洲av片在线观看爽| 欧美日韩一区二区视频在线观看视频在线 | 99riav亚洲国产免费| 国产片特级美女逼逼视频| 三级毛片av免费| 亚洲一区高清亚洲精品| 人妻久久中文字幕网| 国产综合懂色| 精品人妻一区二区三区麻豆| 99热这里只有是精品50| 99精品在免费线老司机午夜| 啦啦啦观看免费观看视频高清| 联通29元200g的流量卡| 亚洲性久久影院| 伦理电影大哥的女人| 欧美+日韩+精品| 国产精品一区二区三区四区免费观看| 亚洲经典国产精华液单| 你懂的网址亚洲精品在线观看 | 精品不卡国产一区二区三区| 成人亚洲欧美一区二区av| 插逼视频在线观看| 日韩欧美 国产精品| 两个人的视频大全免费| 午夜福利在线在线| 亚洲自拍偷在线| 国产不卡一卡二| 99久久精品热视频| 丰满人妻一区二区三区视频av| 精品久久久久久久久av| 国产精品免费一区二区三区在线| 国产黄a三级三级三级人| 综合色av麻豆| 国产精品女同一区二区软件| 超碰av人人做人人爽久久| 日韩成人av中文字幕在线观看| 日韩成人av中文字幕在线观看| 美女cb高潮喷水在线观看| 国产黄色视频一区二区在线观看 | 99热这里只有是精品在线观看| 蜜桃亚洲精品一区二区三区| 国产91av在线免费观看| 午夜久久久久精精品| 欧美高清性xxxxhd video| 精品人妻一区二区三区麻豆| 久久久a久久爽久久v久久| 亚洲三级黄色毛片| 免费观看精品视频网站| av在线亚洲专区| 网址你懂的国产日韩在线| 欧美一区二区精品小视频在线| 亚洲乱码一区二区免费版| 精品国产三级普通话版| 赤兔流量卡办理| 日韩大尺度精品在线看网址| 成人午夜精彩视频在线观看| av在线亚洲专区| 一级黄色大片毛片| 国内少妇人妻偷人精品xxx网站| 国产极品天堂在线| 久久99蜜桃精品久久| 亚洲综合色惰| 国产成人福利小说| 偷拍熟女少妇极品色| 亚洲精品av麻豆狂野| 亚洲成人手机| 亚洲综合精品二区| 亚洲精品国产av成人精品| 亚洲欧美一区二区三区国产| 国产亚洲一区二区精品| 久久久久久久大尺度免费视频| 老熟女久久久| 日本欧美视频一区| 91精品三级在线观看| 九草在线视频观看| 成人毛片60女人毛片免费| 日韩一本色道免费dvd| 汤姆久久久久久久影院中文字幕| 欧美bdsm另类| 国产一区亚洲一区在线观看| 青春草亚洲视频在线观看| 精品卡一卡二卡四卡免费| 亚洲欧美日韩另类电影网站| 亚洲av不卡在线观看| 久久国产精品男人的天堂亚洲 | 国产日韩欧美在线精品| 精品一品国产午夜福利视频| 亚洲国产日韩一区二区| 免费观看在线日韩| 婷婷色麻豆天堂久久| 美女国产视频在线观看| 日韩强制内射视频| 啦啦啦视频在线资源免费观看| 亚洲欧美清纯卡通| 香蕉精品网在线| 国产淫语在线视频| 亚洲精品美女久久av网站| 黑丝袜美女国产一区| 亚洲欧美一区二区三区黑人 | 嘟嘟电影网在线观看| 亚洲欧美色中文字幕在线| 亚洲,欧美,日韩| 国产白丝娇喘喷水9色精品| 精品国产露脸久久av麻豆| 午夜日本视频在线| 成人国产麻豆网| 狂野欧美激情性xxxx在线观看| 午夜91福利影院| 久久精品久久久久久久性| 国产成人精品一,二区| 午夜视频国产福利| .国产精品久久| 日本与韩国留学比较| 成人国语在线视频| 一级二级三级毛片免费看| 国产在视频线精品| 欧美丝袜亚洲另类| 亚洲第一av免费看| 亚洲无线观看免费| 国产成人精品久久久久久| av播播在线观看一区| 亚洲av成人精品一二三区| 亚洲精品色激情综合| 欧美国产精品一级二级三级| 午夜福利影视在线免费观看| 欧美 日韩 精品 国产| 涩涩av久久男人的天堂| 亚洲人成网站在线观看播放| av有码第一页| 国产黄频视频在线观看| 777米奇影视久久| 少妇猛男粗大的猛烈进出视频| 美女国产高潮福利片在线看| 久久久久久久久久久丰满| 国产免费一区二区三区四区乱码| 免费av不卡在线播放| 亚洲性久久影院| 国产午夜精品一二区理论片| 国产午夜精品一二区理论片| 尾随美女入室| 精品一品国产午夜福利视频| 热99久久久久精品小说推荐| 80岁老熟妇乱子伦牲交| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 99久久人妻综合| 99久久综合免费| 久久久久久久精品精品| 一区二区三区四区激情视频| 校园人妻丝袜中文字幕| 男女边摸边吃奶| 两个人免费观看高清视频| 亚洲图色成人| 精品人妻一区二区三区麻豆| 欧美一级a爱片免费观看看| 久久久久久久亚洲中文字幕| 极品少妇高潮喷水抽搐| 午夜精品国产一区二区电影| 狠狠精品人妻久久久久久综合| 亚洲美女搞黄在线观看| 极品少妇高潮喷水抽搐| 男女无遮挡免费网站观看| 欧美老熟妇乱子伦牲交| 欧美老熟妇乱子伦牲交| 赤兔流量卡办理| 色婷婷久久久亚洲欧美| 中文字幕精品免费在线观看视频 | 夜夜看夜夜爽夜夜摸| 国产高清不卡午夜福利| 日本黄大片高清| 亚洲精品日韩在线中文字幕| 欧美日韩在线观看h| 草草在线视频免费看| 国产黄频视频在线观看| 久久精品国产自在天天线| 久久精品国产自在天天线| 日本与韩国留学比较| 欧美人与性动交α欧美精品济南到 | 人妻夜夜爽99麻豆av| 久久免费观看电影| 制服诱惑二区| 日韩精品有码人妻一区| 国产男人的电影天堂91| 久久影院123| 久久久国产一区二区| 国产精品一区二区在线不卡| 国模一区二区三区四区视频| 美女cb高潮喷水在线观看| 精品亚洲成国产av| 日韩熟女老妇一区二区性免费视频| 免费人妻精品一区二区三区视频| 美女主播在线视频| 黑人欧美特级aaaaaa片| 午夜福利视频精品| 国语对白做爰xxxⅹ性视频网站| 18禁在线播放成人免费| 久久久久网色| 国产亚洲一区二区精品| 乱人伦中国视频| 成年女人在线观看亚洲视频| 精品一品国产午夜福利视频| 三级国产精品欧美在线观看| 午夜日本视频在线| 99热国产这里只有精品6| 七月丁香在线播放| 欧美少妇被猛烈插入视频| 久久精品久久精品一区二区三区| 国产精品国产三级国产专区5o| 香蕉精品网在线| 午夜影院在线不卡| 午夜福利,免费看| 大香蕉久久网| 精品熟女少妇av免费看| 不卡视频在线观看欧美| av卡一久久| 日日撸夜夜添| 亚州av有码| 免费不卡的大黄色大毛片视频在线观看| 啦啦啦在线观看免费高清www| 欧美另类一区| 日韩大片免费观看网站| 精品久久久久久久久av| 成人影院久久| 嫩草影院入口| 国产片特级美女逼逼视频| 丝袜美足系列| 一级毛片我不卡| 国产黄片视频在线免费观看| 丰满迷人的少妇在线观看| 9色porny在线观看| 亚洲国产av影院在线观看| 91aial.com中文字幕在线观看| 少妇丰满av| 国产亚洲欧美精品永久| 国产男女内射视频| 国产高清国产精品国产三级| 美女cb高潮喷水在线观看| 晚上一个人看的免费电影| 女的被弄到高潮叫床怎么办| 国产在视频线精品| 韩国高清视频一区二区三区| 欧美日韩亚洲高清精品| 亚洲av在线观看美女高潮| 九九爱精品视频在线观看| 国产淫语在线视频| 看免费成人av毛片| 91久久精品电影网| 亚洲熟女精品中文字幕| 男女边吃奶边做爰视频| 国产在视频线精品| 又黄又爽又刺激的免费视频.| 亚洲精品美女久久av网站| 爱豆传媒免费全集在线观看| 国产成人aa在线观看| 欧美人与善性xxx| 纵有疾风起免费观看全集完整版| 久久99一区二区三区| 成年美女黄网站色视频大全免费 | 特大巨黑吊av在线直播| 欧美日韩综合久久久久久| 啦啦啦在线观看免费高清www| 国产爽快片一区二区三区| 欧美日韩综合久久久久久| 97超视频在线观看视频| 又黄又爽又刺激的免费视频.| 亚洲国产av新网站| 精品亚洲乱码少妇综合久久| 亚洲伊人久久精品综合| 涩涩av久久男人的天堂| 成人国产av品久久久| 天天操日日干夜夜撸| 一本大道久久a久久精品| 国产男女超爽视频在线观看| 午夜老司机福利剧场| 国产探花极品一区二区| 母亲3免费完整高清在线观看 | 国产精品人妻久久久影院| 丝瓜视频免费看黄片| 国产一区二区在线观看日韩| 成人二区视频| 观看av在线不卡| 99热这里只有是精品在线观看| 男女边摸边吃奶| av有码第一页| 热99久久久久精品小说推荐| 亚洲欧洲国产日韩| 简卡轻食公司| 亚洲国产精品一区二区三区在线| 色视频在线一区二区三区| 中国美白少妇内射xxxbb| 一区二区av电影网| 香蕉精品网在线| 午夜福利网站1000一区二区三区| 国产综合精华液| 亚洲无线观看免费| 18禁裸乳无遮挡动漫免费视频| 黄色视频在线播放观看不卡| 亚洲精品视频女| 久久 成人 亚洲| 国产日韩欧美视频二区| 日产精品乱码卡一卡2卡三| 在线观看免费高清a一片| 赤兔流量卡办理| 欧美bdsm另类| 一级黄片播放器| 久久精品人人爽人人爽视色| 男男h啪啪无遮挡| 久久国产精品男人的天堂亚洲 | 永久网站在线| 免费大片18禁| 一级毛片aaaaaa免费看小| 国产黄频视频在线观看| 久久久精品区二区三区| 国产熟女欧美一区二区| 亚洲精品久久午夜乱码| 在线观看一区二区三区激情| 亚州av有码| 91精品国产九色| 免费观看在线日韩| 校园人妻丝袜中文字幕| 国产成人freesex在线| 爱豆传媒免费全集在线观看| 国产欧美日韩一区二区三区在线 | 一级黄片播放器| 亚洲中文av在线| 日韩中文字幕视频在线看片| av在线老鸭窝| 黄片无遮挡物在线观看| 久久精品夜色国产| 国产精品久久久久久久电影| 亚洲av电影在线观看一区二区三区| 亚洲av在线观看美女高潮| 亚洲精品乱码久久久久久按摩| 午夜福利视频精品| 久久久精品免费免费高清| 黑人猛操日本美女一级片| 99热这里只有精品一区| 18在线观看网站| 午夜福利视频精品| 精品久久久噜噜| 亚洲内射少妇av| 成年av动漫网址| 三级国产精品欧美在线观看| 免费观看无遮挡的男女| 成人亚洲欧美一区二区av| 丰满少妇做爰视频| 日本免费在线观看一区| 免费观看的影片在线观看| 26uuu在线亚洲综合色| 老司机影院成人| 欧美精品一区二区大全| 视频中文字幕在线观看| 欧美变态另类bdsm刘玥| 日本色播在线视频| 18禁动态无遮挡网站| 男女边摸边吃奶| 国产视频内射| 美女大奶头黄色视频| 亚洲精品第二区| 插逼视频在线观看| 女性被躁到高潮视频| 黄色一级大片看看| 性色avwww在线观看| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 欧美xxxx性猛交bbbb| 亚洲精品视频女| 免费观看的影片在线观看| 寂寞人妻少妇视频99o| 美女内射精品一级片tv| 精品亚洲乱码少妇综合久久| 午夜激情久久久久久久| 亚洲欧美一区二区三区国产| 成人毛片a级毛片在线播放| 久久午夜福利片| 成年美女黄网站色视频大全免费 | 十八禁高潮呻吟视频| 99久久中文字幕三级久久日本| 亚洲成人一二三区av| 午夜福利网站1000一区二区三区| 亚洲av电影在线观看一区二区三区| 一个人看视频在线观看www免费| 夫妻午夜视频| tube8黄色片| 自拍欧美九色日韩亚洲蝌蚪91| 99久久人妻综合| 国产 一区精品| 久久97久久精品| 成人亚洲欧美一区二区av| 久久综合国产亚洲精品| 一区二区日韩欧美中文字幕 | 最新的欧美精品一区二区| 啦啦啦啦在线视频资源| 母亲3免费完整高清在线观看 | 亚洲第一区二区三区不卡| 亚洲国产欧美在线一区| 成年女人在线观看亚洲视频| 一区二区av电影网| 欧美97在线视频| 亚洲欧洲国产日韩| 亚洲国产毛片av蜜桃av| 亚洲怡红院男人天堂| 老女人水多毛片| 精品一区二区免费观看| 国产亚洲一区二区精品| 精品一区二区免费观看| 香蕉精品网在线| 99热网站在线观看| 久久久久久久大尺度免费视频| 99热网站在线观看| 亚洲国产成人一精品久久久| av免费在线看不卡| 亚洲av男天堂| 少妇 在线观看| 亚洲成色77777| 桃花免费在线播放| 欧美xxxx性猛交bbbb| 日本av免费视频播放| 国产精品一区www在线观看| 免费观看在线日韩| 99热这里只有是精品在线观看| 久久99一区二区三区| 插阴视频在线观看视频| 成人午夜精彩视频在线观看| 久久99热这里只频精品6学生| 久久毛片免费看一区二区三区| 黑人高潮一二区| 80岁老熟妇乱子伦牲交| 国产男女内射视频| 涩涩av久久男人的天堂| 亚洲熟女精品中文字幕| 国产日韩欧美视频二区| 亚洲av中文av极速乱| 9色porny在线观看| 国产高清不卡午夜福利| 91精品国产国语对白视频| 亚洲精华国产精华液的使用体验| av卡一久久| 一个人看视频在线观看www免费| 日本av免费视频播放| 国产精品国产三级国产专区5o| 亚洲精品日韩在线中文字幕| 老司机亚洲免费影院| 日本-黄色视频高清免费观看| 国产精品国产三级国产av玫瑰| 国产黄色视频一区二区在线观看| 欧美bdsm另类| 国产日韩欧美在线精品| 啦啦啦在线观看免费高清www| 人人妻人人爽人人添夜夜欢视频| 国产精品女同一区二区软件| 九九在线视频观看精品| videossex国产| 久久免费观看电影| 欧美日韩av久久| 在线看a的网站| 日本av手机在线免费观看| 国产精品国产av在线观看| 777米奇影视久久| 啦啦啦啦在线视频资源| 大香蕉久久成人网| 蜜桃国产av成人99| 久久久久久伊人网av| 伊人久久国产一区二区| 国产视频首页在线观看| kizo精华| 一区二区三区乱码不卡18| 久久久国产精品麻豆| 丝袜脚勾引网站| 插逼视频在线观看| 国产精品国产三级国产专区5o| 伊人亚洲综合成人网| 午夜精品国产一区二区电影| 成人黄色视频免费在线看| 热re99久久精品国产66热6| 国产免费现黄频在线看| 一级毛片我不卡| 91精品一卡2卡3卡4卡| 国产成人a∨麻豆精品| videosex国产| 国产高清有码在线观看视频| 肉色欧美久久久久久久蜜桃| 日韩欧美精品免费久久| 日韩熟女老妇一区二区性免费视频| 美女xxoo啪啪120秒动态图| 亚洲无线观看免费| 色视频在线一区二区三区| 国产高清不卡午夜福利| 久久国内精品自在自线图片| 十八禁网站网址无遮挡| 在线 av 中文字幕| 久久99一区二区三区| 夜夜骑夜夜射夜夜干| 黄片无遮挡物在线观看| 啦啦啦啦在线视频资源| 2022亚洲国产成人精品| 韩国av在线不卡| 亚洲精品日韩在线中文字幕| 色婷婷久久久亚洲欧美| 国产一区二区在线观看日韩| 久久久久久久大尺度免费视频| 视频在线观看一区二区三区| 亚洲精品乱码久久久v下载方式| 精品一区二区免费观看| 99热这里只有精品一区| 又大又黄又爽视频免费| 亚洲国产精品国产精品| 国产日韩一区二区三区精品不卡 | 22中文网久久字幕| 国产老妇伦熟女老妇高清| 欧美日韩综合久久久久久| a级毛色黄片| 免费人妻精品一区二区三区视频| 国产视频内射| 在线观看三级黄色| 精品一区二区三卡| 成人国产av品久久久| 亚洲人成77777在线视频| 91久久精品国产一区二区成人| 亚洲人成77777在线视频| 建设人人有责人人尽责人人享有的| 大香蕉久久成人网| 国产精品一区www在线观看| 亚洲av日韩在线播放| 欧美成人精品欧美一级黄| 欧美日韩视频精品一区| 国产无遮挡羞羞视频在线观看| 又粗又硬又长又爽又黄的视频| 大又大粗又爽又黄少妇毛片口| 亚洲国产欧美在线一区| 精品人妻一区二区三区麻豆| 日韩制服骚丝袜av| av线在线观看网站| 国产 一区精品| av专区在线播放| 色哟哟·www| 久久婷婷青草| 一级a做视频免费观看| 日韩在线高清观看一区二区三区| 人妻人人澡人人爽人人| 99国产综合亚洲精品| 日韩欧美精品免费久久| 日产精品乱码卡一卡2卡三| 99热全是精品| 日产精品乱码卡一卡2卡三| 韩国高清视频一区二区三区| 人人妻人人爽人人添夜夜欢视频| 校园人妻丝袜中文字幕| 亚洲成色77777| 另类精品久久| 26uuu在线亚洲综合色| 亚洲丝袜综合中文字幕| 欧美国产精品一级二级三级| 青春草视频在线免费观看| 国产一区亚洲一区在线观看| 女性生殖器流出的白浆| 能在线免费看毛片的网站| 最近中文字幕2019免费版| 国内精品宾馆在线| 一级,二级,三级黄色视频| 一个人免费看片子| 国产精品欧美亚洲77777| 久久精品熟女亚洲av麻豆精品| 成人无遮挡网站| 国产成人精品久久久久久| 热99国产精品久久久久久7|