亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        筒體直徑對旋風分離器性能的影響

        2017-08-31 12:22:28孫國剛周發(fā)戚孫占朋
        石油學報(石油加工) 2017年4期
        關鍵詞:旋風筒體分離器

        袁 怡, 孫國剛, 周發(fā)戚, 孫占朋

        (中國石油大學 化學工程學院 過程流體過濾與分離北京市重點實驗室, 北京 102249)

        筒體直徑對旋風分離器性能的影響

        袁 怡, 孫國剛, 周發(fā)戚, 孫占朋

        (中國石油大學 化學工程學院 過程流體過濾與分離北京市重點實驗室, 北京 102249)

        筒體直徑是影響旋風分離器效率和壓降的重要結構參數(shù)。以硅微粉為原料,在保證幾何結構相似的基礎上,對筒體直徑分別為200、300、400 mm的Stairmand型旋風分離器進行冷態(tài)對比實驗,考察了不同情形下筒徑對旋風分離器分離性能的影響。結果表明,旋風分離器按幾何相似放大,筒徑增大,在同一入口氣速下,分離效率下降而壓降升高;在同一處理氣量下,分離效率和壓降都下降。如將入口尺寸與處理氣量固定,其它尺寸按幾何相似放大,筒徑增大則旋風分離器的分離效率增加、壓降降低。進而又用幾個較新的分離模型對實驗進行了計算比較,發(fā)現(xiàn)模型計算結果與實驗結果的趨勢大體一致,基本上能預測筒徑變化對旋風分離器分離效率與壓降的影響趨勢,但都不能定量計算筒徑對其分離性能的影響。

        旋風分離器;筒體直徑;效率;壓降;模型計算

        旋風分離器是氣-固分離過程中一種常用的設備,因其結構簡單,無運轉部件,適用于高溫高壓等條件,目前在石油、化工、冶金、能源、環(huán)境等諸多領域得到廣泛應用。要將在小模型實驗研究中開發(fā)出的高效旋風分離器應用到工業(yè)中去,必須經過放大設計和工業(yè)設計,但小型實驗設計經過放大后其性能如何變化,迄今并無成熟的計算方法[1]。在保證幾何相似的基礎上,對旋風分離器尺寸放大可通過增大筒體直徑來實現(xiàn)。一般認為:隨著筒體直徑的增大,旋風分離器的分離效率降低,壓降升高[2]。然而,筒體直徑對旋風分離器性能的影響取決于多種因素。

        金有海等[3]對蝸殼式旋風分離器進行相似放大的試驗研究,認為筒體直徑越小,離心力場越強,從而粒子被分離的效率越高。Jin等[4]采用CFD數(shù)值模擬的方法,發(fā)現(xiàn)當質量流率和入口截面一定,增大筒體直徑,旋風分離器的總分離效率提高,壓降幾乎不變,流場更穩(wěn)定。袁惠新等[5]、方麒先[6]認為在入口速度一定的條件下,旋風分離器的壓降隨著筒徑的增大而增大。劉書賢等[7]對催化裂化裝置沉降器粗旋在尺寸放大后的壓降進行數(shù)值模擬研究,認為尺寸放大后粗旋的壓降減小。宋健斐等[8]對旋風分離器直徑尺寸變化時的氣相流場進行數(shù)值模擬研究,認為隨著筒體直徑的增大,流場流態(tài)保持相似性。迄今為止,前人雖然做了不少關于筒徑對旋風分離器性能影響的研究工作,但大多集中于采用數(shù)值模擬的研究手段,且其得出的結論并不全面,很難對工程實踐形成有效指導。因此,尋找筒體直徑對旋風分離器性能的影響規(guī)律,就成為迫切需要解決的課題。

        在本研究中,以Stairmand型旋風分離器[9]為基準模型,選擇結構參數(shù)和入口氣速作為變量,分別以理論模型計算及實驗測量的手段考察旋風分離器幾何放大后性能的變化規(guī)律,以期獲得更為全面的認識,從而為旋風分離器的優(yōu)化設計及工業(yè)應用提供有益參考。

        1 試驗裝置與方法

        采用Stairmand型旋風分離器為基準模型,筒徑(D)分別為200 mm、300 mm、400 mm,下接料腿(L)為1500 mm,用有機玻璃制造,試驗所用旋風分離器尺寸結構如圖1所示。圖1(a)中入口截面系數(shù)KA是1個與筒徑及入口高寬有關的無量綱數(shù)[10],KA=7.85,各部分尺寸隨筒徑等比例變化;圖1(b)中入口高×寬為144 mm×58 mm,排氣管插入深度為144 mm,其他各部分尺寸隨筒徑等比例變化。

        實驗介質為常溫空氣,粉料為硅微粉,顆粒密度約為2600 kg/m3,采用BT-9300S 型激光粒度分布儀進行粒度分析,中位粒徑為10.94 μm,粒徑分布dp如圖2所示。試驗在吸風負壓狀態(tài)下操作,試驗流程如圖3所示,入口氣速由標準畢托管測量;風量通過閥門調節(jié);粉料從分離器入口處的加料斗人工加料,加料質量濃度為10 g/m3,通過加料時間來控制加料濃度;分離器壓降用U型管測量,分離效率采用加料、收料稱重法計算。分離器入口段有一定長度,使粉料和氣體在分離器入口段充分混合,混合后的攜塵氣體經過旋風分離器的分離后由排氣管排出到風機外排系統(tǒng)。

        圖1 旋風分離器結構尺寸Fig.1 Critical dimensions of cyclone separator(a) Stairmand cyclone separator with inlet cross-section coefficient KA=7.85;(b) Stairmand cyclone separator with inlet size 144 mm×58 mm

        圖2 實驗用硅微粉粒徑分布Fig.2 Experimental silica powder size distribution

        圖3 實驗裝置示意圖Fig.3 Schematic drawing of the experimental setup

        2 試驗結果及分析

        2.1 筒體直徑對旋風分離器效率的影響

        采用圖1(a)Stairmand型分離器,入口截面系數(shù)相同(KA=7.85),筒體直徑D分別為200、300、400 mm。(1)選取10、15、20、25 m/s 4個氣速,由試驗得到該組旋風分離器在不同氣速下的分離效率如圖4所示。從圖4可以看出,入口氣速相同,筒徑等幾何尺寸等比例放大,旋風分離器分離效率降低;還可以發(fā)現(xiàn),隨著入口速度的增加,旋風分離器效率會出現(xiàn)轉折點,存在1個最大效率速度,且筒徑的變化對該最大效率入口速度的大小幾乎沒有影響。這是因為入口氣速增大,離心力場增強,效率增加;當氣速增大到一定值后,湍流及顆粒碰撞彈跳等導致沉積在器壁處的一些顆粒重新被卷揚起來,加上灰斗夾帶增多等不利因素的綜合影響,導致旋風分離器效率有所降低。(2)選取處理氣量0.167 m3/s,由試驗得到此時該組旋風分離器在不同筒徑下的效率如圖5所示。從圖5可以看出,同一入口氣量時,幾何尺寸按同一比例放大,隨著筒徑增大,入口截面尺寸增大,入口氣速下降,旋風分離器內離心力場減弱,從而導致分離器的分離效率降低。

        圖4 分離效率(η)-入口氣速(Vin)曲線(KA=7.85)Fig.4 Efficiency(η)-gas velocity(Vin) curve (KA=7.85)

        圖5 不同筒徑(D)下的旋風分離器分離效率(η)Fig.5 The efficiency(η) of cyclone under varied cylinder diameters(D)

        2.2 筒體直徑對旋風分離器壓降的影響

        筒體直徑對旋風分離器壓降的影響如圖6所示。從圖6可以看出,當旋風分離器幾何尺寸按比例放大,同一入口氣速下,隨著筒體直徑的增加壓降增大;同一處理氣量下,隨著筒體直徑的增加壓降降低;當固定入口截面尺寸大小,旋風分離器其他尺寸按比例放大,則入口氣速及處理氣量均相同,旋風分離器的壓降隨著筒徑的增大而減小。

        圖6 不同筒徑(D)下旋風分離器的總壓降(Δp)Fig.6 The pressure drop(Δp) of the cyclone under varied cylinder diameters(D) (a) The same inlet velocity; (b) The same consumption flow

        3 模型預測結果與分析

        3.1 旋風分離器分離模型的預測結果分析

        作為旋風分離器分離效率的預測模型即分離模型,從20世紀50年代以后在國內外取得廣泛的研究。各研究者基于不同假設,從不同的角度闡述了各自理論,歸納起來主要有“平衡軌道”模型、“停留時間”模型以及介于兩者間的混合模型。表1列出了幾種近年來較具代表性的旋風分離器分離模型具體效率計算方法。

        表1 旋風分離器分離效率的計算模型Table 1 Calculation model of separation efficiency for cyclones

        筆者以實驗測量的方法對不同筒徑下的Stairmand型旋風分離器的分離效率進行了研究,分離模型也可以預測分離器的效率。在工業(yè)過程中,總分離效率通常是1個最常用的評價指標,將表1中模型預測的切割粒徑或粒級效率換算成總效率,總效率是通過這樣方式計算的:首先按照顆粒粒徑分布劃分為N個粒度級的質量組分,每個質量組分之和構成了全部顆粒質量;然后,用每個質量分數(shù)乘以該組分平均粒徑下的粒級效率,所以N個粒度級組分的總和就得到總效率。數(shù)學表達式為:

        (1)

        式中,ηi為第i個組分的粒級效率,wi為第i個組分的質量分數(shù)。

        圖7為旋風分離器總效率的實驗值與模型計算值對比。圖7(a)、(b)分別為入口氣速15 m/s及20 m/s 下,幾何尺寸相似放大的旋風分離器分離效率隨筒徑的變化關系;圖7(c)為同一處理氣量下,幾何尺寸相似放大的旋風分離器分離效率隨筒徑的變化關系;圖7(d)同一處理氣量下,保持入口截面尺寸大小不變,其它尺寸相似放大的旋風分離器分離效率隨筒徑的變化關系。將不同模型計算結果與實驗結果對比,可以看出模型預測結果和實驗結果隨筒徑變化的趨勢大體上一致,但二者之間還存在一定的差距。這是因為旋風分離器內部三維強旋湍流運動的復雜性,各理論模型均為研究者們基于各自的假設,通過對其工作過程的簡化而得,雖具有一定的準確性,但與實際分離過程相比,仍存在一定的局限性。

        圖7 旋風分離器總效率(η)的試驗值與模型計算值對比Fig.7 Comparison of overall efficiency(η) between experiments and model prediction(a) Overall efficiency with inlet violet 15 m/s (KA=7.85); (b) Overall efficiency with inlet violet 20 m/s (KA=7.85);(c) Overall efficiency with consumption flow 0.167 m3/s(KA=7.85);(d) Overall efficiency with consumption flow 0.167 m3/s (a×b=144×58 mm2)

        3.2 旋風分離器壓降模型的預測結果分析

        目前的旋風分離器壓降計算公式很多是經驗關聯(lián)式,且多用阻力系數(shù)ξ來表達。表2列出了幾種旋風分離器阻力系數(shù)的具體計算方法。

        在本研究中,選擇了7個經驗模型來預測旋風分離器的壓力損失,圖8為本實驗條件下不同筒體直徑的旋風分離器實驗壓降結果與模型壓降計算值的對比。從圖8(a)中可知,當KA=7.85且入口氣速相同時,袁怡等[21]模型計算值與試驗值有相同的變化趨勢且偏差較小,其余模型均忽視了筒徑對旋風分離器壓降的影響,總是保持一定值;從圖8(b)中可知,當KA=7.85且處理氣量相同時,各個模型計算值與試驗值隨筒徑變化的趨勢一致,即隨著筒體直徑的增加,旋風分離器壓降降低;從圖8(c)中可知,當入口高寬保持一定值時,除Barth等[18]模型外,其余模型預測結果與實驗結果變化趨勢一致,即隨著筒體直徑的增加,旋風分離器壓降降低。

        表2 旋風分離器阻力系數(shù)的計算模型Table 2 Calculation model of pressure drop coefficient for cyclones

        4 結 論

        以Stairmand型旋風分離器為基準模型,通過實驗測量和理論模型計算研究了筒體直徑對旋風分離器性能的影響。

        (1)分離器結構尺寸按幾何相似放大,入口氣速相同時,隨著筒體直徑的增大分離效率降低,壓降升高;處理氣量一定時,隨著筒體直徑的增大分離效率降低,壓降降低。

        (2)固定入口截面尺寸,其它結構尺寸按幾何相似放大,即處理氣量及入口氣速都相同時,隨著筒體直徑的增大分離效率升高,旋風分離器壓降降低。

        (3)模型計算與實驗值對比發(fā)現(xiàn),近年文獻中的幾個較新的分離模型基本上能夠預測旋風分離器分離效率和壓降隨筒徑的變化趨勢,但模型預測值都和實驗結果有一定差距。

        圖8 旋風分離器壓降(Δp)的試驗值與模型計算值對比Fig.8 Comparison of pressure drop(Δp) between experimental value and model prediction(a) Pressure drop with inlet violet 20 m/s (KA=7.85); (b) Pressure drop with consumption flow 0.167 m3/s (KA=7.85);(c) Pressure drop with consumption flow 0.167 m3/s (a×b=144×58 mm2)

        符號說明:

        a——矩形入口高度,m;

        ad——入口寬度修正系數(shù);

        a0r——平均流量寬度,m;

        b——矩形入口寬度,m;

        B——排塵口直徑,m;

        Cc——坎寧修正系數(shù);

        CD——曳力系數(shù);

        Cin——入口濃度,kg/m3;

        D——分離器筒體直徑,m;

        De——排氣管直徑,m;

        dp——顆粒直徑,m;

        Ds——摩擦面的無量綱平均直徑;

        d50——切割粒徑,m;

        D0——筒徑修正系數(shù),D0=0.3 m;

        f——含塵氣流與壁面摩擦系數(shù);

        f0——純凈氣流與壁面摩擦系數(shù);

        Fs——旋風分離器器壁的總面積,m2;

        g——重力加速度,m/s2;

        h——旋風分離器筒體高度,m;

        hc——旋風分離器錐體高度,m;

        H——旋風分離器總高度,m;

        KA——入口截面系數(shù),KA=πD2/4ab;

        Ks——絕對粗糙度,m;

        LVS——自然旋風長,m;

        n——旋流的旋渦指數(shù);

        Qin——入口氣體流量,m3/s;

        R——分離器筒體半徑,m;

        Re——排氣管半徑,m;

        Rin——入口半徑,m;

        Re——氣體雷諾數(shù);

        r0r——顆粒軌跡的平均半徑,m;

        Rw——分離器筒體半徑,m;

        rt——內外旋流分界半徑,m;

        S——排氣管插入深度,m;

        tres——總氣體停留時間,s;

        Uθw——分離器器壁顆粒切向速度分量,m/s;

        Ve——排氣管平均軸向氣速,m/s;

        Vin——入口氣體速度,m/s;

        V0——入口平均速度,m/s;

        Vθe——排氣管壁面氣體切向速度,m/s;

        Vθw——分離器器壁氣體切向速度,m/s;

        wi——第i個組分的質量分數(shù)

        希臘字母:

        ηi——粒級效率;

        η——總分離效率;

        φ——慣性碰撞系數(shù);

        ζ——進口損失;

        ζe——器內旋流損失;

        ρg——氣體密度,kg/m3;

        ρp——顆粒密度,kg/m3;

        ρstr——灰?guī)芏?,kg/m3;

        δ——顆粒軌道半徑與出口半徑之比;

        v——氣體運動黏度,m2/s;

        μg——氣體動力黏度,Pa·s;

        ξ——純氣流阻力系數(shù);

        ξst——純氣流阻力系數(shù)(Stairmand模型)。

        [1] 時銘顯, 吳小林. 旋風分離器的大型冷模試驗研究[J].化工機械, 1993, 20(4): 188-189. (SHI Mingxian, WU Xiaolin. Experimental study on large-scale cold mold cyclone[J].Chemical Machinery, 1993, 20(4): 188-189.)

        [2] 金國淼. 除塵設備[M].北京: 化學工業(yè)出版社, 2002: 40-41.

        [3] 金有海, 時銘顯. 旋風分離器相似放大實驗研究[J].石油大學學報, 1990, 14(5): 46-55. (JIN Youhai, SHI Mingxian. Experimental studies on scale-up of separator cyclone[J].Journal of the University of Petroleum, China, 1990, 14(5): 46-55.)

        [4] JIN W L, YANG H J, DONG Y L. Effect of the cylinder shape of a long-coned cyclone on the stable flow-field establishment[J].Powder Technology, 2006, 165(1): 30-38.

        [5] 袁惠新, 李麗麗. 物性和結構參數(shù)對旋風分離器壓降影響的數(shù)值模擬[J].礦山機械, 2010, 5(38): 91-94. (YUAN Huixin, LI Lili. Numerical simulation of material properties and structural parameters’ effects on pressure drop of the cyclone separator[J].Mining Machinery, 2010, 5(38): 91-94.)

        [6] 方麒先. 旋風分離器結構尺寸對壓降影響的實驗研究[D].北京: 中國石油大學(北京), 2015.

        [7] 劉書賢, 孫國剛, 時銘顯. 催化裂化裝置沉降器粗旋結構設計探討[J].煉油技術與工程, 2008, 10(38): 25-30. (LIU Shuxian, SUN Guogang, SHI Mingxian. Structure design of a rough-cut cyclone in fluidized catalytic cracking disengager[J].Petroleum Refinery Engineering, 2008, 10(38): 25-30.)

        [8] 宋健斐, 楊光福, 陳建義, 等. 旋風分離器內氣相流場的相似?;治?Ⅱ)尺寸參數(shù)[J].化工學報, 2010, 9(9): 2274-2279. (SONG Jianfei, YANG Guangfu, CHEN Jianyi, et al. Similarity analysis of modeling of gas phase flow field in cyclone separator(Ⅱ) Size parameters[J].CIESC Journal, 2010, 9(9): 2274-2279.)

        [9] OSCAR L, SGROTT J, DIRCEU N, et al. Cyclone optimization by COMPLEX method and CFD simulation[J].Powder Technology, 2015, 277(1): 11-21.

        [10] YANG J X, SUN G G, GAO C Z. Effect of the inlet dimensions on the maximum-efficiency cyclone height[J].Separation and Purification Technology, 2013, 105(1): 15-23.

        [11] ZHANG R, BASU P. A simple model for prediction of solid collection efficiency of a gas-solid separator[J].Powder Technology, 2004, 147(1): 86-93.

        [12] 錢付平, 張吉光, 楊豐, 等. 靜電旋風分離器結構尺寸的優(yōu)化設計[J].鍋爐技術, 2004, 35(2): 10-13. (QIAN Fuping, ZHANG Jiguang, YANG Feng, et al. Optimization design of configuration size on electro cyclone[J].Boiler Technology, 2004, 35(2): 10-13.)

        [13] ZHAO B T. Prediction of gas-particle separation efficiency for cyclones: A time-of-flight model[J].Separation and Purification Technology. 2012, 85(2): 171-177.

        [14] AVCI A, KARAGOZ I. Effects of flow and geometrical parameters on the collection efficiency in cyclone separators[J].Journal of Aerosol Science, 2003, 34(7): 937-955.

        [15] SHEPHERD C B, LAPPLE C E. Flow pattern and pressure drop in cyclone dust collectors cyclone without intel vane[J].Industrial & Engineering Chemistry, 1940, 32(9): 1246-1248.

        [16] DIRGO J. Relationships between cyclone dimensions and performance[D].Boston: Univ of Havarad, 1988.

        [17] 陳建義. 切流返轉式旋風分離器分離理論和優(yōu)化設計方法的研究[D].北京: 中國石油大學(北京), 2007.

        [18] BARTH W, LEINEWEBER L. Beurteilung und auslegung von zyklonabscheidern[J].Staub, 1964, 24(2): 41-55.

        [19] MUSCHELKNAUTZ E. Die berechnung von zyklonabscheidern für gase[J].Chemie Ingenieur Technik, 1972, 44(1-2): 63-71.

        [20] DEMIR S. A practical model for estimating pressure drop in cyclone separators: An experimental study[J].Powder Technology, 2014, 268: 329-338.

        [21] 袁怡, 孫國剛, 周發(fā)戚. 旋風分離器的壓力降計算[J].煉油技術與工程, 2016, 46(10): 39-44. (YUAN Yi, SUN Guogang, ZHOU Faqi. Pressure drop calculation of cyclone separator[J].Petroleum Refinery Engineering, 2016, 46(10): 39-44.)

        Effects of the Cylinder Diameter on Cyclone Performance

        YUAN Yi, SUN Guogang, ZHOU Faqi, SUN Zhanpeng

        (BeijingKeyLaboratoryofProcessFluidFiltrationandSeparation,CollegeofChemicalEngineering,ChinaUniversityofPetroleum,Beijing102249,China)

        Cylinder diameter was an important structural parameter of the cyclone. In this paper, a series of laboratory cyclone model experiments were carried out to investigate the effect of the cylinder diameter on Stairmand cyclone performances.The cylinder diameter was 200 mm, 300 mm and 400 mm, and the powder was microsilica. The results showed that, according to the geometric similarity to enlarge cyclone, if the inlet velocity was the same, with the increase of the cylinder diameters, the collection efficiency of the cyclone decreased and the pressure drop increased. Moreover, if gas flow rate was the same, with the increase of the cylinder diameters, the collection efficiency of the cyclone decreased and the pressure drop decreased. As the inlet size was fixed, if enlarging other dimensions with geometric similarity, namely inlet velocity and gas flow rate remained unchanged, with the increase of the cylinder diameters, the collection efficiency of the cyclone increased and the pressure drop decreased. In addition, under the experimental conditions, the comparison between the computations of several recent performance prediction models and the experimental date, showed that the theoretical models could predict the tendency of cylinder diameter effects on cyclone performances.

        cyclone; cylinder diameter; collection efficiency; pressure drop; model calculations

        2016-08-02

        國家自然科學基金(21276274)和國家重點基礎研究發(fā)展計劃“973”項目(2014CB744304)資助

        袁怡,女,碩士,從事氣-固分離的數(shù)值模擬及實驗研究;Tel:010-89734820;E-mail:m15201106422@163.com

        孫國剛,男,教授,博士,從事氣-固分離及流態(tài)化工程研究;Tel:010-89734820;E-mail:ggsunbj@163.com

        1001-8719(2017)04-0738-08

        TQ021.1

        A

        10.3969/j.issn.1001-8719.2017.04.019

        猜你喜歡
        旋風筒體分離器
        b型管板與筒體溫差應力的分析計算和評定
        化工管理(2021年7期)2021-05-13 00:46:04
        回轉窯筒體對接操作方法
        水泥技術(2021年2期)2021-04-20 12:37:26
        一種臥式筒體糞污發(fā)酵裝置的筒體設計與分析
        強制循環(huán)流化床物料分離效率的研究及應用
        旋風分離器切向進料接管焊縫修理與改進建議
        腦力急旋風
        腦力急旋風
        腦力急旋風
        腦力急旋風
        HAZOP理論分析在分離器作業(yè)中的應用探討
        安全(2015年1期)2016-01-19 06:19:07
        国产精选自拍视频网站| 大地资源在线播放观看mv| 亚洲日韩精品欧美一区二区 | 人人玩人人添人人澡| 亚洲国产日韩在线人成蜜芽| 官网A级毛片| 国产中文字幕亚洲综合| 日韩人妻大奶子生活片| 亚洲啪啪视频一区二区| 午夜爽爽爽男女污污污网站| 插b内射18免费视频| 18禁美女裸身无遮挡免费网站| 亚洲 欧美 激情 小说 另类| 亚洲国产精品嫩草影院久久av | 日本高清色一区二区三区| 亚洲男人的天堂av一区| 奇米影视7777久久精品| 中文字幕亚洲乱码熟女一区二区| 亚洲天堂免费视频| 精品无人区无码乱码大片国产| 伊人影院成人在线观看| 亚洲一二三区免费视频| 亚洲国产精品无码专区影院| 欧美黄色免费看| 亚洲中文字幕无码不卡电影| 羞涩色进入亚洲一区二区av| 成人大片免费视频播放一级| 国产精品人人做人人爽人人添| 色噜噜av亚洲色一区二区| 国产在线成人精品| 亚洲精品国产福利在线观看| 亚洲男人免费视频网站| 国产aⅴ无码专区亚洲av麻豆| 成人免费毛片内射美女-百度| 国产精品丝袜在线不卡| 伊人久久婷婷综合五月97色| 久久一区二区av毛片国产| 日韩午夜理论免费tv影院| 久久视频在线| 吸咬奶头狂揉60分钟视频| 九九精品国产99精品|