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

    葉片結(jié)構(gòu)對(duì)軸流導(dǎo)葉式旋風(fēng)分離器性能的影響分析

    2021-01-08 08:24:18楊維旺王燕云于洪傳
    流體機(jī)械 2020年12期
    關(guān)鍵詞:切向速度圈數(shù)軸流

    楊維旺,王燕云,于洪傳

    (1.中車戚墅堰機(jī)車車輛工藝研究所有限公司,江蘇常州 213000;2.青島威奧軌道股份有限公司,山東青島 266108;3.中車青島四方機(jī)車車輛股份有限公司,山東青島 266108)

    0 引言

    旋風(fēng)分離器是一種常見的氣液分離設(shè)備,具有結(jié)構(gòu)簡(jiǎn)單、無運(yùn)動(dòng)部件、分離效率高、維修方便等優(yōu)點(diǎn),被廣泛應(yīng)用在石油、化工、環(huán)保等領(lǐng)域[1]。隨著全球能源短缺和對(duì)環(huán)境保護(hù)的重視,提高分離效率并降低壓降已經(jīng)成為石化行業(yè)新型分離器研究的熱點(diǎn)。

    目前國(guó)內(nèi)外對(duì)旋風(fēng)分離器性能的研究主要集中在切向入口的旋風(fēng)分離器上,而對(duì)軸向入口的旋風(fēng)分離器的研究則相對(duì)較少。軸向入口的旋風(fēng)分離器的結(jié)構(gòu)緊湊,徑向尺寸小,分離效率高,壓力損失小,近年來受到越來越多研究人員的關(guān)注[2-3]。王振波等[4]通過試驗(yàn)研究了一種新型分離器的葉道出口角度、錐體角度等因素對(duì)分離性能的影響情況。Win field等[5]對(duì)單切向入口和三切向入口的旋風(fēng)分離器進(jìn)行了數(shù)值模擬,結(jié)果表明,與單切向入口比較,三切向入口的旋風(fēng)分離器的效率不變,壓降更低。Lee等[6]通過試驗(yàn)測(cè)定了單管旋風(fēng)分離器和雙管旋風(fēng)分離器對(duì)于分離固體顆粒的分離效率和壓降的差異,試驗(yàn)結(jié)果表明,雙管旋風(fēng)分離器的效率較高,但壓降較低。金向紅等[7]通過試驗(yàn)對(duì)比研究了管錐式分離器與管柱式分離器的分離效率與壓降情況,結(jié)果表明管柱式分離器分離效率較高,管錐式的壓降損失較小。Yang等[8]通過試驗(yàn)和數(shù)值模擬的方法,研究分離效率最大值時(shí)的分離器高度,結(jié)果表明最大效率旋風(fēng)分離器高度隨著入口尺寸減小而降低。

    目前,在石油催化裂化工藝過程中,大量被應(yīng)用于油氣過濾和凈化的旋風(fēng)分離器結(jié)構(gòu)陳舊,進(jìn)出口的壓降很大,分離效率也不理想。本文利用FLUENT軟件對(duì)軸流導(dǎo)葉式旋風(fēng)分離器的分離特性進(jìn)行數(shù)值模擬,同時(shí)進(jìn)行試驗(yàn)驗(yàn)證。文中分析了軸流導(dǎo)葉式旋風(fēng)分離器內(nèi)部流場(chǎng)的分布規(guī)律,模擬計(jì)算了不同葉片結(jié)構(gòu)的分離器內(nèi)部壓力場(chǎng)、速度場(chǎng)、分離效率和壓降值,進(jìn)而分析了旋風(fēng)分離器葉片結(jié)構(gòu)對(duì)其分離效率和壓降的影響規(guī)律,給出了具體的葉片高度和葉片圈數(shù)的優(yōu)化參考建議。

    1 幾何模型

    本文的研究是基于一項(xiàng)發(fā)明專利和實(shí)用新型專利[9-10],對(duì)專利中的軸流導(dǎo)葉式旋風(fēng)分離器進(jìn)行結(jié)構(gòu)上的優(yōu)化研究。本文利用Solidworks軟件建立旋風(fēng)分離器的三維模型如圖1所示。

    圖1 軸流導(dǎo)葉式氣液旋風(fēng)分離器結(jié)構(gòu)

    軸流導(dǎo)葉式旋風(fēng)分離器由排氣管、導(dǎo)流葉片、外筒體和儲(chǔ)液槽構(gòu)成,其中外筒體可劃分為3部分:柱筒部分、錐筒部分和排液管。排液管底部外壁面與儲(chǔ)液槽內(nèi)壁面用連接板焊接相連,可進(jìn)行液封,保證氣流全部從排氣管排出。如圖1(a)所示,分離器的總高度為1 230 mm,其中外筒體高度H1=1 120 mm,柱筒高度H2=120 mm,錐筒高度H3=150 mm。排液管的高度為850 mm,柱筒外徑D1=76 mm,排液管外徑D2=34 mm,壁厚均為3 mm。排氣管高度H4=190 mm,排氣管插入柱筒段深度H5=100 mm,排氣管外徑D3=51 mm,壁厚為3 mm。儲(chǔ)液槽高度H7=80 mm,排液管插入儲(chǔ)液槽深度H8=60 mm,連接儲(chǔ)液槽與排液管的連接板高度H9=35 mm,儲(chǔ)液槽外徑D4=90 mm,連接板和儲(chǔ)液槽的壁厚均為3 mm。圖1(b)(c)分別示出分離器內(nèi)部含導(dǎo)流葉片的分離空間橫截面和儲(chǔ)液槽的橫截面。圖中所示的旋風(fēng)分離器的導(dǎo)流葉片數(shù)量為4,葉片高度H6=50 mm,葉片厚度為2 mm,葉片的結(jié)構(gòu)如圖1(d)(e)所示,葉片螺旋角為 60°,葉片旋轉(zhuǎn)角為 90°。

    2 數(shù)值模擬

    2.1 計(jì)算模型

    目前對(duì)分離器內(nèi)部流場(chǎng)的模擬計(jì)算,較多使用標(biāo)準(zhǔn)k-ε模型、RNG k-ε模型和雷諾應(yīng)力模型(RSM)。對(duì)這3種模型,前人進(jìn)行過大量的研究,其中RSM的計(jì)算結(jié)果與實(shí)際測(cè)量結(jié)果最吻合[11-13]。因?yàn)?,RSM摒棄了基于各向同性渦黏性的Boussinesq假設(shè),包含更多物理過程的影響,考慮了湍流各向異性,特別是旋轉(zhuǎn)效應(yīng)、浮力效應(yīng)、曲率效應(yīng)等[14]。

    RSM模型是直接對(duì)Reynolds應(yīng)力建立微分方程式,并對(duì)連續(xù)性方程、動(dòng)量方程、雷諾應(yīng)力方程及湍動(dòng)能方程建立方程組進(jìn)行求解:

    Reynolds應(yīng)力方程:

    Cij——對(duì)流項(xiàng);

    Dij——擴(kuò)散項(xiàng);

    DTij——湍流擴(kuò)散項(xiàng);

    DLij——分子黏性擴(kuò)散項(xiàng);

    Pij——雷諾剪應(yīng)力產(chǎn)生項(xiàng);

    φij——壓力應(yīng)變項(xiàng);

    εij——黏性耗散項(xiàng)。

    在上式各項(xiàng)中,Cij,DLij,Pij均只包含二階關(guān)聯(lián)項(xiàng),不必進(jìn)行模化處理。DTij,φij,εij包含未知的關(guān)聯(lián)項(xiàng),需模擬為表示的式子,才能完成對(duì)雷諾應(yīng)力輸運(yùn)方程的封閉。式中各項(xiàng)的?;椒?、計(jì)算參數(shù)選擇及湍動(dòng)方程可參考文獻(xiàn)[15-16]。

    軸流導(dǎo)葉式旋風(fēng)分離器內(nèi)的氣體流動(dòng)是三維強(qiáng)旋轉(zhuǎn)湍流流動(dòng),運(yùn)動(dòng)形式比較復(fù)雜,在計(jì)算中采用SIMPLEC算法及QUICK差分格式[17-18]。對(duì)動(dòng)量方程和湍流強(qiáng)度方程,是采用基于壓力速度耦合的方法進(jìn)行求解的[19]。

    2.2 邊界條件

    2.2.1 入口邊界

    在常溫常壓狀態(tài)下,入口空氣的密度ρ為1.237 kg/m3,空氣動(dòng)力黏度μ為 1.807×10-5Pa·s,空氣摩爾質(zhì)量為28.96 g/mol,氣流進(jìn)入方向?yàn)榇怪庇诜蛛x空間的環(huán)形斷面。模擬計(jì)算中,進(jìn)氣速度依次設(shè)置為 5,7.5,10,12.5 m/s。

    2.2.2 出口邊界

    在模擬計(jì)算中,氣液兩相流從分離器入口進(jìn)入,大量液滴在分離空間被捕集,沿分離器內(nèi)壁流入排液管中。排液管底部有液封結(jié)構(gòu),氣流無法通過,氣體全部由排氣管排出,部分液滴隨氣流由排氣管逃逸,未被分離。排氣管的出口邊界設(shè)置為out flow,DPM選項(xiàng)設(shè)為escape,即為氣液兩相流計(jì)算終止處。排液管口設(shè)為out flow,DPM選項(xiàng)設(shè)為trap。排液管口處有液體密封,沒有氣流流出,所以取排液管口處的流量權(quán)重為0,排氣管口處的流量權(quán)重為1。

    2.2.3 壁面邊界條件

    模擬計(jì)算中,壁面效應(yīng)是湍流和旋渦的重要來源,因此在近壁區(qū)的處理對(duì)數(shù)值計(jì)算結(jié)果的準(zhǔn)確性有顯著的影響[20-23]。本文模擬計(jì)算前對(duì)筒壁面和導(dǎo)流葉片各處進(jìn)行合理劃分,設(shè)置了不同的碰撞恢復(fù)情況,見表1。

    表1 不同壁面位置的碰撞恢復(fù)情況

    2.3 網(wǎng)格劃分

    本文運(yùn)用Gambit軟件生成網(wǎng)格模型,如圖2所示。為了實(shí)現(xiàn)旋風(fēng)分離器網(wǎng)格的整體優(yōu)化,將旋風(fēng)分離器內(nèi)部的流場(chǎng)分為排氣管內(nèi)、分離空間、錐筒內(nèi)和儲(chǔ)液槽空間4個(gè)部分。在分離空間和錐筒內(nèi)區(qū)域,氣流運(yùn)動(dòng)比較劇烈,壁面作用使流體存在較大的速度梯度,為了能夠準(zhǔn)確模擬壁面附近流場(chǎng),采用邊界層網(wǎng)格和COOPER網(wǎng)格生成技術(shù)。其它區(qū)域,為保證網(wǎng)格質(zhì)量同時(shí)提高計(jì)算的經(jīng)濟(jì)性,采用四面體網(wǎng)格,圖中的網(wǎng)格模型包含約13萬個(gè)網(wǎng)格單元。

    圖2 旋風(fēng)分離器的網(wǎng)格模型

    3 模擬結(jié)果分析

    3.1 壓力分布

    為了研究旋風(fēng)分離器內(nèi)部不同位置的流場(chǎng)分布情況,在主筒區(qū)域即排氣管底部到錐筒頂部之間流場(chǎng)區(qū)域選取6個(gè)不同高度的橫截面,對(duì)計(jì)算結(jié)果進(jìn)行對(duì)比分析。6個(gè)橫截面與分離器最底端的距離分別為h1=0.890 m,h2=0.895 m,h3=0.900 m,h4=0.905 m,h5=0.910 m,h6=0.915 m,如圖3所示。

    圖3 不同高度的橫截面

    圖4示出旋風(fēng)分離器內(nèi)部流場(chǎng)在Y=0截面上的靜壓力分布云圖。

    圖4 旋風(fēng)分離器Y=0截面靜壓值分布云圖

    由圖可以看出,對(duì)于軸流導(dǎo)葉式旋風(fēng)分離器,內(nèi)部流場(chǎng)分為內(nèi)旋流和外旋流兩部分,外旋流的靜壓值較高,內(nèi)旋流的靜壓值較低。沿徑向上靜壓值由內(nèi)向外逐漸增大,靠近分離器壁面上的靜壓值最大。在中心軸線上的壓力很小,在排氣管入口中心處出現(xiàn)了負(fù)壓區(qū),達(dá)到-300 Pa左右。沿軸向上靜壓值變化不大,在柱筒和錐筒接合斷面處分離器壁面存在較大的靜壓值,達(dá)到700 Pa。

    對(duì)不同截面的靜壓分布情況進(jìn)行分析,由圖5可以看出,軸流導(dǎo)葉式旋風(fēng)分離器內(nèi)部氣流靜壓力場(chǎng)在一個(gè)橫截面上沿徑向上的分布曲線呈V型。這表明,在同一橫截面上,中心處的靜壓值最小,沿徑向向外逐漸增大。在不同高度橫截面上,靠近壁面的靜壓值變化很小,在中心位置的靜壓值變化比較明顯。高度越高,壓力越低。

    圖5 旋風(fēng)分離器不同橫截面的靜壓分布

    3.2 速度分布

    圖6示出旋風(fēng)分離器內(nèi)部氣流在Y=0截面上的軸向速度、切向速度分布云圖。

    圖6 旋風(fēng)分離器Y=0截面軸向速度、切向速度分布云圖

    由圖6(a)可知,軸流導(dǎo)葉式旋風(fēng)分離器的內(nèi)部氣流的軸向速度分布沿中心軸線左右對(duì)稱,在排氣管入口中心處的軸向速度達(dá)到最大值,在排氣管入口外側(cè)的軸向速度值最小,且二者方向相反。由圖6(b)中可以看出,軸流導(dǎo)葉式旋風(fēng)分離器內(nèi)部的氣流切向速度分布沿中心軸線上左右對(duì)稱。在排氣管入口與錐筒上半段存在一個(gè)切向速度較大的區(qū)段,排氣管入口處內(nèi)側(cè)切向速度達(dá)到最大值,而外側(cè)切向速度最小。綜上表明,在排氣管入口區(qū)域,氣流發(fā)生轉(zhuǎn)向,由向下運(yùn)動(dòng)轉(zhuǎn)為向上流動(dòng),分離器內(nèi)部氣流形成了二次渦流。這是由于進(jìn)入旋風(fēng)分離器的氣流,在葉片導(dǎo)流作用下,由初始的向下軸向運(yùn)動(dòng)轉(zhuǎn)化為旋轉(zhuǎn)運(yùn)動(dòng),形成外層的渦旋運(yùn)動(dòng),在到達(dá)分離器底部后,又轉(zhuǎn)而向上運(yùn)動(dòng),形成內(nèi)部向上的渦流運(yùn)動(dòng)。因此,軸流導(dǎo)葉式旋風(fēng)分離器內(nèi)的氣流分為2個(gè)部分,即外部的準(zhǔn)自由渦和內(nèi)部的準(zhǔn)強(qiáng)制渦。內(nèi)外漩渦在排氣管入口處交匯,氣流運(yùn)動(dòng)較為劇烈且發(fā)生轉(zhuǎn)向,同時(shí)存在返混現(xiàn)象和旋進(jìn)渦核現(xiàn)象。

    圖7示出旋風(fēng)分離器內(nèi)部氣流在不同截面上的軸向速度、切向速度分布曲線。從圖7(a)看出,氣流的軸向速度沿徑向上呈駝峰狀分布,且內(nèi)外旋流軸向流動(dòng)的方向相反。從圖7(b)可知,切向速度沿徑向上呈雙峰狀,由內(nèi)向外,切向速度先升高后降低。切向速度始終同一方向,表明在同一截面上的內(nèi)旋流和外旋流的旋轉(zhuǎn)方向是同向的。此外,圖中切向速度值最大處為排氣管入口附近,即是內(nèi)外旋流的交匯處。而在分離器壁面處的切向速度最小,接近于0,這說明在壁面處的氣流幾乎沒有旋轉(zhuǎn)流動(dòng)。

    圖7 旋風(fēng)分離器不同截面的軸向速度、切向速度分布

    4 試驗(yàn)驗(yàn)證

    為了驗(yàn)證數(shù)值模擬的可信性,作者搭建了試驗(yàn)裝置,進(jìn)行了試驗(yàn)驗(yàn)證。試驗(yàn)中測(cè)定了軸流導(dǎo)葉式旋風(fēng)分離器在不同霧化液滴粒徑下的分離效率和不同進(jìn)氣速度下的壓降值,然后與數(shù)值計(jì)算結(jié)果對(duì)比。搭建了試驗(yàn)裝置。

    在試驗(yàn)中,為了產(chǎn)生氣液混合兩相流,對(duì)液體進(jìn)行霧化,液體選用癸二酸二辛酯。采用孔徑為0.46 mm的單流體霧化噴嘴進(jìn)行液體的霧化,霧化噴嘴孔徑一定時(shí),影響霧化液滴粒徑的因素是液體泵的操作壓力[24]。通過調(diào)節(jié)泵的操作壓力,進(jìn)而控制霧化產(chǎn)生液滴的平均粒徑。本試驗(yàn)中,霧化液滴的平均粒徑與液體泵的操作壓力關(guān)系如圖8所示。

    圖8 液體泵操作壓力與霧化液滴平均粒徑關(guān)系

    本試驗(yàn)中旋風(fēng)分離器葉片高度為50 mm,結(jié)構(gòu)尺寸與模擬計(jì)算中完全相同。在壓降測(cè)定試驗(yàn)中,調(diào)節(jié)風(fēng)機(jī)來控制進(jìn)氣速度,進(jìn)氣速度依次為2.5,5.0,7.5,10.0,12.5 m/s。在效率測(cè)定試驗(yàn)中,調(diào)節(jié)泵的壓力來控制霧化液滴粒徑,依次為5,6,7,8,9,10 μm。

    圖9示出試驗(yàn)測(cè)定的分離效率結(jié)果與模擬計(jì)算結(jié)果的對(duì)比。

    圖9 分離效率對(duì)比

    由圖可以看出,二者相差很小,平均誤差為2.32%。試驗(yàn)測(cè)定的壓降結(jié)果與模擬計(jì)算結(jié)果對(duì)比如圖10所示,二者也比較吻合,平均誤差為2.17%。通過以上試驗(yàn)與模擬的對(duì)比結(jié)果可以說明,本文采用的模擬計(jì)算方法合理,結(jié)果準(zhǔn)確。

    圖10 壓降對(duì)比

    5 葉片結(jié)構(gòu)影響分析

    在驗(yàn)證模擬計(jì)算的準(zhǔn)確性之后,對(duì)不同葉片結(jié)構(gòu)的旋風(fēng)分離器內(nèi)部流場(chǎng)進(jìn)行模擬計(jì)算,并對(duì)計(jì)算結(jié)果進(jìn)行對(duì)比分析。本文構(gòu)建了2組不同的葉片高度和葉片圈數(shù)的旋風(fēng)分離器,分別進(jìn)行數(shù)值模擬計(jì)算,葉片結(jié)構(gòu)參數(shù)見表2。

    表2 葉片結(jié)構(gòu)參數(shù)

    數(shù)值模擬計(jì)算時(shí),分別計(jì)算不同葉片結(jié)構(gòu)的旋風(fēng)分離器在不同液滴粒徑下的分離效率和不同進(jìn)氣速度下的壓降。然后,保持單一變量,對(duì)不同葉片結(jié)構(gòu)、相同液滴粒徑和進(jìn)氣速度下的分離效率和壓降計(jì)算結(jié)果進(jìn)行對(duì)比分析,獲得葉片結(jié)構(gòu)對(duì)分離器性能的影響規(guī)律。

    5.1 葉片高度的影響

    圖11示出軸流導(dǎo)葉式旋風(fēng)分離器的分離效率隨葉片高度變化的曲線。由圖可知,葉片高度增加,旋風(fēng)分離器的分離效率呈下降趨勢(shì)。這是由于其它參數(shù)不變,葉片高度的增加使得葉片的螺旋角減小,受葉片導(dǎo)流的旋轉(zhuǎn)氣流的軸向速度將增加,切向速度減小,因而氣流中的液滴的離心力減小,更容易通過分離空間,經(jīng)排氣管逃逸。圖12示出旋風(fēng)分離器的壓降隨葉片高度變化的曲線,由圖可知,葉片高度增加,壓降隨之降低。這是由于葉片高度越小,受葉片導(dǎo)流的氣流在分離空間內(nèi)旋轉(zhuǎn)流動(dòng)的時(shí)間則越短,氣流的速度提升速率更快,氣流與導(dǎo)流葉片的作用力也更大,排氣口的壓力增量卻較小,所以壓降值變大。

    圖12 壓降隨葉片高度變化曲線

    綜合以上葉片高度對(duì)旋風(fēng)分離器性能的影響情況,為保證較高的分離效率,且盡量減小壓降,葉片高度應(yīng)在40~45 mm范圍內(nèi)選擇。

    5.2 葉片圈數(shù)的影響

    圖13示出軸流導(dǎo)葉式旋風(fēng)分離器的分離效率隨葉片圈數(shù)變化的曲線。

    圖13 分離效率隨葉片圈數(shù)變化曲線

    由圖可知,隨著葉片圈數(shù)的增加,分離效率先增大后減小。這是由于葉片圈數(shù)增加,葉片的其它結(jié)構(gòu)參數(shù)不變,則葉片的螺旋角將增大,對(duì)氣流的導(dǎo)流作用將增強(qiáng),氣流中的液滴的旋轉(zhuǎn)角度也將增大,更多的液滴將受離心力作用被分離出來,分離效率增大。然而,當(dāng)葉片圈數(shù)過大,氣流在分離空間的行程增加,運(yùn)動(dòng)氣流對(duì)攜帶液滴的作用力也將增加,液滴更易于被氣流攜帶經(jīng)排氣管逃逸,將不利于液滴的分離。

    圖14示出旋風(fēng)分離器的壓降隨葉片圈數(shù)變化的曲線。由圖可以看出,分離器的壓降隨著葉片圈數(shù)的增加逐漸升高。因?yàn)槿~片圈數(shù)增加,葉片對(duì)氣流的旋流作用增大,葉片的壁面對(duì)氣流的作用力也增加,通過分離空間的氣流的壓降即增大。另外,進(jìn)氣速度增大,氣流的動(dòng)能增加,葉片壁面對(duì)氣流的作用力也增大,壓降也隨之增大。

    圖14 壓降隨葉片圈數(shù)變化曲線

    葉片圈數(shù)對(duì)旋風(fēng)分離器性能的影響情況,為保證較高的分離效率,且較低的壓降,葉片圈數(shù)應(yīng)選擇在0.25~0.30范圍內(nèi)。

    6 結(jié)論

    (1)分離器內(nèi)柱筒和錐筒接合區(qū)段的外旋流具有很高的切向速度,是分離器的高壓區(qū),排氣管入口前后切向速度很低,形成一個(gè)低壓區(qū),是壓降產(chǎn)生的主要區(qū)段。軸流導(dǎo)葉式分離器的分離效率的提高往往伴隨著壓降的增大,應(yīng)用中需要權(quán)衡利弊。

    (2)隨著葉片高度增大,軸流導(dǎo)葉式旋風(fēng)分離器的分離效率逐漸降低,壓降平穩(wěn)減小。軸流導(dǎo)葉式旋風(fēng)分離器的葉片高度建議在40~45 mm內(nèi)選擇。

    (3)隨著葉片圈數(shù)增大,軸流導(dǎo)葉式旋風(fēng)分離器的分離效率先增大后略有減小,壓降一直持續(xù)增大。軸流導(dǎo)葉式旋風(fēng)分離器的葉片圈數(shù)選擇范圍為 0.25~0.30。

    猜你喜歡
    切向速度圈數(shù)軸流
    一起來了解風(fēng)向和風(fēng)速吧
    旋風(fēng)分離器內(nèi)氣相旋轉(zhuǎn)流不穩(wěn)定性的實(shí)驗(yàn)研究
    軸流壓氣機(jī)效率評(píng)定方法
    通過絞車鋼絲繩計(jì)算井深
    雙旋流耦合式旋流反應(yīng)器內(nèi)切向速度分布研究
    晨起轉(zhuǎn)腰改善便秘
    超大型軸流風(fēng)機(jī)在十天高速西秦嶺隧道中的應(yīng)用
    空中顯示計(jì)數(shù)的跳繩
    微型軸流風(fēng)機(jī)低溫失效分析
    高爐用軸流壓縮機(jī)的控制與實(shí)踐
    河南科技(2014年5期)2014-02-27 14:08:33
    亚洲av片天天在线观看| 成年美女黄网站色视频大全免费| 俄罗斯特黄特色一大片| 欧美在线黄色| 两个人免费观看高清视频| 日本欧美视频一区| 老熟女久久久| 欧美97在线视频| 黑丝袜美女国产一区| 69精品国产乱码久久久| 色视频在线一区二区三区| 久久人人爽人人片av| 成年动漫av网址| 欧美日韩亚洲综合一区二区三区_| 欧美日韩视频精品一区| 欧美另类一区| 可以免费在线观看a视频的电影网站| 亚洲第一av免费看| 欧美精品高潮呻吟av久久| 极品少妇高潮喷水抽搐| 国产日韩欧美在线精品| 老汉色av国产亚洲站长工具| 国产精品九九99| av有码第一页| 亚洲第一青青草原| 国产精品久久久人人做人人爽| 久久精品熟女亚洲av麻豆精品| 80岁老熟妇乱子伦牲交| 国产三级黄色录像| 国产精品秋霞免费鲁丝片| www.精华液| 一本久久精品| 午夜福利乱码中文字幕| 久久久久国内视频| 国产高清国产精品国产三级| 黄色视频在线播放观看不卡| 一个人免费看片子| 香蕉丝袜av| 精品久久久久久电影网| 亚洲欧美一区二区三区久久| www.自偷自拍.com| 国产有黄有色有爽视频| 麻豆国产av国片精品| 亚洲国产av影院在线观看| 中文欧美无线码| 免费观看a级毛片全部| 如日韩欧美国产精品一区二区三区| 麻豆av在线久日| 久久99一区二区三区| 亚洲国产欧美一区二区综合| 中国美女看黄片| 亚洲熟女精品中文字幕| 菩萨蛮人人尽说江南好唐韦庄| av福利片在线| 午夜老司机福利片| 人人妻人人澡人人看| 桃花免费在线播放| 中文字幕最新亚洲高清| 国产1区2区3区精品| 亚洲成人免费电影在线观看| 成年人免费黄色播放视频| 制服诱惑二区| 久久影院123| 国产精品99久久99久久久不卡| 精品一区二区三卡| 国产成人精品久久二区二区91| 中文字幕色久视频| 色老头精品视频在线观看| 久久久国产欧美日韩av| 99精国产麻豆久久婷婷| 最近最新免费中文字幕在线| 欧美日韩成人在线一区二区| 久久九九热精品免费| 日韩一卡2卡3卡4卡2021年| 18禁黄网站禁片午夜丰满| 91国产中文字幕| 青草久久国产| 国产欧美日韩一区二区精品| 国产成人一区二区三区免费视频网站| 大片电影免费在线观看免费| 美女国产高潮福利片在线看| 国产av一区二区精品久久| 人成视频在线观看免费观看| 久久国产精品人妻蜜桃| 夜夜骑夜夜射夜夜干| 亚洲午夜精品一区,二区,三区| 一区二区三区乱码不卡18| 欧美人与性动交α欧美软件| 欧美激情久久久久久爽电影 | 男人操女人黄网站| 91精品国产国语对白视频| 纯流量卡能插随身wifi吗| 捣出白浆h1v1| 极品人妻少妇av视频| 久久精品国产综合久久久| 99国产极品粉嫩在线观看| 久久久国产精品麻豆| 免费在线观看黄色视频的| 亚洲 欧美一区二区三区| 中文字幕制服av| 国产欧美亚洲国产| 亚洲欧美成人综合另类久久久| 午夜福利视频精品| 日韩欧美国产一区二区入口| 久久人妻熟女aⅴ| 国产男人的电影天堂91| 亚洲精品第二区| 亚洲第一青青草原| 久久久精品国产亚洲av高清涩受| 高清欧美精品videossex| 少妇精品久久久久久久| 99国产精品一区二区蜜桃av | 桃花免费在线播放| 最近最新免费中文字幕在线| 高清在线国产一区| 在线亚洲精品国产二区图片欧美| 90打野战视频偷拍视频| 老司机午夜福利在线观看视频 | 久热爱精品视频在线9| 少妇被粗大的猛进出69影院| 99精品久久久久人妻精品| 真人做人爱边吃奶动态| 精品人妻1区二区| 中文字幕另类日韩欧美亚洲嫩草| 久久影院123| 999精品在线视频| 欧美日本中文国产一区发布| 两人在一起打扑克的视频| 女人爽到高潮嗷嗷叫在线视频| videos熟女内射| 在线看a的网站| 免费久久久久久久精品成人欧美视频| 欧美成人午夜精品| 国产一区二区三区在线臀色熟女 | 欧美变态另类bdsm刘玥| 欧美xxⅹ黑人| 欧美大码av| 丝袜在线中文字幕| 少妇裸体淫交视频免费看高清 | 国产成人欧美| 中文字幕另类日韩欧美亚洲嫩草| 国产高清videossex| 亚洲三区欧美一区| 女性生殖器流出的白浆| 国产成人欧美在线观看 | 高清在线国产一区| 国产高清视频在线播放一区 | 国产精品免费视频内射| 五月开心婷婷网| 日本一区二区免费在线视频| 欧美国产精品va在线观看不卡| 午夜福利乱码中文字幕| 精品国产乱码久久久久久小说| 精品久久久精品久久久| 国产不卡av网站在线观看| 久久精品亚洲av国产电影网| 美国免费a级毛片| 久久精品熟女亚洲av麻豆精品| 汤姆久久久久久久影院中文字幕| 午夜免费鲁丝| 女人高潮潮喷娇喘18禁视频| 最近最新免费中文字幕在线| 在线观看舔阴道视频| 国产精品久久久av美女十八| 欧美亚洲 丝袜 人妻 在线| av在线播放精品| 国产日韩欧美亚洲二区| 中国美女看黄片| 看免费av毛片| 成年av动漫网址| 国产有黄有色有爽视频| av天堂久久9| 狠狠精品人妻久久久久久综合| 90打野战视频偷拍视频| 多毛熟女@视频| 激情视频va一区二区三区| 亚洲精品粉嫩美女一区| 免费久久久久久久精品成人欧美视频| 亚洲专区中文字幕在线| 黑丝袜美女国产一区| 美女扒开内裤让男人捅视频| 免费在线观看影片大全网站| 国产成人精品久久二区二区91| 丰满迷人的少妇在线观看| 国产高清视频在线播放一区 | 另类精品久久| 免费女性裸体啪啪无遮挡网站| 高清在线国产一区| 天天躁狠狠躁夜夜躁狠狠躁| 大香蕉久久成人网| 人妻人人澡人人爽人人| 亚洲精品美女久久av网站| 亚洲伊人久久精品综合| 亚洲欧美日韩另类电影网站| 精品国产一区二区三区久久久樱花| 亚洲国产中文字幕在线视频| 99久久人妻综合| 高清黄色对白视频在线免费看| 欧美精品av麻豆av| 免费黄频网站在线观看国产| 麻豆国产av国片精品| 日韩免费高清中文字幕av| 欧美老熟妇乱子伦牲交| 精品国产一区二区久久| 欧美日韩国产mv在线观看视频| 伊人久久大香线蕉亚洲五| 黄片小视频在线播放| 日韩欧美免费精品| 国产成人影院久久av| 少妇的丰满在线观看| 午夜免费鲁丝| 国产日韩欧美视频二区| 成年动漫av网址| 欧美精品一区二区大全| av网站在线播放免费| 亚洲va日本ⅴa欧美va伊人久久 | 色94色欧美一区二区| 国产熟女午夜一区二区三区| 精品国产乱子伦一区二区三区 | 桃红色精品国产亚洲av| 国产亚洲午夜精品一区二区久久| 亚洲精品一卡2卡三卡4卡5卡 | 色精品久久人妻99蜜桃| 久久午夜综合久久蜜桃| 久久久久久免费高清国产稀缺| 久久毛片免费看一区二区三区| 中文字幕人妻丝袜制服| 亚洲一区二区三区欧美精品| 亚洲国产av影院在线观看| 久久国产精品人妻蜜桃| 岛国在线观看网站| 精品一区二区三卡| 亚洲欧美一区二区三区黑人| av片东京热男人的天堂| 美女福利国产在线| 18禁观看日本| 午夜福利在线免费观看网站| 精品国产一区二区久久| 国产成人精品久久二区二区免费| 狠狠精品人妻久久久久久综合| 视频区欧美日本亚洲| 亚洲人成电影免费在线| 精品人妻熟女毛片av久久网站| 男女边摸边吃奶| 熟女少妇亚洲综合色aaa.| 亚洲一区中文字幕在线| 久久久久久久国产电影| 精品一品国产午夜福利视频| 亚洲精品中文字幕一二三四区 | 久久精品成人免费网站| 久久久久久亚洲精品国产蜜桃av| 脱女人内裤的视频| 一区二区三区激情视频| 亚洲男人天堂网一区| 日本av免费视频播放| 青青草视频在线视频观看| 蜜桃国产av成人99| 岛国毛片在线播放| av片东京热男人的天堂| 国产成+人综合+亚洲专区| 亚洲国产毛片av蜜桃av| 亚洲七黄色美女视频| 2018国产大陆天天弄谢| 国产精品99久久99久久久不卡| 午夜福利一区二区在线看| 亚洲av国产av综合av卡| 窝窝影院91人妻| av天堂在线播放| 国产有黄有色有爽视频| 日本wwww免费看| 美女福利国产在线| 亚洲av成人不卡在线观看播放网 | 中国国产av一级| 亚洲欧美激情在线| 一区福利在线观看| 亚洲国产精品一区二区三区在线| 日韩大片免费观看网站| 王馨瑶露胸无遮挡在线观看| 中文字幕制服av| 女人久久www免费人成看片| 亚洲欧美成人综合另类久久久| 午夜日韩欧美国产| av在线老鸭窝| 欧美亚洲 丝袜 人妻 在线| 欧美乱码精品一区二区三区| 老汉色av国产亚洲站长工具| 人人澡人人妻人| 精品视频人人做人人爽| 国产一区二区在线观看av| 国产亚洲av高清不卡| 俄罗斯特黄特色一大片| 国产1区2区3区精品| 亚洲av美国av| 黑人猛操日本美女一级片| 国产成人系列免费观看| 精品久久久久久电影网| 久久久久久久大尺度免费视频| 1024视频免费在线观看| 国产老妇伦熟女老妇高清| 99精国产麻豆久久婷婷| 精品少妇内射三级| 欧美激情极品国产一区二区三区| 两个人看的免费小视频| 80岁老熟妇乱子伦牲交| 男女高潮啪啪啪动态图| 精品国产一区二区三区四区第35| 一本色道久久久久久精品综合| 在线亚洲精品国产二区图片欧美| 日本a在线网址| 99久久综合免费| 久久av网站| 久久久国产一区二区| 在线观看www视频免费| 久久久久久久久久久久大奶| 欧美激情久久久久久爽电影 | 麻豆国产av国片精品| 97精品久久久久久久久久精品| 高清黄色对白视频在线免费看| 精品卡一卡二卡四卡免费| 久久av网站| 各种免费的搞黄视频| 纯流量卡能插随身wifi吗| 制服诱惑二区| 丰满饥渴人妻一区二区三| 高清黄色对白视频在线免费看| 亚洲精品国产av成人精品| 在线天堂中文资源库| av一本久久久久| 精品国产一区二区三区四区第35| av又黄又爽大尺度在线免费看| 国产深夜福利视频在线观看| 日本一区二区免费在线视频| 一区二区三区激情视频| 久久九九热精品免费| 久久久久国产精品人妻一区二区| 老鸭窝网址在线观看| 黄色毛片三级朝国网站| 久久精品aⅴ一区二区三区四区| 欧美日韩视频精品一区| 亚洲成人手机| 天天影视国产精品| 欧美亚洲 丝袜 人妻 在线| 淫妇啪啪啪对白视频 | av视频免费观看在线观看| 久久香蕉激情| 国产在视频线精品| 少妇被粗大的猛进出69影院| 视频区图区小说| 一区二区三区乱码不卡18| 久久精品国产亚洲av香蕉五月 | 男人添女人高潮全过程视频| 国产av又大| 又黄又粗又硬又大视频| 久久中文看片网| 亚洲成国产人片在线观看| 亚洲精品第二区| 午夜福利在线免费观看网站| 久久热在线av| 久久av网站| 成人影院久久| 另类亚洲欧美激情| 中文字幕最新亚洲高清| 宅男免费午夜| svipshipincom国产片| 老熟女久久久| 亚洲精品久久久久久婷婷小说| 99国产精品99久久久久| 天天躁狠狠躁夜夜躁狠狠躁| 国产黄色免费在线视频| 看免费av毛片| 午夜福利视频在线观看免费| 色精品久久人妻99蜜桃| 窝窝影院91人妻| 日本欧美视频一区| 国产熟女午夜一区二区三区| 欧美国产精品va在线观看不卡| 一本一本久久a久久精品综合妖精| 国产麻豆69| 久久国产精品影院| av有码第一页| 欧美少妇被猛烈插入视频| 黑人巨大精品欧美一区二区mp4| 老司机影院毛片| 黄片播放在线免费| 在线观看人妻少妇| 久久久久久久久久久久大奶| 精品高清国产在线一区| 欧美午夜高清在线| 桃花免费在线播放| 亚洲国产成人一精品久久久| 最黄视频免费看| www日本在线高清视频| 在线精品无人区一区二区三| 自拍欧美九色日韩亚洲蝌蚪91| 中文字幕制服av| 电影成人av| 丁香六月欧美| 老司机福利观看| 久久久久精品国产欧美久久久 | 久热爱精品视频在线9| 国产精品九九99| 国产精品亚洲av一区麻豆| 国产高清视频在线播放一区 | 黑人巨大精品欧美一区二区mp4| 欧美日韩亚洲综合一区二区三区_| 色视频在线一区二区三区| 欧美+亚洲+日韩+国产| 国产精品 国内视频| 黄色a级毛片大全视频| 久久狼人影院| 各种免费的搞黄视频| 欧美日韩黄片免| 亚洲一码二码三码区别大吗| 久久午夜综合久久蜜桃| 人人妻,人人澡人人爽秒播| 老鸭窝网址在线观看| 91老司机精品| 乱人伦中国视频| 超色免费av| 91精品国产国语对白视频| 国产欧美日韩一区二区三 | 丰满迷人的少妇在线观看| 新久久久久国产一级毛片| 国产精品成人在线| 色老头精品视频在线观看| 成人免费观看视频高清| 熟女少妇亚洲综合色aaa.| 极品少妇高潮喷水抽搐| 69精品国产乱码久久久| 亚洲欧美精品综合一区二区三区| 亚洲九九香蕉| 久久 成人 亚洲| 精品亚洲成国产av| 亚洲一区二区三区欧美精品| 久久国产精品人妻蜜桃| 亚洲精品av麻豆狂野| 欧美变态另类bdsm刘玥| 国产精品亚洲av一区麻豆| 91字幕亚洲| a在线观看视频网站| 欧美国产精品va在线观看不卡| 一二三四社区在线视频社区8| 国产男人的电影天堂91| av天堂在线播放| 亚洲欧美精品自产自拍| 视频区欧美日本亚洲| 美女扒开内裤让男人捅视频| 亚洲少妇的诱惑av| 在线观看www视频免费| 丝袜美足系列| 蜜桃在线观看..| 大片电影免费在线观看免费| 老熟女久久久| 午夜福利在线观看吧| 一区二区三区乱码不卡18| 亚洲久久久国产精品| 亚洲精品国产一区二区精华液| 新久久久久国产一级毛片| 亚洲色图 男人天堂 中文字幕| 97人妻天天添夜夜摸| 日韩欧美免费精品| 久久久久网色| 欧美亚洲 丝袜 人妻 在线| 亚洲国产欧美日韩在线播放| 日本vs欧美在线观看视频| 亚洲第一av免费看| av天堂久久9| 国产亚洲欧美在线一区二区| 久久精品国产综合久久久| 日本a在线网址| 国产视频一区二区在线看| 成人手机av| 久久人人97超碰香蕉20202| 亚洲伊人久久精品综合| 免费少妇av软件| 午夜成年电影在线免费观看| 精品少妇黑人巨大在线播放| 夜夜骑夜夜射夜夜干| 一级片免费观看大全| 婷婷成人精品国产| 女人精品久久久久毛片| kizo精华| 啦啦啦中文免费视频观看日本| 女人久久www免费人成看片| 久久精品亚洲熟妇少妇任你| a在线观看视频网站| 男人舔女人的私密视频| 久久久久久亚洲精品国产蜜桃av| 午夜福利一区二区在线看| av在线播放精品| 人人妻人人澡人人爽人人夜夜| 热99国产精品久久久久久7| 电影成人av| 婷婷色av中文字幕| 国产成人欧美在线观看 | 99国产精品一区二区三区| 亚洲激情五月婷婷啪啪| 午夜视频精品福利| 日韩中文字幕视频在线看片| 少妇被粗大的猛进出69影院| 美女福利国产在线| 大型av网站在线播放| 97人妻天天添夜夜摸| 欧美日韩视频精品一区| 国产精品久久久av美女十八| 青春草亚洲视频在线观看| 女性生殖器流出的白浆| 97人妻天天添夜夜摸| 欧美日韩视频精品一区| 搡老岳熟女国产| 中国国产av一级| 精品久久久久久久毛片微露脸 | 久久人妻福利社区极品人妻图片| 97精品久久久久久久久久精品| 久9热在线精品视频| 国产在线观看jvid| 中文字幕最新亚洲高清| 黄色视频,在线免费观看| 欧美激情久久久久久爽电影 | 美女福利国产在线| 大香蕉久久网| 黄色a级毛片大全视频| 国产欧美日韩精品亚洲av| 纵有疾风起免费观看全集完整版| 国产成人啪精品午夜网站| 国产在线一区二区三区精| 日本av免费视频播放| 国产亚洲一区二区精品| 中亚洲国语对白在线视频| 精品少妇一区二区三区视频日本电影| 伊人亚洲综合成人网| 欧美日韩av久久| 啪啪无遮挡十八禁网站| 久久国产精品影院| 久久精品成人免费网站| 午夜两性在线视频| 欧美日韩视频精品一区| 国产xxxxx性猛交| 精品一区在线观看国产| 人人妻人人澡人人爽人人夜夜| 一级a爱视频在线免费观看| 日本一区二区免费在线视频| 考比视频在线观看| 午夜日韩欧美国产| 国产成人影院久久av| 亚洲 国产 在线| 淫妇啪啪啪对白视频 | 桃花免费在线播放| 国产一卡二卡三卡精品| 99九九在线精品视频| 欧美日韩成人在线一区二区| 新久久久久国产一级毛片| 国产人伦9x9x在线观看| 香蕉丝袜av| 一边摸一边做爽爽视频免费| 亚洲av男天堂| 亚洲第一青青草原| 一本久久精品| 波多野结衣一区麻豆| 国产又爽黄色视频| 啦啦啦视频在线资源免费观看| 成年美女黄网站色视频大全免费| 久久久精品区二区三区| 少妇的丰满在线观看| 国产精品二区激情视频| 欧美性长视频在线观看| 91国产中文字幕| tube8黄色片| 一区二区三区四区激情视频| 男女免费视频国产| 亚洲精品中文字幕在线视频| 亚洲精品国产一区二区精华液| 黄色视频不卡| 欧美激情久久久久久爽电影 | 美国免费a级毛片| 一区二区三区精品91| 免费黄频网站在线观看国产| 淫妇啪啪啪对白视频 | 国产在视频线精品| 国产精品 国内视频| av线在线观看网站| av一本久久久久| 99精国产麻豆久久婷婷| 婷婷色av中文字幕| 人妻久久中文字幕网| 久久久久精品人妻al黑| 亚洲va日本ⅴa欧美va伊人久久 | 久久国产精品影院| 免费不卡黄色视频| 考比视频在线观看| av国产精品久久久久影院| 日韩欧美免费精品| 性色av一级| 色综合欧美亚洲国产小说| 久久精品人人爽人人爽视色| 亚洲欧美日韩另类电影网站| 老司机午夜十八禁免费视频| 久久久欧美国产精品| 午夜福利影视在线免费观看| 欧美日韩视频精品一区| 国产日韩欧美在线精品| 中文字幕av电影在线播放| 国产成人影院久久av| 新久久久久国产一级毛片| 国产日韩一区二区三区精品不卡| 丝袜美腿诱惑在线| 丝袜在线中文字幕| 精品一区在线观看国产| 亚洲性夜色夜夜综合| 999久久久精品免费观看国产| 人妻久久中文字幕网| 久久国产精品大桥未久av| 热re99久久精品国产66热6| 中文字幕色久视频| 久久久久精品人妻al黑|