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

    旋風(fēng)分離器內(nèi)置導(dǎo)流葉片結(jié)構(gòu)參數(shù)優(yōu)化研究

    2024-12-31 00:00:00曾云陳宓魏軻李美求李偉偉董吉寧
    石油機(jī)械 2024年8期
    關(guān)鍵詞:效率優(yōu)化模型

    為研究旋風(fēng)分離器內(nèi)置導(dǎo)流葉片各個結(jié)構(gòu)參數(shù)之間的相互作用,基于CFD數(shù)值模擬方法和BBD試驗(yàn)設(shè)計(jì),采用二階多項(xiàng)式基函數(shù)建立了Stairmand旋風(fēng)分離器的葉片軸向位置、葉片長度及葉片偏轉(zhuǎn)角度與分離效率及壓降間的數(shù)學(xué)模型。擬合結(jié)果表明:響應(yīng)目標(biāo)的決定系數(shù)均在0.99以上,表明相關(guān)性和回歸模型效果較好。使用Design Expert軟件處理數(shù)據(jù)后選擇了最佳推薦點(diǎn),對分離器結(jié)構(gòu)進(jìn)行了優(yōu)化。優(yōu)化后的模型與原結(jié)構(gòu)對比結(jié)果如下:優(yōu)化模型的壓降略有上升;顆粒直徑在1~20 μm范圍內(nèi)時,優(yōu)化前后分離器的分離效率變化明顯,尤其在1~15 μm范圍內(nèi)時,優(yōu)化前分離效率為6.8%~31.0%,優(yōu)化后則達(dá)到了17.0%~54.0%,提升約20個百分點(diǎn);而當(dāng)粒徑在15~40 μm區(qū)間時,優(yōu)化前后的分離效率變化不明顯,尤其是在粒徑20~40" μm時,優(yōu)化前的分離效率為71.0%~98.0%,而優(yōu)化后的分離效率只有77%~99%,提升約10個百分點(diǎn);優(yōu)化后的模型性能符合“高效低阻”的設(shè)計(jì)目標(biāo)。所得結(jié)論可為旋風(fēng)分離器內(nèi)置葉片的改進(jìn)提供設(shè)計(jì)指導(dǎo)。

    旋風(fēng)分離器;Fluent;結(jié)構(gòu)優(yōu)化;響應(yīng)面法;BBD試驗(yàn)設(shè)計(jì);分離效率

    TE926

    A

    012

    Optimization on Structural Parameters of Guide

    Vanes Built in Cyclone Separators

    Zeng Yun1 Chen Mi1 Wei Ke2 Li Meiqiu1 Li Weiwei3 Dong Jining3

    (1.School of Mechanical Engineering,Yangtze University;2.Shandong Natural Gas Pipeline Co.,Ltd.;3.Shandong Weiyuan Technology Co.,Ltd.)

    In order to study the interaction among structural parameters of the guide vanes built in cyclone separator,based on CFD numerical simulation method and BBD test design,a second-order polynomial basis function was used to build a mathematical model for the relationships among the axial position,length and deflection angle of the blade,the separation efficiency and the pressure drop of the Stairmand cyclone separator.The fitting results show that the determination coefficients of the response target are all above 0.99,indicating a good effect of correlation and regression model.Moreover,after having processed the data using Design Expert software,the best recommendation point was selected to conduct optimization on the structure of the separator.The comparison between the optimized model and the original structure reveals that the pressure drop of the optimized model slightly increases.When the particle diameter is within the range of 1-20 μm,there is a significant change in the separation efficiency of the separator before and after optimization,especially within the range of 1-15 μm,the separation efficiency is 6.8%-31.0% before optimization,and reaches 17.0%-54.0% after optimization,an improvement of about 20 percentage points.When the particle diameter is within the range of 15-40 μm,there is no significant change in separation efficiency before and after optimization,especially within the range of 20-40 μm,the separation efficiency is 71.0%-98.0% before optimization,and reaches 77%-99% after optimization,an improvement of about 10 percentage points.The performance of the optimized model conforms to the design goal of “high efficiency and low resistance”.The study conclusions provide design guidance for the improvement of the guide vanes built in the cyclone separator.

    cyclone separator;Fluent;structural optimization;response surface method;BBD test design;separation efficiency

    0 引 言

    旋風(fēng)分離器作為鉆井平臺固控系統(tǒng)中的重要設(shè)備,利用氣流的旋轉(zhuǎn)產(chǎn)生離心力來分離氣流中的顆粒。與濾膜法、重力沉降等分離設(shè)備相比,旋風(fēng)分離器因其性能穩(wěn)定、分離效率高、結(jié)構(gòu)簡單、操作簡便等優(yōu)點(diǎn)被廣泛應(yīng)用于環(huán)保及油氣處理等領(lǐng)域[1-2。長久以來,國內(nèi)外學(xué)者致力于旋風(fēng)分離器的性能研究,其中以經(jīng)典的Stairmand旋風(fēng)分離器為主要研究對象。Stairmand分離器內(nèi)部主要分為外旋流和內(nèi)旋流2部分,除此之外還有短路流、偏心環(huán)流等二次流現(xiàn)象,內(nèi)部流場較為復(fù)雜,對分離效率的影響較大。針對分離器結(jié)構(gòu)參數(shù)與操作參數(shù)的優(yōu)化研究較多,具體有:封躍鵬等[3將傳統(tǒng)出口結(jié)構(gòu)設(shè)計(jì)成帶夾層的結(jié)構(gòu),使從入口進(jìn)入的氣流加速旋轉(zhuǎn),改善了流場狀況,提高了分離效率;解凱等4基于標(biāo)準(zhǔn)Stairmand型旋風(fēng)分離器,設(shè)計(jì)了2種外導(dǎo)流管式旋風(fēng)分離器,優(yōu)化了分離器內(nèi)二次渦的分布;趙新學(xué)等[5基于CFD模擬研究了旋風(fēng)分離器顆粒出口直徑對其壁面磨損的影響,研究發(fā)現(xiàn),隨著顆粒出口直徑的減小,分離空間壁面的磨損增加,在圓錐底部磨損增加程度最大;付烜等[6采用數(shù)值模擬方法,研究不同入口形式對旋風(fēng)分離器分離效率的影響,研究發(fā)現(xiàn),分離器在環(huán)形管入口處的分離效率略高于兩側(cè)入口的分離效率,壓降損失降低了25%;凌光磊等[7探究了旋風(fēng)分離器蝸殼半徑對分離性能的影響,研究發(fā)現(xiàn),隨著蝸殼半徑增大,分離效率逐漸提高,并且蝸殼半徑增大,物料分離性能越好;趙洋等8針對分離器入口面積對旋風(fēng)分離器性能的影響進(jìn)行了研究,發(fā)現(xiàn)在進(jìn)氣流量減小后,入口面積減小有利于提高分離器的分離效率。除了傳統(tǒng)的優(yōu)化結(jié)構(gòu)參數(shù)與操作參數(shù)外,學(xué)者們對分離器內(nèi)部構(gòu)件的研究也越來越重視,并提出了創(chuàng)新性的結(jié)構(gòu):楊昌智等[9在旋風(fēng)分離器內(nèi)部添加旋流葉片,研究發(fā)現(xiàn),在入口速度較小時分離效率顯著提高,在入口速度較大時分割粒徑顯著減??;ZHANG S.等10在分離器排塵口附近添加防混反錐,驗(yàn)證了防混反錐與排塵口的最佳距離為0.15倍的筒體直徑;梁嘉豪等[11則在分離器進(jìn)口處設(shè)置導(dǎo)流板,研究表明,分離器的分離效率明顯提高。與此同時,國內(nèi)外許多學(xué)者針對分離器內(nèi)部構(gòu)件的不同結(jié)構(gòu)參數(shù)進(jìn)行了優(yōu)化研究。但大多數(shù)都是僅針對單一因素進(jìn)行研究,未將其他影響因素綜合考慮,在實(shí)際工程意義上仍有欠缺。

    響應(yīng)面法是一種將試驗(yàn)設(shè)計(jì)與數(shù)學(xué)建模相結(jié)合的試驗(yàn)設(shè)計(jì)方法,通過在局部范圍內(nèi)的代表性試驗(yàn)點(diǎn)上進(jìn)行試驗(yàn),可以全局?jǐn)M合設(shè)計(jì)變量和響應(yīng)目標(biāo)之間的函數(shù)關(guān)系[12。并且,響應(yīng)面法考慮了試驗(yàn)的隨機(jī)誤差,將復(fù)雜未知的函數(shù)關(guān)系在一定的試驗(yàn)區(qū)域內(nèi)用簡單的多項(xiàng)式模型來擬合,計(jì)算簡便13-14。本文在已有的內(nèi)置導(dǎo)流葉片的旋風(fēng)分離器基礎(chǔ)上,采用響應(yīng)面法建立了葉片軸向位置、葉片長度及葉片偏轉(zhuǎn)角度與壓降和分離效率之間的數(shù)學(xué)模型,基于模型對不同參數(shù)之間的耦合進(jìn)行了研究,進(jìn)而對旋風(fēng)分離器進(jìn)行了結(jié)構(gòu)優(yōu)化,以期為內(nèi)置導(dǎo)流葉片的旋風(fēng)分離器設(shè)計(jì)提供理論指導(dǎo)。

    1 旋風(fēng)分離器內(nèi)部結(jié)構(gòu)及工作原理

    在經(jīng)典旋風(fēng)分離器中,氣流在其內(nèi)部被分為速度較快的外旋流和速度較慢的內(nèi)旋流。分離器在分離過程中,會出現(xiàn)部分粗顆粒被卷入速度較慢的內(nèi)旋流的情況。當(dāng)部分粗顆粒被卷入內(nèi)旋流后,較低的切向速度無法提供足夠的離心力使顆粒擺脫向上的螺旋氣流,從而隨著細(xì)顆粒被帶出排氣口,造成不充分分離。同時,在圓錐體的下端存在偏心環(huán)流的現(xiàn)象,這會導(dǎo)致已被分離至排塵口的粗顆粒重新被卷入內(nèi)旋流。

    為了阻擋部分粗顆粒被卷入內(nèi)旋流中,并使顆粒在分離器內(nèi)部的停留時間延長,同時使其內(nèi)部的氣壓分布具有更好的對稱性,削弱底部偏心環(huán)流和旋渦擺尾對分離效率的影響,這里對Stairmand旋風(fēng)分離器添加了內(nèi)置旋流葉片,將導(dǎo)流葉片置于內(nèi)外旋交界面之間(54.87~67.30 mm),距溢流管下端50 mm處?;跅畈堑?sup>[9提出的內(nèi)置導(dǎo)流葉片深入研究葉片參數(shù)的選取。同時,考慮到過大的葉片偏轉(zhuǎn)角度和葉片厚度會使壓降大幅增大,此處初選葉片偏轉(zhuǎn)角度為10°,葉片厚度為2 mm。葉片數(shù)量會影響分離器對氣體的導(dǎo)流效果,初選葉片數(shù)量為4個。具體生成過程及尺寸如圖1所示。

    2 數(shù)值模擬

    2.1 連續(xù)相方程

    采用CFD方法研究旋風(fēng)分離器的分離特性。數(shù)值模擬所采用的連續(xù)性方程如下:

    N-S方程為:

    式中:u、u′分別為流體時均速度和脈動速度,m/s;ui和uj為速度分量,m/s;ρ為流體密度,kg/m3;p為壓力,Pa;xi和xj為坐標(biāo)分量,m;μ為動力黏度,Pa·s;g為重力加速度,m/s2;u′iu′j為雷諾應(yīng)力分量。

    動量守恒方程如下:

    式中:u、v和w分別為流體在t時刻在點(diǎn)(x、y、z)處的速度分量,m/s;fx、fy、fz為控制體各方向上的質(zhì)量力,m/s2

    2.2 湍流方程

    由于雷諾應(yīng)力在湍流流動中具有明顯的各向異性,所以在旋風(fēng)分離器的數(shù)值模擬中采用雷諾應(yīng)力模型(RSM)。該模型對旋風(fēng)分離器具有強(qiáng)旋轉(zhuǎn)運(yùn)動的問題研究具有較高精度。雷諾應(yīng)力輸運(yùn)方程如下:

    式中:Di,j為擴(kuò)散項(xiàng),Pa/s;pi,j為應(yīng)力產(chǎn)生項(xiàng),Pa/s;Gi,j為浮力產(chǎn)生項(xiàng),Pa/s;φi,j為壓力應(yīng)變再分配項(xiàng),Pa/s;εi,j為離散項(xiàng),Pa/s;Fi,j為旋轉(zhuǎn)系統(tǒng)產(chǎn)生項(xiàng),Pa/s;Suser為自定義源項(xiàng),Pa/s。

    2.3 多相流模型

    對于多相流問題,主要采用歐拉-歐拉方法和歐拉-拉格朗日方法,考慮所研究的旋風(fēng)分離器中固相體積分?jǐn)?shù)不超過10%,因此選用歐拉-拉格朗日方法,也就是離散相模型(DPM)模擬。其顆粒運(yùn)動時的平衡方程為:

    式中:mp為顆粒質(zhì)量,kg;FD 為顆粒所受曳力,N;up為顆粒速度,m/s;ρp為顆粒密度,kg/m3;dp為顆粒直徑,m;Re為相對雷諾數(shù);F為顆粒所受到的其他作用力,包括虛擬質(zhì)量力、薩夫曼升力、布朗力等,N。

    在這里,粉煤灰為亞觀尺寸顆粒,因此需要考慮薩夫曼升力,而ρlt;ρp,虛擬質(zhì)量力可忽略。同時,由于DPM模型不考慮顆粒之間的相互作用力,顆粒所受的布朗力也可忽略。

    2.4 參數(shù)設(shè)置及邊界條件

    將旋風(fēng)分離器的入口邊界條件設(shè)置為速度入口,出口邊界條件設(shè)置為壓力出口。入射的流體為空氣。入口風(fēng)速設(shè)置為15 m/s,湍流強(qiáng)度為4%,入射顆粒為體積分?jǐn)?shù)小于5%的粉煤灰。由于顆粒對氣體流動的影響可以不計(jì),所以采用基于歐拉-拉格朗日方法的DPM離散相模型模擬氣固兩相流動。旋風(fēng)分離器內(nèi)部流動較復(fù)雜,因此采用湍流脈動的隨機(jī)軌道模型[15。入射顆粒物粉煤灰的密度為1 850 kg/m3,顆粒的初始速度設(shè)置為與氣體速度相同為15 m/s。

    2.5 網(wǎng)格劃分與無關(guān)性驗(yàn)證

    由于內(nèi)置導(dǎo)流葉片處流域較為復(fù)雜,考慮網(wǎng)格劃分所需要的時間和模擬的計(jì)算精度,采用適應(yīng)性較好的非結(jié)構(gòu)化四面體網(wǎng)格進(jìn)行網(wǎng)格劃分。在網(wǎng)格數(shù)為53萬、58萬、60萬和68萬下進(jìn)行模擬,并選取某截面處切向速度進(jìn)行對比,以此進(jìn)行網(wǎng)格無關(guān)性驗(yàn)證,結(jié)果如圖2所示。

    由圖2可以看出:網(wǎng)格數(shù)量對計(jì)算結(jié)果有顯著影響,當(dāng)網(wǎng)格數(shù)為53萬時,切向速度的分布明顯較平緩;當(dāng)網(wǎng)格數(shù)為58萬時,切向速度分布在靠近中心位置,較前者已經(jīng)有較大差異;當(dāng)網(wǎng)格數(shù)為60萬時,計(jì)算結(jié)果與網(wǎng)格數(shù)量68萬的結(jié)果無太大差異。因此,可選擇60萬網(wǎng)格數(shù)作為離散標(biāo)準(zhǔn)。3 結(jié)構(gòu)參數(shù)響應(yīng)面優(yōu)化

    3.1 單因素參數(shù)選擇

    導(dǎo)流葉片上端到溢流管下端之間的距離為導(dǎo)流葉片的軸向位置,考慮到葉片軸向位置過小會使溢流管附近的“節(jié)流效應(yīng)”增強(qiáng),選定軸向位置d=50~250 mm。葉片偏轉(zhuǎn)角度過大會使壓降增大并增加粗顆粒進(jìn)入內(nèi)旋流的概率,選定葉片偏轉(zhuǎn)角度α=10°~15°。旋風(fēng)分離器的溢流管下端距離排塵口為670 mm,用字母E來表示。選定葉片長度y=0.4E~0.6E,即268~402 mm之間。在不改變其他結(jié)構(gòu)參數(shù)和工況的前提下,通過研究單因素參數(shù)對分離器分離性能的影響,進(jìn)而選取合適的參數(shù)范圍進(jìn)行響應(yīng)面優(yōu)化。

    圖3為不同軸向位置時的壓力云圖。由圖3可知,在不同軸向位置下,中心空氣柱均呈現(xiàn)較好的柱狀分布,外部輪廓明顯且無明顯彎曲。在緊靠溢流管的下端產(chǎn)生小范圍的收縮現(xiàn)象,這是因?yàn)樯闲械臍饬鲝膶?dǎo)流葉片上部出口流出時,管徑突然增大使此處的靜壓降低,從而導(dǎo)致氣體被壓縮。

    圖4為不同葉片長度時的壓力云圖。由圖4可以看出:隨葉片長度的增大,分離器中心處負(fù)壓連續(xù)性逐漸下降;在葉片長度為268~301 mm時,中心處的負(fù)壓區(qū)均呈現(xiàn)出良好的連續(xù)性;當(dāng)葉片長度為335 mm時,負(fù)壓連續(xù)性變差,且最大負(fù)壓出現(xiàn)在溢流管下端區(qū)域。這是因?yàn)楫?dāng)導(dǎo)流葉片距離溢流管下端較大時,分離器內(nèi)的壓力差變大、湍流強(qiáng)度增強(qiáng),在排塵口附近產(chǎn)生明顯擾動。

    圖5為不同偏轉(zhuǎn)角度時的壓力云圖。由圖5可以看出:當(dāng)偏轉(zhuǎn)角度為13°時,最大靜壓值最小,且內(nèi)旋區(qū)和外旋區(qū)的壓差明顯,靜壓值從兩側(cè)中心處遞減;當(dāng)偏轉(zhuǎn)角度在10°~15°的范圍內(nèi),負(fù)壓值隨角度增大而增大。這是因?yàn)閷?dǎo)流葉片的存在,兩葉片之間的縫隙處空間狹小,空氣氣流增大,從而使得靜壓值降低。

    通過在局部范圍內(nèi)對具有代表性的點(diǎn)進(jìn)行試驗(yàn),分別選取葉片長度為268~301 mm,葉片軸向位置為100~250 mm,偏轉(zhuǎn)角度為12°~14°。

    3.2 BBD試驗(yàn)設(shè)計(jì)

    響應(yīng)面優(yōu)化方法是使用簡單的多項(xiàng)式模型擬合某個試驗(yàn)區(qū)域內(nèi)的復(fù)雜未知函數(shù)關(guān)系。在試驗(yàn)條件范圍內(nèi)分析時,可以對不同水平的不同因素進(jìn)行連續(xù)分析,從而獲得準(zhǔn)確的最優(yōu)解[16。BBD設(shè)計(jì)在每個因素中取3個水平值,分別為0、-1、1,以(0,0,0)為中心點(diǎn)進(jìn)行編碼,-1和1分別代表每個因素的最大值和最小值。由于BBD設(shè)計(jì)的試驗(yàn)次數(shù)較少,且優(yōu)化得出的參數(shù)值不會超過參數(shù)的取值范圍,所以采用BBD方法進(jìn)行試驗(yàn)樣本點(diǎn)的設(shè)計(jì)。這里對葉片軸向位置、葉片長度、葉片偏轉(zhuǎn)角度3個結(jié)構(gòu)參數(shù)進(jìn)行優(yōu)化,優(yōu)化參數(shù)及范圍如表1所示。

    采用Design-Expert軟件進(jìn)行試驗(yàn)點(diǎn)設(shè)計(jì),共進(jìn)行15組數(shù)值模擬試驗(yàn),如表2所示。設(shè)計(jì)變量分別為葉片軸向位置P1、葉片偏轉(zhuǎn)角度P2及葉片長度P3,目標(biāo)變量為壓降P4和分離效率P5。

    通過這15組試驗(yàn)可以看出:當(dāng)設(shè)計(jì)變量同時發(fā)生變化時,分離效率最低值為64.94%,最高值為82.42%;而壓降最低值為719.5 Pa,最高值為763.8 Pa。由此可見,在邊界條件不變的情況下,分離效率和壓降均發(fā)生明顯變化,原因可能在于結(jié)構(gòu)參數(shù)對分離效率和壓降存在明顯影響。

    3.3 基于kriging模型的響應(yīng)面擬合分析

    根據(jù)上述試驗(yàn)結(jié)果,采用kriging建立模型,其表達(dá)形式如下:

    式中:x為維度為d的設(shè)計(jì)參數(shù)[x1,x2,…,xd];y(x)為多項(xiàng)式模型;z(x)為正態(tài)分布函數(shù),其不獨(dú)立但同分布;f(x)為類似響應(yīng)面模型;m為多項(xiàng)式的數(shù)目。

    將表2中的數(shù)據(jù)導(dǎo)入Design-expert軟件中,建立響應(yīng)目標(biāo)(分離效率和壓降)分別與P1(葉片軸向位置)、P2(葉片偏轉(zhuǎn)角度)及P3(葉片長度)的回歸方程,并對以分離效率和壓降為基礎(chǔ)建立的kriging模型進(jìn)行擬合分析[17。

    為了驗(yàn)證回歸方程的預(yù)測精度,對其進(jìn)行誤差統(tǒng)計(jì)分析,得到如表3所示的分析結(jié)果。表3中Y1為壓降,Y2為分離效率。

    由表3可知,Y1和Y2的決定系數(shù)均在0.99以上,并且校正決定系數(shù)和預(yù)測決定系數(shù)這2個值高且接近(校正決定系數(shù)-預(yù)測決定系數(shù)<2)。這表明相關(guān)性和回歸模型效果較好,擬合得到的回歸方程可以較好地滿足誤差的要求[18,預(yù)測的內(nèi)置導(dǎo)流葉片的旋風(fēng)分離器壓力損失模型與分離效率模型合理。

    3.4 響應(yīng)面分析

    為了直觀地分析幾何參數(shù)對性能目標(biāo)的影響,通過將其他參數(shù)設(shè)置為常量,將所研究的2個設(shè)計(jì)變量作為因變量,分析兩兩參數(shù)的交互作用以獲得最佳參數(shù)[19。這里根據(jù)回歸方程進(jìn)行葉片軸向位置、葉片偏轉(zhuǎn)角度、葉片長度3個因素的響應(yīng)曲面分析。

    3.4.1 單因素對響應(yīng)目標(biāo)的影響

    圖6為各單因素設(shè)計(jì)變量對目標(biāo)變量的影響。由圖6可以看出,分離器分離效率與壓降在各單因素變量之間均呈先上升后下降的趨勢,并且在達(dá)到頂端后,下降速度逐漸加快。如圖6a所示,分離效率及壓降均在P1為153.57 mm,即葉片距離溢流管底端為153.57 mm時達(dá)到最高,之后便呈下降趨勢;如圖6b所示,分離效率在P2(葉片偏轉(zhuǎn)角度)為13.14°時達(dá)到最大,而壓降則在P2為12.57°時達(dá)到頂峰,但總體來說P2對壓降的影響較微弱;如圖6c所示,分離器分離效率在P3為286.28 mm時達(dá)到最高,而壓降則在P3為281.71 mm時達(dá)到頂峰。

    3.4.2 雙因素變量對響應(yīng)目標(biāo)的影響

    圖7為不同雙變量對分離器壓降的影響規(guī)律。從圖7a、圖7c可以看出,P4隨著P1、P2和P3的增大呈先上升后下降的趨勢,但相比于P1和P3,P2對P4的影響較微弱,這與單因素分析中得出的結(jié)論一致。而從圖6b可以看出,隨著P1和P3不斷增加,壓降P4呈明顯的先上升后下降趨勢,這同樣與單因素分析中得出的結(jié)論一致。

    圖8為不同雙變量對分離器分離效率的影響規(guī)律。從圖8a可以看出,當(dāng)P1不變時,隨著P2的增加,分離效率先升高后下降。當(dāng)葉片軸向位置增大,分離效率隨葉片偏轉(zhuǎn)角度增加的變化趨勢由快速上升后趨于平穩(wěn)逐漸變?yōu)槠骄徤仙罂焖傧陆?。從圖8b、圖8c可以看出,當(dāng)P1、P2保持不變時,隨著P3的增加,分離器分離效率逐漸升高達(dá)到頂峰后又呈下降趨勢,但相比于P1和P2,P3對P5的影響明顯較弱,這與前文單因素影響規(guī)律的結(jié)論一致。

    3.5 優(yōu)化結(jié)果及驗(yàn)證

    以“高效低阻”為設(shè)計(jì)目標(biāo)[20,需要考慮在滿足提高分離效率的同時保持壓降不大幅度升高。將分離效率和壓降的優(yōu)化比重均設(shè)置為1,以最大分離效率和最小壓降為優(yōu)化目標(biāo)。通過Design-expert軟件對數(shù)據(jù)進(jìn)行處理,在Solution中選擇最佳推薦點(diǎn)。經(jīng)過優(yōu)化后的設(shè)計(jì)變量取值為:P1(葉片軸向位置)=100 mm,P2(葉片偏轉(zhuǎn)角度)=14°,P3(葉片長度)=300 mm。優(yōu)化前、后旋風(fēng)分離器內(nèi)置導(dǎo)流葉片的結(jié)構(gòu)參數(shù)如表4所示。

    為進(jìn)一步驗(yàn)證優(yōu)化后的性能得到提升,根據(jù)優(yōu)化所得設(shè)計(jì)變量的參數(shù)值進(jìn)行建模。在邊界條件不變的情況下,對比分析優(yōu)化前后的分離效率變化,結(jié)果如圖9所示。由圖9a可以看出,在粒徑為1~20 μm內(nèi)時,優(yōu)化前、后分離器的分離效率變化明顯,尤其在1~15 μm范圍內(nèi)時,優(yōu)化前分離效率為6.8%~31.0%,優(yōu)化后則達(dá)到了17.0%~54.0%,提升約20個百分點(diǎn)。而當(dāng)粒徑在15~40 μm內(nèi)時,優(yōu)化前、后的分離效率變化不明顯,尤其是在20~40 μm時,優(yōu)化前的分離效率為71%~98%,而優(yōu)化后的分離效率只有77%~99%,提升約10個百分點(diǎn)。

    由圖9b可見,雖然內(nèi)置導(dǎo)流葉片的分離器壓降相比經(jīng)典Stairmand分離器略有上升,但相比于壓降升高,分離效率有更大幅度的提高,優(yōu)化所得參數(shù)仍符合“高效低阻”設(shè)計(jì)目標(biāo)。

    4 結(jié)論及建議

    (1)使用BBD方法設(shè)計(jì)進(jìn)行了15組試驗(yàn),并驗(yàn)證了使用kriging模型建立的分離效率和壓降響應(yīng)面函數(shù)的擬合精度。結(jié)果表明,響應(yīng)目標(biāo)的相關(guān)系數(shù)均在0.99以上,并且回歸模型效果較好,回歸方程符合精度要求。

    (2)內(nèi)置導(dǎo)流葉片的旋風(fēng)分離器內(nèi)部流場相比經(jīng)典的Stairmand旋風(fēng)分離器更加穩(wěn)定,并且阻擋了部分粗顆粒被卷入內(nèi)旋流中,使顆粒在旋風(fēng)分離器內(nèi)部停留的時間延長,提高了分離效率。

    (3)對建立的響應(yīng)面函數(shù)進(jìn)行了響應(yīng)面分析,結(jié)果表明,葉片軸向位置和葉片偏轉(zhuǎn)角度對分離效率的影響更加明顯,但葉片偏轉(zhuǎn)角度對壓降的影響較小,并且分離效率和壓降的變化趨勢均隨3個結(jié)構(gòu)參數(shù)的增大而先增大后減小。

    (4)通過Design-expert軟件優(yōu)化,得到了旋流葉片的最佳結(jié)構(gòu)參數(shù),并進(jìn)一步驗(yàn)證得出:優(yōu)化后分離器的分離效率在粒徑為1~15 μm時得到明顯提高,相比于優(yōu)化前提高了約20個百分點(diǎn)。同時壓降也有一定程度的升高,但相比于壓降升高,分離效率有更大幅度的提高,優(yōu)化所得參數(shù)仍符合“高效低阻”設(shè)計(jì)目標(biāo)。

    (5)建議后續(xù)開展旋風(fēng)分離器結(jié)構(gòu)參數(shù)之間的耦合仿真,以完善研究結(jié)果,為優(yōu)化旋風(fēng)分離器奠定理論基礎(chǔ)、提供數(shù)據(jù)支持。

    [1] 梁延璽,孫寶福,吳健懋.入口位置對旋風(fēng)分離器性能的影響[J].輕金屬,2023(5):15-20.

    LIANG Y X,SUN B F,WU J M.Effect of inlet position on the performance of the cyclone separator[J].Light Metals,2023(5): 15-20.

    [2] KUMAR A,DHAWAN S,KUMAR M V,et al.Detection and identification of shape,size,and concentration of particulate matter in ambient air using bright field microscopy-based system[J].Atmospheric Pollution Research,2023,14(11): 101913.

    [3] 封躍鵬,王宇虹,姜大志.一種新型的旋風(fēng)分離器[J].中國粉體技術(shù),2008,14(4):43-46.

    FENG Y P,WANG Y H,JIANG D Z.A new-style cyclone separator[J].China Powder Science and Technology,2008,14(4): 43-46.

    [4] 解凱,王娟,鄒槊,等.外導(dǎo)流管對旋風(fēng)分離器流場的調(diào)控[J].石油學(xué)報(石油加工),2022,38(3):571-581.

    XIE K,WANG J,ZOU S,et al.Cyclone separator flow field control through outer diversion tube[J].Acta Petrolei Sinica (Petroleum Processing Section),2022,38(3): 571-581.

    [5] 趙新學(xué),金有海.排塵口直徑對旋風(fēng)分離器壁面磨損影響的數(shù)值模擬[J].機(jī)械工程學(xué)報,2012,48(6):142-148.

    ZHAO X X,JIN Y H.Effect of dust discharge diameter on wall erosion in cyclone separator[J].Journal of Mechanical Engineering,2012,48(6): 142-148.

    [6] 付烜,孫國剛,劉佳,等.旋風(fēng)分離器進(jìn)口渦旋感生速度場的減阻增效作用[J].化工學(xué)報,2011,62(7):1927-1932.

    FU X,SUN G G,LIU J,et al.Effect of induced velocity on separation efficiency and pressure drop of cyclones caused by vortex in vortex-tube inlet pipe[J].CIESC Journal,2011,62(7): 1927-1932.

    [7] 凌光磊,劉云峰,黃濤,等.旋風(fēng)分離器蝸殼半徑對物料分離性能的影響[J].礦山機(jī)械,2023,51(1):45-51.

    LING G L,LIU Y F,HUANG T,et al.Influence of radius of volute of cyclone separator on material separation performance[J].Mining amp; Processing Equipment,2023,51(1): 45-51.

    [8] 趙洋,陳建義,曹鳴謙,等.入口面積可變式旋風(fēng)分離器的性能[J].石油學(xué)報(石油加工),2021,37(5):1031-1039.

    ZHAO Y,CHEN J Y,CAO M Q,et al.Performance of the cyclone separator with variable inlet area[J].Acta Petrolei Sinica (Petroleum Processing Section),2021,37(5): 1031-1039.

    [9] 楊昌智,彭紅,商宇軒,等.一種內(nèi)置旋流葉片的新型旋風(fēng)分離器[J].流體機(jī)械,2019,47(12):33-38.

    YANG C Z,PENG H,SHANG Y X,et al.A new cyclone separator with built-in swirling blades[J].Fluid Machinery,2019,47(12): 33-38.

    [10] ZHANG S,SHIN W G.The effect of a counter cone on the performance of an axial cyclone separator[J].Journal of Mechanical Science and Technology,2023,37(9): 4889-4898.

    [11] 梁嘉豪,李陽,孔祥功.導(dǎo)流板對旋風(fēng)分離器內(nèi)氣固兩相分離性能的影響[J].過程工程學(xué)報,2014,14(1):36-41.

    LIANG J H,LI Y,KONG X G.Effect of deflector on gas-solid separation performance in cyclone separators[J].The Chinese Journal of Process Engineering,2014,14(1): 36-41.

    [12] 李莉,張賽,何強(qiáng),等.響應(yīng)面法在試驗(yàn)設(shè)計(jì)與優(yōu)化中的應(yīng)用[J].實(shí)驗(yàn)室研究與探索,2015,34(8):41-45.

    LI L,ZHANG S,HE Q,et al.Application of response surface methodology in experiment design and optimization[J].Research and Exploration in Laboratory,2015,34(8): 41-45.

    [13] 邢雷,李金煜,趙立新,等.基于響應(yīng)面法的井下旋流分離器結(jié)構(gòu)優(yōu)化[J].中國機(jī)械工程,2021,32(15):1818-1826.

    XING L,LI J Y,ZHAO L X,et al.Structural optimization of downhole hydrocyclones based on response surface methodology[J].China Mechanical Engineering,2021,32(15): 1818-1826.

    [14] 熊攀,鄢曙光,劉瑋寅.基于響應(yīng)曲面法的旋風(fēng)分離器結(jié)構(gòu)優(yōu)化[J].化工學(xué)報,2019,70(1):154-160.

    XIONG P,YAN S G,LIU W Y.Structure optimization of cyclone based on response surface method[J].CIESC Jorunal,2019,70(1): 154-160.

    [15] 周培源,林建忠,魏中磊.二維自由剪切湍流的實(shí)驗(yàn)結(jié)果及其與理論計(jì)算的比較[J].水動力學(xué)研究與進(jìn)展,1991(2):1-9.

    ZHOU P Y,LIN J Z,WEI Z L.the experimental results of 2-D free shear turbulent flows and the comparison with the theoretical computation[J].Journal of Hydrodynamics,1991(2): 1-9.

    [16] 熊攀.基于響應(yīng)曲面法的旋風(fēng)分離器結(jié)構(gòu)優(yōu)化研究[D].武漢:武漢科技大學(xué),2019.

    XIONG P.Structure optimization of cyclone based on response surface method[D].Wuhan: Wuhan University of Science and Technology,2019.

    [17] 李亞寧,郝亞超,李亮,等.響應(yīng)曲面法優(yōu)化納濾去除廢水中重金屬離子的研究[J].膜科學(xué)與技術(shù),2020,40(6):111-117.

    LI Y N,HAO Y C,LI L,et al.Optimization of removal heavy metal ions from high-salt wastewater by nanofiltration based on response surface methodology[J].Membrane Science and Technology,2020,40(6): 111-117.

    [18] 李貝,張倩.基于響應(yīng)面法的挖掘機(jī)動臂優(yōu)化設(shè)計(jì)[J].工程機(jī)械,2021,52(5):79-85.

    LI B,ZHANG Q.Optimization design of excavator boom based on response surface method[J].Construction Machinery and Equipment,2021,52(5): 79-85.

    [19] 張煥梅,楊瑞剛.基于響應(yīng)面法的起重機(jī)結(jié)構(gòu)可靠性靈敏度分析[J].中國工程機(jī)械學(xué)報,2020,18(2):131-136.

    ZHANG H M,YANG R G.Reliability sensitivity analysis of crane structure based on response surface method[J].Chinese Journal of Construction Machinery,2020,18(2): 131-136.

    [20] 趙盼婷.超高壓除砂器旋流筒結(jié)構(gòu)對除砂性能的影響[D].西安:西安石油大學(xué),2020.

    ZHAO P T.Effect of the structure of the super-high pressure sand remover cyclone on the performance of sand removal[D].Xi’an: Xi’an Shiyou University,2020.第一

    曾云,講師,生于1991年,2020年畢業(yè)于長江大學(xué)石油礦場機(jī)械專業(yè),獲博士學(xué)位,現(xiàn)從事壓裂設(shè)備方面的科研和教學(xué)工作。地址:(434023)湖北省荊州市。email:mechanicszy@163.com。

    2024-03-25

    王剛慶

    猜你喜歡
    效率優(yōu)化模型
    一半模型
    超限高層建筑結(jié)構(gòu)設(shè)計(jì)與優(yōu)化思考
    民用建筑防煙排煙設(shè)計(jì)優(yōu)化探討
    關(guān)于優(yōu)化消防安全告知承諾的一些思考
    一道優(yōu)化題的幾何解法
    重要模型『一線三等角』
    提升朗讀教學(xué)效率的幾點(diǎn)思考
    甘肅教育(2020年14期)2020-09-11 07:57:42
    重尾非線性自回歸模型自加權(quán)M-估計(jì)的漸近分布
    3D打印中的模型分割與打包
    跟蹤導(dǎo)練(一)2
    久久久精品94久久精品| 中文字幕最新亚洲高清| 国产成人a∨麻豆精品| 国产精品熟女久久久久浪| 久久人人97超碰香蕉20202| 国产成人午夜福利电影在线观看| 成年美女黄网站色视频大全免费| 精品国产乱码久久久久久小说| 欧美日韩亚洲国产一区二区在线观看 | 久久韩国三级中文字幕| 亚洲欧美日韩另类电影网站| 国产一区有黄有色的免费视频| 欧美激情极品国产一区二区三区| 青草久久国产| 亚洲精品久久午夜乱码| 香蕉精品网在线| 国产不卡av网站在线观看| 人成视频在线观看免费观看| 亚洲综合色网址| 亚洲成人手机| 男女午夜视频在线观看| 伦精品一区二区三区| 成人午夜精彩视频在线观看| 久久久久国产网址| 在线天堂最新版资源| 极品人妻少妇av视频| 久久 成人 亚洲| 蜜桃国产av成人99| 欧美+日韩+精品| 这个男人来自地球电影免费观看 | 9191精品国产免费久久| 午夜福利视频在线观看免费| 寂寞人妻少妇视频99o| 成人黄色视频免费在线看| 久久久久精品久久久久真实原创| 久久精品人人爽人人爽视色| 狠狠精品人妻久久久久久综合| 少妇 在线观看| 天堂俺去俺来也www色官网| 乱人伦中国视频| 亚洲精品久久成人aⅴ小说| 久久鲁丝午夜福利片| 国产高清不卡午夜福利| 老鸭窝网址在线观看| 一区在线观看完整版| 成人国产av品久久久| 高清欧美精品videossex| 欧美日韩精品网址| 久热久热在线精品观看| 亚洲av在线观看美女高潮| 国产免费又黄又爽又色| 午夜福利一区二区在线看| 老汉色av国产亚洲站长工具| 大陆偷拍与自拍| 亚洲一区二区三区欧美精品| 91久久精品国产一区二区三区| 免费日韩欧美在线观看| 人妻 亚洲 视频| 老鸭窝网址在线观看| 国产精品三级大全| 啦啦啦啦在线视频资源| 久久久久精品人妻al黑| 亚洲伊人久久精品综合| 中文乱码字字幕精品一区二区三区| 中文字幕制服av| 亚洲国产看品久久| 赤兔流量卡办理| 国产人伦9x9x在线观看 | 国产精品久久久久成人av| 少妇熟女欧美另类| 十八禁高潮呻吟视频| 老汉色∧v一级毛片| 老熟女久久久| 国产精品蜜桃在线观看| 热99久久久久精品小说推荐| 性高湖久久久久久久久免费观看| 国产精品二区激情视频| 久久久久国产网址| 亚洲精品久久成人aⅴ小说| 久久精品久久久久久久性| 欧美成人午夜精品| 久久人人97超碰香蕉20202| av视频免费观看在线观看| 亚洲精品美女久久av网站| 免费观看av网站的网址| 边亲边吃奶的免费视频| 国产精品偷伦视频观看了| 少妇人妻久久综合中文| 国产人伦9x9x在线观看 | 亚洲,欧美精品.| 亚洲av日韩在线播放| 一级毛片黄色毛片免费观看视频| 下体分泌物呈黄色| 国产一级毛片在线| 十八禁高潮呻吟视频| 精品亚洲成国产av| 日韩视频在线欧美| 国产日韩欧美在线精品| 人妻 亚洲 视频| 国产男女超爽视频在线观看| 男的添女的下面高潮视频| 免费人妻精品一区二区三区视频| 咕卡用的链子| 亚洲精品成人av观看孕妇| 伊人久久国产一区二区| 国产成人午夜福利电影在线观看| 尾随美女入室| 国产免费现黄频在线看| 电影成人av| 精品一区二区三区四区五区乱码 | 另类亚洲欧美激情| 国产极品粉嫩免费观看在线| 亚洲国产最新在线播放| 一区福利在线观看| 哪个播放器可以免费观看大片| 最近中文字幕高清免费大全6| 国产无遮挡羞羞视频在线观看| 亚洲,欧美精品.| 99香蕉大伊视频| 男男h啪啪无遮挡| 国产亚洲精品第一综合不卡| 肉色欧美久久久久久久蜜桃| 亚洲av电影在线观看一区二区三区| 日韩一区二区三区影片| 看十八女毛片水多多多| 亚洲综合色网址| 久久人人爽人人片av| 亚洲精华国产精华液的使用体验| 69精品国产乱码久久久| 欧美 亚洲 国产 日韩一| 精品一品国产午夜福利视频| 制服诱惑二区| 老司机亚洲免费影院| 80岁老熟妇乱子伦牲交| 午夜av观看不卡| 热re99久久国产66热| tube8黄色片| 天天躁夜夜躁狠狠久久av| 激情五月婷婷亚洲| 狠狠精品人妻久久久久久综合| 捣出白浆h1v1| 国产精品人妻久久久影院| 亚洲av免费高清在线观看| 久久这里有精品视频免费| 亚洲国产精品一区二区三区在线| 男女国产视频网站| av线在线观看网站| 亚洲激情五月婷婷啪啪| 久久99热这里只频精品6学生| 晚上一个人看的免费电影| 成人免费观看视频高清| 日韩中文字幕视频在线看片| 人妻少妇偷人精品九色| 国产av一区二区精品久久| 黄片无遮挡物在线观看| 国产免费一区二区三区四区乱码| 久久婷婷青草| 欧美bdsm另类| a级毛片黄视频| 巨乳人妻的诱惑在线观看| 久久久久国产网址| 久久久精品区二区三区| 九九爱精品视频在线观看| 在线观看人妻少妇| 日韩,欧美,国产一区二区三区| 国产黄色免费在线视频| 久久影院123| 亚洲精品成人av观看孕妇| 久久久欧美国产精品| 大话2 男鬼变身卡| 女人精品久久久久毛片| 丝袜美足系列| 久久久久久免费高清国产稀缺| 亚洲国产毛片av蜜桃av| 久久久久国产网址| 国产成人精品久久二区二区91 | 啦啦啦在线免费观看视频4| 国产人伦9x9x在线观看 | av不卡在线播放| 久久久久久久国产电影| 在线精品无人区一区二区三| 日韩人妻精品一区2区三区| 亚洲成国产人片在线观看| 日韩成人av中文字幕在线观看| 美女国产高潮福利片在线看| 亚洲少妇的诱惑av| 成人亚洲精品一区在线观看| 五月伊人婷婷丁香| 多毛熟女@视频| 亚洲婷婷狠狠爱综合网| 丰满少妇做爰视频| 午夜日韩欧美国产| 中文乱码字字幕精品一区二区三区| av女优亚洲男人天堂| 国产日韩欧美亚洲二区| 欧美日本中文国产一区发布| 欧美日韩视频精品一区| av天堂久久9| 亚洲国产欧美网| 精品国产超薄肉色丝袜足j| 男女午夜视频在线观看| 日韩精品免费视频一区二区三区| 亚洲欧美色中文字幕在线| 国产精品一区二区在线观看99| 在线看a的网站| 精品久久久久久电影网| 日本爱情动作片www.在线观看| 日本爱情动作片www.在线观看| 水蜜桃什么品种好| 在线观看人妻少妇| 中文字幕另类日韩欧美亚洲嫩草| 欧美亚洲 丝袜 人妻 在线| 久久精品久久久久久噜噜老黄| 久久久久国产一级毛片高清牌| 男女边吃奶边做爰视频| 国产av国产精品国产| 女性被躁到高潮视频| 在线亚洲精品国产二区图片欧美| 久久精品国产a三级三级三级| 亚洲欧美色中文字幕在线| 久久久精品免费免费高清| 国产精品av久久久久免费| www.精华液| 黑人猛操日本美女一级片| 色播在线永久视频| 免费黄色在线免费观看| 亚洲欧美色中文字幕在线| 午夜福利在线免费观看网站| 中文字幕人妻熟女乱码| 久久午夜综合久久蜜桃| 欧美国产精品va在线观看不卡| 免费日韩欧美在线观看| 91成人精品电影| 国产白丝娇喘喷水9色精品| 国产精品无大码| 中文字幕最新亚洲高清| 天堂俺去俺来也www色官网| 女性生殖器流出的白浆| 人人妻人人添人人爽欧美一区卜| 国产成人精品一,二区| 免费观看a级毛片全部| 一本久久精品| 久久久久精品人妻al黑| 国产亚洲午夜精品一区二区久久| 国产在线一区二区三区精| 亚洲一区二区三区欧美精品| 国产片特级美女逼逼视频| 亚洲精品美女久久久久99蜜臀 | 国产男女超爽视频在线观看| 免费观看性生交大片5| 黄网站色视频无遮挡免费观看| 亚洲欧美成人精品一区二区| 日日撸夜夜添| 丝袜美足系列| 亚洲欧美成人精品一区二区| 精品一区二区免费观看| 水蜜桃什么品种好| 国产成人精品久久二区二区91 | 夜夜骑夜夜射夜夜干| 人人妻人人爽人人添夜夜欢视频| 国产成人精品婷婷| 日日啪夜夜爽| 又大又黄又爽视频免费| 婷婷成人精品国产| 只有这里有精品99| av电影中文网址| 十八禁网站网址无遮挡| 在线观看国产h片| 免费少妇av软件| 日韩精品有码人妻一区| 日韩中文字幕欧美一区二区 | 一区二区日韩欧美中文字幕| 国产97色在线日韩免费| 午夜精品国产一区二区电影| 国产黄色免费在线视频| 美女国产高潮福利片在线看| 色网站视频免费| xxxhd国产人妻xxx| 久久久精品免费免费高清| 成人国产av品久久久| 精品亚洲成国产av| 日韩一卡2卡3卡4卡2021年| 亚洲av综合色区一区| 天天躁夜夜躁狠狠久久av| 色播在线永久视频| 免费大片黄手机在线观看| 国产高清国产精品国产三级| 国产精品麻豆人妻色哟哟久久| 成人黄色视频免费在线看| 久久精品夜色国产| 午夜福利在线观看免费完整高清在| 最近中文字幕高清免费大全6| 国产淫语在线视频| 成人国产av品久久久| 精品一区二区三卡| 亚洲国产欧美在线一区| 九九爱精品视频在线观看| 午夜日本视频在线| 亚洲欧美精品综合一区二区三区 | 久久久久久久久免费视频了| 欧美老熟妇乱子伦牲交| 精品一区在线观看国产| 午夜免费观看性视频| 美女国产视频在线观看| 国产欧美日韩一区二区三区在线| 大香蕉久久成人网| 欧美成人午夜精品| 日产精品乱码卡一卡2卡三| 99久久精品国产国产毛片| 99re6热这里在线精品视频| 国产亚洲午夜精品一区二区久久| 极品少妇高潮喷水抽搐| 精品一区二区三区四区五区乱码 | 国产一区二区 视频在线| 精品一品国产午夜福利视频| 国产男女内射视频| 亚洲精品久久成人aⅴ小说| av免费在线看不卡| 韩国精品一区二区三区| 久久精品国产自在天天线| 妹子高潮喷水视频| 只有这里有精品99| 少妇 在线观看| 久久精品国产亚洲av高清一级| 国精品久久久久久国模美| 亚洲成av片中文字幕在线观看 | 精品人妻一区二区三区麻豆| 少妇被粗大猛烈的视频| 一区二区三区激情视频| 免费高清在线观看视频在线观看| 亚洲欧美中文字幕日韩二区| 午夜福利在线免费观看网站| 亚洲国产欧美网| 女人久久www免费人成看片| av一本久久久久| 午夜福利乱码中文字幕| 欧美中文综合在线视频| 精品人妻在线不人妻| 777米奇影视久久| 女人高潮潮喷娇喘18禁视频| 日本wwww免费看| 精品99又大又爽又粗少妇毛片| 精品人妻偷拍中文字幕| 午夜日韩欧美国产| 女人被躁到高潮嗷嗷叫费观| 亚洲欧美日韩另类电影网站| 捣出白浆h1v1| 一本色道久久久久久精品综合| 欧美精品国产亚洲| 久久久久国产精品人妻一区二区| 深夜精品福利| 亚洲av男天堂| 亚洲一区中文字幕在线| 在线精品无人区一区二区三| 一本色道久久久久久精品综合| 男人爽女人下面视频在线观看| 尾随美女入室| 日本vs欧美在线观看视频| 在线免费观看不下载黄p国产| 高清视频免费观看一区二区| 十分钟在线观看高清视频www| 秋霞在线观看毛片| 成人亚洲精品一区在线观看| 国产人伦9x9x在线观看 | 日本-黄色视频高清免费观看| 99久久综合免费| 国产成人精品在线电影| 性色avwww在线观看| 亚洲av电影在线进入| 寂寞人妻少妇视频99o| 国产综合精华液| 久久av网站| 精品第一国产精品| 亚洲成av片中文字幕在线观看 | av在线app专区| 国产成人免费观看mmmm| tube8黄色片| 国产97色在线日韩免费| 精品一区二区三卡| av一本久久久久| kizo精华| 精品第一国产精品| 久久久久国产网址| 国产成人免费无遮挡视频| 大香蕉久久网| 黄色毛片三级朝国网站| 国产成人精品在线电影| 我要看黄色一级片免费的| 精品国产乱码久久久久久男人| 精品国产国语对白av| 色网站视频免费| 亚洲精品国产av成人精品| 色网站视频免费| 亚洲欧美清纯卡通| 你懂的网址亚洲精品在线观看| 丰满少妇做爰视频| 国产福利在线免费观看视频| 97精品久久久久久久久久精品| 国产精品 国内视频| 亚洲欧美日韩另类电影网站| av在线观看视频网站免费| 久久人人爽人人片av| 久久免费观看电影| 激情视频va一区二区三区| 香蕉国产在线看| 18禁动态无遮挡网站| 一区二区av电影网| av在线老鸭窝| 在线观看三级黄色| 久久久欧美国产精品| 丰满乱子伦码专区| 1024香蕉在线观看| 久久青草综合色| 精品酒店卫生间| videossex国产| 久久久久久久久久人人人人人人| 亚洲综合精品二区| 电影成人av| 国产精品一二三区在线看| 国产黄频视频在线观看| 晚上一个人看的免费电影| 欧美日韩精品成人综合77777| 欧美亚洲日本最大视频资源| 成人国产麻豆网| 999久久久国产精品视频| 日本vs欧美在线观看视频| 亚洲欧美色中文字幕在线| 有码 亚洲区| 欧美激情 高清一区二区三区| 国产一区二区 视频在线| 欧美日本中文国产一区发布| 咕卡用的链子| 伊人亚洲综合成人网| 亚洲三区欧美一区| 免费播放大片免费观看视频在线观看| 伊人久久国产一区二区| 大片电影免费在线观看免费| 黄片小视频在线播放| 亚洲色图综合在线观看| 熟女电影av网| 狂野欧美激情性bbbbbb| 国产精品一国产av| 极品少妇高潮喷水抽搐| 男女边吃奶边做爰视频| 国产伦理片在线播放av一区| 免费观看无遮挡的男女| 国产有黄有色有爽视频| av天堂久久9| 一级片免费观看大全| 下体分泌物呈黄色| 久久精品亚洲av国产电影网| 999久久久国产精品视频| 桃花免费在线播放| 日韩视频在线欧美| 国产精品偷伦视频观看了| 丝袜在线中文字幕| 性少妇av在线| 成年女人在线观看亚洲视频| 国产成人av激情在线播放| 成人免费观看视频高清| 女性被躁到高潮视频| 免费高清在线观看视频在线观看| 另类精品久久| 日日撸夜夜添| 久久狼人影院| 国产精品一区二区在线不卡| 久久久久久久久久久免费av| 国产在视频线精品| 国产黄色免费在线视频| 伊人久久大香线蕉亚洲五| 亚洲精品乱久久久久久| 亚洲精品一二三| 国产成人免费无遮挡视频| 免费观看av网站的网址| 日韩中文字幕欧美一区二区 | 国产伦理片在线播放av一区| 男男h啪啪无遮挡| 亚洲av日韩在线播放| 青春草国产在线视频| 最新中文字幕久久久久| 亚洲国产看品久久| 青青草视频在线视频观看| 久久精品国产亚洲av涩爱| 王馨瑶露胸无遮挡在线观看| 国产精品 国内视频| 亚洲国产欧美日韩在线播放| 日韩大片免费观看网站| 亚洲欧美成人精品一区二区| 国产男女内射视频| 欧美97在线视频| 99香蕉大伊视频| 国产一区亚洲一区在线观看| 亚洲欧美一区二区三区国产| 香蕉丝袜av| 黄色毛片三级朝国网站| 亚洲人成77777在线视频| 久久久久精品性色| 日本午夜av视频| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲av电影在线观看一区二区三区| 亚洲久久久国产精品| 精品国产国语对白av| 天天躁夜夜躁狠狠久久av| 激情视频va一区二区三区| 日韩一卡2卡3卡4卡2021年| 久热这里只有精品99| 久久精品久久精品一区二区三区| 久久久国产精品麻豆| 欧美人与性动交α欧美软件| 亚洲一区二区三区欧美精品| 乱人伦中国视频| 国产激情久久老熟女| 伦精品一区二区三区| av国产精品久久久久影院| 天堂中文最新版在线下载| 在线观看www视频免费| 久久久久久伊人网av| 精品一品国产午夜福利视频| 黄色毛片三级朝国网站| 亚洲精品aⅴ在线观看| 欧美 日韩 精品 国产| 亚洲av福利一区| 十八禁网站网址无遮挡| 国产精品嫩草影院av在线观看| 一级,二级,三级黄色视频| 精品99又大又爽又粗少妇毛片| 少妇的丰满在线观看| 国产伦理片在线播放av一区| 亚洲精品国产av成人精品| 日韩av在线免费看完整版不卡| 丝袜在线中文字幕| 午夜激情久久久久久久| 黄频高清免费视频| 超碰97精品在线观看| 国产成人免费观看mmmm| 国产伦理片在线播放av一区| 一本大道久久a久久精品| 成人手机av| 在线观看人妻少妇| 欧美日韩视频高清一区二区三区二| 国产免费现黄频在线看| 侵犯人妻中文字幕一二三四区| 欧美日韩精品成人综合77777| 天天躁狠狠躁夜夜躁狠狠躁| 观看美女的网站| 成年美女黄网站色视频大全免费| 青春草国产在线视频| 成人国产av品久久久| 精品一品国产午夜福利视频| 久久 成人 亚洲| 国产精品 国内视频| 韩国精品一区二区三区| 男女国产视频网站| 国产精品av久久久久免费| 一二三四在线观看免费中文在| 亚洲天堂av无毛| 久久国产精品男人的天堂亚洲| 美女xxoo啪啪120秒动态图| 亚洲人成77777在线视频| 国产乱来视频区| 最近最新中文字幕免费大全7| 少妇被粗大的猛进出69影院| 亚洲欧美一区二区三区久久| 亚洲精品久久久久久婷婷小说| 黑人巨大精品欧美一区二区蜜桃| √禁漫天堂资源中文www| 久久女婷五月综合色啪小说| 侵犯人妻中文字幕一二三四区| 精品国产一区二区久久| 亚洲精品一二三| 国产精品久久久久成人av| 亚洲熟女精品中文字幕| 久久狼人影院| 99久久综合免费| 久久99精品国语久久久| videos熟女内射| 日本av手机在线免费观看| 午夜福利在线观看免费完整高清在| 国产精品一二三区在线看| 亚洲av男天堂| 黄片小视频在线播放| 超碰成人久久| 最黄视频免费看| 亚洲人成网站在线观看播放| 丝袜美足系列| 街头女战士在线观看网站| 精品国产一区二区三区四区第35| 人人妻人人添人人爽欧美一区卜| 欧美xxⅹ黑人| 自线自在国产av| 两个人看的免费小视频| 高清黄色对白视频在线免费看| 亚洲av综合色区一区| 国产精品免费大片| 色网站视频免费| 老汉色∧v一级毛片| 午夜久久久在线观看| 久久久久久久久久人人人人人人| 一边摸一边做爽爽视频免费| 精品久久蜜臀av无| 女人久久www免费人成看片| 欧美成人午夜精品| 在线观看三级黄色| 90打野战视频偷拍视频| 久久精品久久久久久久性| 精品久久蜜臀av无| 久久精品人人爽人人爽视色| 韩国高清视频一区二区三区| 中文精品一卡2卡3卡4更新| 亚洲国产日韩一区二区| 在线免费观看不下载黄p国产| 99热网站在线观看| 亚洲男人天堂网一区| 人妻少妇偷人精品九色| 超色免费av|