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

    基于過程模擬軟件的管殼式換熱器優(yōu)化設(shè)計

    2010-12-27 06:01:06李保紅張玨庭
    大連民族大學(xué)學(xué)報 2010年1期
    關(guān)鍵詞:優(yōu)化設(shè)計

    李保紅,張玨庭

    (1.大連民族學(xué)院生命科學(xué)學(xué)院,遼寧大連 116605;2.成功大學(xué)化學(xué)工程系,臺灣臺南 70101)

    基于過程模擬軟件的管殼式換熱器優(yōu)化設(shè)計

    李保紅1,張玨庭2

    (1.大連民族學(xué)院生命科學(xué)學(xué)院,遼寧大連 116605;2.成功大學(xué)化學(xué)工程系,臺灣臺南 70101)

    提出基于國際流行的 ASPEN PLUS模擬軟件,通過與必要的手工計算相結(jié)合,高效地設(shè)計出符合中國相關(guān)標準管殼式換熱器的步驟和方法。復(fù)雜而且繁瑣的能量平衡和壓力降計算由軟件來完成,設(shè)計者只需依照相關(guān)的標準,通過簡單的手工計算確定出離散變量的取值,再基于模擬軟件的計算結(jié)果,在壓力降和標準許可的范圍內(nèi),調(diào)整離散變量的取值以便進一步提高總傳熱系數(shù),從而節(jié)省傳熱面積。通過重新設(shè)計一個文獻實例來演示所提方法的簡單性和有效性,所得的換熱面積比報道值節(jié)省了 66.7%。關(guān)鍵詞:管殼式換熱器;優(yōu)化設(shè)計;過程模擬軟件

    換熱器是一種實現(xiàn)物料之間熱量傳遞的設(shè)備,廣泛應(yīng)用于化工、冶金、電力、食品等行業(yè)。在化工裝置中換熱設(shè)備占設(shè)備數(shù)量的 40%左右,占總投資的 35%~46%[1]。目前,在換熱設(shè)備中,使用量最大的是管殼式換熱器,尤其在高溫、高壓和大型換熱設(shè)備中占有絕對優(yōu)勢[2]。一般來講,管殼式換熱器具有易于加工制造、成本低、可靠性高,且能適應(yīng)高溫高壓的特點。隨著新型高效傳熱管的不斷出現(xiàn),使得管殼式換熱器的應(yīng)用范圍得以不斷擴大,更增添了管殼式換熱器的生命力[3]。

    如何根據(jù)不同的生產(chǎn)工藝條件設(shè)計出投資省、能耗低、傳熱效率高、維修方便的換熱器,是工藝設(shè)計人員重要的工作[4],也是化工類專業(yè)學(xué)生必修的課程設(shè)計項目之一。換熱器的工藝設(shè)計主要包括傳熱和阻力計算兩個方面。由于換熱器的設(shè)計方法比較煩雜[5-7],且需要迭代計算,故借助于日益普及的計算機軟件進行優(yōu)化設(shè)計則可以極大地提高工作效率。

    目前,工程上已大量使用商業(yè)軟件進行換熱器的計算。最著名的專業(yè)換熱器計算軟件主要有成立于 1962年的美國傳熱研究公司 (Heat Transfer Research Inc.,即 HTR I)開發(fā)的 Xchanger Suite軟件;成立于 1967年的英國傳熱及流體服務(wù)中心(Heat Transfer and Fluid Flow Service,即 HTFS)開發(fā)的 HTFS系列軟件[1]和 B-JAC軟件[8]。換熱器計算軟件發(fā)展到今天,在功能上已經(jīng)可以向制造廠商提供設(shè)備條件[4]。

    為了便于組織工業(yè)生產(chǎn),換熱器的設(shè)計要盡可能符合相關(guān)的行業(yè)標準。對于管殼式換熱器,國外主要標準有 TEMA(Tubular ExchangersManufacturers Association)和AS ME(American Society ofMechanical Engineers);國內(nèi)主要標準有國標GB151-1999(管殼式換熱器標準),行業(yè)標準 JB/T 4715-92(固定管板式換熱器形式與基本參數(shù))和 HG 21503-92(鋼制固定式薄管板換熱器)。正是這些標準的制約,使得設(shè)計變量如殼體直徑和管長是非連續(xù)變化的,因而進一步增加了優(yōu)化設(shè)計的難度。

    隨著中國科技與經(jīng)濟實力的不斷增強,愈來愈多的科研單位和高校引進了國際流行的化工過程模擬系統(tǒng)如 ASPEN PLUS、ProII和 HYSYS軟件,這些軟件都具有功能強大的物性計算系統(tǒng)和嚴格的換熱器單元計算模型。但到目前為止,采用這些功能強大的模擬軟件進行換熱器優(yōu)化設(shè)計的研究還鮮有報道[9]。另外,雖然這些軟件所附帶的說明文檔介紹了軟件功能和使用方法[10],但如何充分利用這些功能高效地解決實際問題仍是使用者需要仔細思考的問題。還有從國外引進的軟件所用的默認值與國內(nèi)相關(guān)的標準有差異,如何設(shè)計出符合中國相關(guān)標準的換熱器也是需要解決的重要問題。

    本文選用在中國較為流行的 ASPEN PLUS模擬軟件 (2006.5版)作為管殼式換熱器優(yōu)化設(shè)計的工具,提出高效地設(shè)計出符合中國相關(guān)工業(yè)標準管殼式換熱器的詳細步驟和方法。應(yīng)當指出,由于其他的模擬軟件與 ASPEN PLUS在功能上是相通的,因而本文所提出的設(shè)計步驟和方法也可為使用其他模擬軟件進行換熱器優(yōu)化設(shè)計提供十分有益的借鑒。

    1 問題定義

    已知要被加熱或冷卻的工藝物流的流量、壓力、組成、初始和目標溫度,以及與之匹配換熱的流股組成、初始和目標溫度。要求設(shè)計出符合相關(guān)行業(yè)標準的管殼式換熱器,且該換熱器能夠完成指定的熱交換負荷并符合給定的壓力降要求。優(yōu)化的目標是在完成任務(wù)條件下所需的換熱面積最小。

    2 設(shè)計步驟

    基于研究者的設(shè)計經(jīng)驗和ASPEN PLUS(2006.5版)現(xiàn)有的功能,主要的設(shè)計步驟總結(jié)如下:

    (1)根據(jù)換熱流股所涉及的組分、操作溫度和壓力,選擇適合的物性計算方法 (PropertyMethod)。

    (2)選用 ASPEN PLUS中的“Heater”單元模型(只需輸入一個流股數(shù)據(jù),進行能量平衡計算),輸入工藝物流的相關(guān)數(shù)據(jù),計算出換熱器的負荷 Q。

    (3)選用“HeatX”模型 (換熱器嚴格計算模型)替換“Heater”單元模型,并選用模型中的“Shortcut”計算類型和“Design”模式,以確定匹配熱流股的流量。

    (4)參考相關(guān)的國家和行業(yè)標準,根據(jù)工藝流股與匹配流股的物性以及操作條件,選定殼程與管程流股;選用“HeatX”模型中的“Detailed”計算類型和“Rating”模式,通過手工計算來選取離散變量和軟件運算來進行能量平衡和嚴格壓力降的計算;然后再基于所得結(jié)果進行調(diào)優(yōu)直至滿意為止,調(diào)優(yōu)的原則是在壓力降和標準許可的范圍內(nèi),調(diào)整離散變量的取值以便提高總傳熱系數(shù),從而節(jié)省傳熱面積。

    (5)選用“HeatX”模型中的“Detailed”計算類型和“Simulation”模式進行核算與驗證。

    (6)若上一步所得設(shè)計結(jié)果不符合面積裕度或者壓力降約束,則返回至第 4步。

    3 設(shè)計實例

    3.1 給定條件

    換熱器工藝流程如圖 1,二氟二氯甲烷 (氟里昂 -12)(Freon-12,CCl2F2)作為工藝冷流股,流量是 10 560 kg·h-1,壓力為 7.58×105Pa(絕壓,下同),其溫度需要從 240 K升至 300 K;與之匹配換熱的熱流股是乙二醇 (ethylene glycol,C2H6O2),其初始溫度是 350 K,壓力為 2.02×105Pa。要求設(shè)計管殼式換熱器完成上述任務(wù)。另外,現(xiàn)場工程師推薦乙二醇的出口溫度應(yīng)當至少比冷流股的出口溫度高 10 K,且使用碳鋼傳熱管,要求殼程和管程的壓力降均不超過 6.8×104Pa。優(yōu)化設(shè)計的目標是,在完成給定任務(wù)條件下,所需的換熱器面積最小[9]。

    圖 1 換熱器工藝流程圖

    3.2 設(shè)計過程

    3.2.1 選擇合適的物性計算方法

    由于匹配換熱所涉及的流股組分為極性、非電解質(zhì)且操作壓力小于 1.01×106Pa,故選擇NRTL類模型進行相關(guān)的物性計算[10]。應(yīng)當指出,基于流股 Freon-12數(shù)據(jù),既使采用同屬于NRTL類模型的不同物性計算方法,所得換熱器的加熱負荷也有明顯差異,具體數(shù)據(jù)見表 1,最大偏差為 6%。由此可見,用可靠的物性數(shù)據(jù)或者實驗數(shù)據(jù)來選擇合適的物性計算方法是十分必要的。這里選擇NRTL-HOC僅用于舉例目的。

    表 1 物性計算方法對熱負荷的影響

    3.2.2 確定換熱負荷

    選用“Heater”單元模型,輸入冷流股 Freon-12的相關(guān)數(shù)據(jù),計算出所需的熱負荷為173 840.5W(見表 1)。應(yīng)當指出,僅用圖 1中 Freon-12流股的數(shù)據(jù)是不足以來計算熱負荷的,還需設(shè)定出口壓力或者氣相分率,這里給定的入口到出口的壓力降是 5.05×104Pa。原因是允許的壓力降是6.8×104Pa,若壓力降大于 7.07×104Pa,則該流股在出口處會發(fā)生氣化 (這一點可以由泡點曲線(如圖 2)或者給定出口的氣相分率來確定),意味著換熱過程中涉及相變。這種情況應(yīng)盡可能避免。

    圖 2 流股 FREON-12的泡點曲線

    3.2.3 確定匹配流股流量

    用 HeatX”模型替換“Heater”單元模型,并輸入熱流股相關(guān)數(shù)據(jù),選用“Shortcut”計算類型進行換熱器熱量平衡計算。此時乙二醇的流量待求,故在模型運行前需輸入其流量的估計值 6 083 kg·h-1,對應(yīng)的出口溫度是 308.2 K,不符合給定條件。通過“Design Spec(設(shè)計給定)”模塊,調(diào)整乙二醇的流量使其離開換熱器時的出口溫度正好為 310 K,計算得到的值為 6 342.18 kg·h-1,其他的流股數(shù)據(jù)見表 2。由于這一步僅涉及物料和能量平衡計算,管程和殼程的壓力降不考慮,暫時均視為零。

    表 2 相關(guān)換熱流股數(shù)據(jù)

    3.2.4 換熱器結(jié)構(gòu)設(shè)計

    選用“HeatX”模型中“Detailed”計算方法和“Rating”模式,進行換熱器的設(shè)計。

    首先,需要選定殼程與管程流股?;谏弦徊接嬎愕玫降奈镄詳?shù)據(jù),乙二醇流股的粘度在3 cp以上,而 Freon-12的粘度均小于 1 cp;另外,Freon-12的壓力為 7.58×105Pa,大于乙二醇流股的 2.02×105Pa;故選擇 Freon-12流股走管程,乙二醇流股走殼程[5,7]。

    其次,設(shè)定以下計算方法:(a)傳熱溫差LMTD(即 △tm)校正系數(shù)基于“Geometry”(幾何尺寸)計算;(b)管程與殼程的壓力降均基于“Geometry”計算;(c)總傳熱系數(shù)U基于“Film Coefficient”(膜系數(shù) ) 計算 ,即

    式中,Uo為以管外側(cè)傳熱面積為基準的總傳熱膜系數(shù),hi,ho分別為管內(nèi)和管外傳熱膜系數(shù);rif,rof,rt分別為管內(nèi)、管外污垢熱阻以及管壁熱阻。應(yīng)當注意,傳熱系數(shù)的計算是以換熱管外表面為基準的,因而式 (1)中管內(nèi)膜系數(shù)和污垢熱阻應(yīng)是乘了換熱管的外徑與內(nèi)徑的比 (do/di)之后的校正值;(d)管內(nèi)和管外的分傳熱膜系數(shù)均基于“Geometry”計算,管內(nèi)外側(cè)的污垢熱阻均取 3.52×10-4m2·K·W-1(依據(jù)國標 GB151-1999)。

    最后,借助手工計算確定換熱器的結(jié)構(gòu)尺寸。

    ①計算平均傳熱溫差

    先按純逆流計算

    ②估算傳熱面積

    根據(jù)熱流股是乙二醇,其粘度在 1.0 cp以上,冷流股也是有機物,粘度在 1.0 cp以下,總傳熱系數(shù)的范圍是 150~300 W·m-2·K-1[6]。這里選取 200 W·m-2·K-1作為總傳熱系數(shù),則所需的傳熱面積為

    由于總傳熱溫差校正系數(shù)尚未考慮,故實際傳熱面積取估計值的 120%,即 14.6×1.2=17.5 m2。

    ③管徑和管內(nèi)流速

    由于冷熱流股的體積流量在 100 L·min-1左右,流量較小,按照固定管板式換熱器形式與基本參數(shù) (JB/T 4715-92標準),選 φ19×2的碳鋼換熱管。考慮到第(2)步確定的管程流股壓力降應(yīng)不大于 7.07×104Pa,故管內(nèi)流速選取較低值,為1 m·s-1(正常范圍為 0.5~3 m·s-1[5])。

    ④管程數(shù)和傳熱管總根數(shù)

    流股 Freon-12的平均體積流量

    ⑤查 JB/T 4715-92標準

    根據(jù)換熱面積和管程數(shù)以及換熱管根數(shù),選最接近的換熱器配置。這里選定換熱面積為 17.9 m2,殼體公稱直徑為 325 mm(壁厚 12.5 mm),管長為4.5 m,管子總數(shù)為 68根的換熱器。采用正三角形排列,管心距取 25 mm;采用單弓形折流板,缺口率為殼體內(nèi)徑的 25%,折流板間距取 200 mm。

    ⑥設(shè)定接管尺寸

    由于管程與殼程流股均為液態(tài),選定接管內(nèi)流速為 1 m·s-1,按求所需的接管內(nèi)徑。作為示例,殼程和管程流股的平均體積流量取近似值 100 mL·min-1,得到內(nèi)徑Di=46 mm,故殼程和管程的接管直徑均選取 φ57×3.5的無縫鋼管。

    ⑦數(shù)據(jù)輸入“HeatX”模型

    將以上手工計算的數(shù)據(jù)輸入“HeatX”模型,運行后得到初步設(shè)計結(jié)果,見表 3。表 3中各個符號所代表的意義見符號說明部分。

    表 3 換熱器設(shè)計方案 1的主要結(jié)果

    3.2.5 對設(shè)計結(jié)果進行調(diào)優(yōu)

    由表 3數(shù)據(jù)可以發(fā)現(xiàn),所得設(shè)計方案雖然是可行的,但殼程錯流速度和壓力降偏小 (因為流速的正常范圍是 0.2~1.5 m·s-1[5]),需要進一步增加流速,從而提高殼程傳熱膜系數(shù)以達到降低所需傳熱面積的目的。此任務(wù),可以通過縮小折流板間距和殼體直徑的方式來完成。

    ①折流板間距縮小至 100 mm。即將折流板數(shù)目由 21塊增加至 44塊。所得主要結(jié)果見表 4,可見所得的方案也是可行的。4

    表 4 換熱器設(shè)計方案 2的主要結(jié)果

    ②由表 4數(shù)據(jù)發(fā)現(xiàn),實際的總傳熱系數(shù)為247.5 W·m-2·K-1,比估計值大 25%,導(dǎo)致面積裕度大于 30%,另外殼程流體速度依然偏小,因而需要重新選擇傳熱面積大于 13.5 m2,管程數(shù)為 4,殼體內(nèi)徑比 325 mm小的換熱器?;谝陨蠗l件,查 JB/T 4715-92標準發(fā)現(xiàn),殼體外徑為273 mm,換熱面積為 14.7 m2,管程數(shù)為 2,換熱管總數(shù)為 56的這一組參數(shù)最為接近。但需要將管程數(shù)調(diào)整為 4,以保證管程流體的速度在 1 m·s-1附近,所付出的代價是換熱管總數(shù)需要減少 16根 (按每多一程少 8根計算)用于放置分程隔板。這樣即減少 28.6%的換熱面積;作為補償,將換熱管的長度增加到 6 m,使傳熱面積增加 33.3%,達到基本平衡。采用單弓形折流板,缺口率為殼體內(nèi)徑的 25%,擋板間距取 100 mm,擋板數(shù)目是 59塊。將相關(guān)數(shù)據(jù)輸入后,所得的主要結(jié)果見表 5。

    表 5 換熱器設(shè)計方案 3的主要結(jié)果

    從表 5結(jié)果可以看出,雖然殼程的壓力降和流速都有進一步提升的空間,即進一步縮小折流板間距或者減少殼體的直徑,但由于目前折流板間距已是 JB/T 4715-92標準中的最小值,另外,殼體直徑若減小到下一個檔 (219 mm),則無法排40根 4管程的換熱管,所以調(diào)優(yōu)過程到此為止。3.2.6 設(shè)計結(jié)果驗證及對比

    選用“HeatX”模型中的“Detailed”計算類型和“Simulation”模式進行核算。結(jié)果發(fā)現(xiàn),實際的換熱負荷為 187 152 W,比要求值增加了 7.7%,而且冷流股在出口處發(fā)生了 2%的氣化,這說明設(shè)計合理。

    與文獻 [9]在忽略管內(nèi)、外側(cè)的污垢熱阻,面積裕度為 13.2%的條件下所得的設(shè)計結(jié)果相比,其所需的實際換熱面積為 43 m2,而本文所得結(jié)果為 14.3 m2,即節(jié)省了 66.7%的換熱面積。

    4 結(jié) 論

    提出了基于 ASPEN PLUS(版本 2006.5)模擬軟件和手工計算的有機結(jié)合,進行管殼式換熱器優(yōu)化設(shè)計的步驟和方法。通過對一個文獻實例的設(shè)計,說明了所提方法的簡單性和有效性。設(shè)計過程表明,借助于 ASPEN PLUS模擬軟件多次運行,完全可以高效地設(shè)計符合中國工業(yè)標準的管殼式換熱器。

    文中符號說明:

    A——傳熱面積 (m2);

    d,D——分別是傳熱管和殼體的直徑 (mm);

    h——分傳熱膜系數(shù) (W·m-2·K-1);

    L——換熱管長度 (m);

    Ntt,Ntp——分別是換熱管總數(shù)和每管程的管子數(shù);

    r——傳熱熱阻 (m2·K·W-1);

    Re——雷諾數(shù);

    u——流動速度 (m·s-1);

    U——總傳熱系數(shù) (W·m-2·K-1);

    Q——換熱負荷 (W);

    V——體積流量 (L·min-1);

    ε——溫差校正系數(shù);

    η——面積裕度;

    VP——壓力降 (Pa);

    △tm——對數(shù)平均傳熱溫差 (K);

    上角標:

    in——入口;

    out——出口;

    r——需要 ;

    s——實際提供;

    下角標:

    i——管內(nèi) ;

    o——管外;

    s——殼程 ;

    t——管程。

    [1]秦叔經(jīng),葉文邦.換熱器[M].北京:化學(xué)工業(yè)出版社,2003.

    [2]王元文.管殼式換熱器優(yōu)化設(shè)計[J].廣東化工,2005(3):43-44.

    [3]錢頌文.換熱器設(shè)計手冊[M].北京:化學(xué)工業(yè)出版社,2002.

    [4]中國石化集團上海有限公司.化工工藝設(shè)計手冊[M].3版.北京:化學(xué)工業(yè)出版社,2003.

    [5]大連理工大學(xué)化工原理教研室.化工原理課程設(shè)計[M].大連:大連理工大學(xué)出版社,1994.

    [6]潘國昌,郭慶豐.化工設(shè)備設(shè)計[M].北京:清華大學(xué)出版社,1996.

    [7]賈紹義,柴誠敬.化工原理課程設(shè)計 [M].天津:天津大學(xué)出版社,2002.

    [8]LEONG K C,TOH K C,LEONG Y C.Shell and tube heat exchanger design sof tware for educational applications[J].Int.J.Eng.Ed.1998,14(3):217-224.

    [9]LANG J.Design procedure for a heat exchanger on the AspenPlus Software[DB/OL](2009-08-18).http://cbe.sdsmt.edu/nsfproj/aspen/no-ph-cgehex.pdf,2000.

    [10]Asoeb Plus 12.1 User guide[M].Cambridge,USA:Aspen Technology,2003.

    OptimalDesign of Shell-and-tube Heat Exchangers Based on Process Simulation Software

    L I Bao-hong1,ZHANG Jue-ting2
    (1.College ofLife Science,Dalian NationalitiesUniversity,Dalian Liaoning 116605,China;
    2.Depar tment of Chemical Engineering,Cheng KungUniversity,Tainan Taiwan 70101,China)

    We proposed procedures and a method for efficiently designing shell-and-tube heat exchangers compliantwith China's national standards based on Aspen Plus processmodeling tool and in combination with necessary manual calculation.Complex and tedious energy balancing and pressure-drop calculation are performed by software.The designer need only establish the values of discrete variables through simple manual calculation under relevant standards.He/she may then adjust those valueswithin the range of pressure drops and allowances specified by the standards based on the calculation result of the simulation software,thereby further increasing the total coefficient of heat transfer and hence reducing the area of heat transfer.We demonstrated the simplicity and effectiveness of our proposed method bymaking an example design in the literature.The resulting area of heat transferwas 66.7%s maller than the reported value.

    shell-and-tube heat exchanger;optimal design;process simulation software

    TQ015.9

    A

    1009-315X(2010)01-0001-05

    2009-09-25

    國家自然科學(xué)基金資助項目 (20806015)。

    李保紅 (1973-),男,陜西武功人,副教授,博士,主要從事化工過程模擬與優(yōu)化研究。

    (責任編輯 鄒永紅)

    猜你喜歡
    優(yōu)化設(shè)計
    導(dǎo)彈舵面的復(fù)合材料設(shè)計與分析
    航空兵器(2016年4期)2016-11-28 21:47:29
    礦井主排水系統(tǒng)的優(yōu)化設(shè)計與改造
    科技資訊(2016年19期)2016-11-15 08:34:13
    數(shù)據(jù)挖掘?qū)虒W(xué)管理的優(yōu)化設(shè)計
    如何實現(xiàn)小學(xué)數(shù)學(xué)課堂練習(xí)設(shè)計優(yōu)化
    淺析人機工程學(xué)在家具創(chuàng)作中的作用
    試析機械結(jié)構(gòu)優(yōu)化設(shè)計的應(yīng)用及趨勢
    汽車行李箱蓋鉸鏈機構(gòu)的分析及優(yōu)化
    東林煤礦保護層開采卸壓瓦斯抽采優(yōu)化設(shè)計
    橋式起重機主梁結(jié)構(gòu)分析和優(yōu)化設(shè)計
    對無線傳感器網(wǎng)絡(luò)MAC層協(xié)議優(yōu)化的研究與設(shè)計
    科技視界(2016年22期)2016-10-18 15:25:08
    日韩大片免费观看网站| 国产精品久久久久成人av| 日韩免费高清中文字幕av| 日本撒尿小便嘘嘘汇集6| 久久中文字幕一级| a级片在线免费高清观看视频| 日本欧美视频一区| 国产一级毛片在线| 在线观看www视频免费| 亚洲中文字幕日韩| 亚洲国产欧美日韩在线播放| 激情视频va一区二区三区| 热re99久久精品国产66热6| 国产一级毛片在线| 欧美 日韩 精品 国产| 蜜桃国产av成人99| www日本在线高清视频| 欧美人与性动交α欧美软件| a 毛片基地| 99九九在线精品视频| 一级毛片女人18水好多| 在线观看免费午夜福利视频| 色精品久久人妻99蜜桃| 亚洲国产日韩一区二区| 亚洲av美国av| 美女高潮到喷水免费观看| 日本黄色日本黄色录像| 一本久久精品| 三级毛片av免费| 91国产中文字幕| 日本av免费视频播放| 亚洲色图 男人天堂 中文字幕| 91国产中文字幕| tocl精华| av福利片在线| 国产精品成人在线| 青草久久国产| 欧美日韩视频精品一区| 久9热在线精品视频| 国产国语露脸激情在线看| 久久久精品国产亚洲av高清涩受| 成人三级做爰电影| 在线观看一区二区三区激情| 又紧又爽又黄一区二区| 狂野欧美激情性bbbbbb| 99九九在线精品视频| 亚洲自偷自拍图片 自拍| 老司机亚洲免费影院| netflix在线观看网站| 午夜福利影视在线免费观看| 亚洲欧美清纯卡通| 欧美 亚洲 国产 日韩一| 美国免费a级毛片| 99热全是精品| av线在线观看网站| 一个人免费看片子| 亚洲视频免费观看视频| 久久久久精品人妻al黑| 在线天堂中文资源库| 天天添夜夜摸| 99国产精品一区二区蜜桃av | 午夜成年电影在线免费观看| av一本久久久久| 亚洲 欧美一区二区三区| 美国免费a级毛片| 亚洲欧美日韩另类电影网站| 婷婷色av中文字幕| 乱人伦中国视频| 日本欧美视频一区| 久久久欧美国产精品| 久久久水蜜桃国产精品网| 日韩熟女老妇一区二区性免费视频| 国产免费福利视频在线观看| 国产成人精品久久二区二区91| 视频在线观看一区二区三区| 最黄视频免费看| 亚洲成国产人片在线观看| 成人三级做爰电影| 精品第一国产精品| 真人做人爱边吃奶动态| 搡老熟女国产l中国老女人| 在线十欧美十亚洲十日本专区| 黑人巨大精品欧美一区二区mp4| 亚洲精品国产区一区二| 久久人人97超碰香蕉20202| 色老头精品视频在线观看| 欧美精品啪啪一区二区三区 | 亚洲av成人不卡在线观看播放网 | 日韩制服骚丝袜av| 一本—道久久a久久精品蜜桃钙片| 免费观看a级毛片全部| 色精品久久人妻99蜜桃| 又大又爽又粗| 超色免费av| 国产一区二区三区av在线| 亚洲国产欧美一区二区综合| 后天国语完整版免费观看| 久久久国产一区二区| 亚洲av电影在线观看一区二区三区| 久久 成人 亚洲| 老司机深夜福利视频在线观看 | 最新的欧美精品一区二区| 欧美精品人与动牲交sv欧美| 精品久久久久久久毛片微露脸 | www.精华液| 淫妇啪啪啪对白视频 | 十八禁高潮呻吟视频| 最黄视频免费看| 免费看十八禁软件| 国产av又大| 男女边摸边吃奶| 欧美精品一区二区免费开放| 国产精品秋霞免费鲁丝片| 夜夜夜夜夜久久久久| 午夜激情久久久久久久| 久久热在线av| 免费日韩欧美在线观看| 婷婷色av中文字幕| 免费少妇av软件| 国产一区二区三区综合在线观看| 亚洲七黄色美女视频| 91国产中文字幕| 桃花免费在线播放| 亚洲国产精品一区三区| 亚洲欧美清纯卡通| 啦啦啦啦在线视频资源| 日韩欧美国产一区二区入口| 国产成人系列免费观看| 十八禁网站免费在线| 黄色视频,在线免费观看| 大型av网站在线播放| 99精品久久久久人妻精品| 一级片'在线观看视频| 久久影院123| 丁香六月欧美| 老司机在亚洲福利影院| 老汉色∧v一级毛片| 巨乳人妻的诱惑在线观看| 亚洲成av片中文字幕在线观看| 99国产综合亚洲精品| 久久性视频一级片| 亚洲精品乱久久久久久| 啦啦啦中文免费视频观看日本| 五月天丁香电影| 在线亚洲精品国产二区图片欧美| 国产精品麻豆人妻色哟哟久久| 高清视频免费观看一区二区| 天堂中文最新版在线下载| 欧美日韩福利视频一区二区| 一区二区av电影网| 50天的宝宝边吃奶边哭怎么回事| 亚洲va日本ⅴa欧美va伊人久久 | 亚洲精品一二三| 久久精品人人爽人人爽视色| 亚洲一区二区三区欧美精品| 精品一区二区三区四区五区乱码| 青青草视频在线视频观看| 狂野欧美激情性bbbbbb| 天天影视国产精品| 亚洲第一青青草原| 亚洲九九香蕉| 精品熟女少妇八av免费久了| 午夜免费鲁丝| 老司机影院成人| 精品人妻一区二区三区麻豆| 日韩精品免费视频一区二区三区| 女人久久www免费人成看片| 汤姆久久久久久久影院中文字幕| 女人精品久久久久毛片| 制服人妻中文乱码| 王馨瑶露胸无遮挡在线观看| 亚洲专区国产一区二区| 国产精品麻豆人妻色哟哟久久| 亚洲成av片中文字幕在线观看| 人成视频在线观看免费观看| 久久香蕉激情| 国产人伦9x9x在线观看| www日本在线高清视频| 人人妻,人人澡人人爽秒播| 99re6热这里在线精品视频| 天堂中文最新版在线下载| 三级毛片av免费| 国产日韩一区二区三区精品不卡| 亚洲精品国产av成人精品| 欧美成人午夜精品| 久热这里只有精品99| 国产福利在线免费观看视频| 精品一区二区三区四区五区乱码| 黄色片一级片一级黄色片| 亚洲国产av影院在线观看| √禁漫天堂资源中文www| 久久毛片免费看一区二区三区| 亚洲av美国av| 国产高清videossex| 操美女的视频在线观看| 青春草亚洲视频在线观看| 午夜91福利影院| 伊人久久大香线蕉亚洲五| 人妻一区二区av| 亚洲精品久久久久久婷婷小说| 免费av中文字幕在线| 不卡av一区二区三区| 丝袜在线中文字幕| 午夜福利,免费看| 中亚洲国语对白在线视频| 在线av久久热| 亚洲 欧美一区二区三区| 久热爱精品视频在线9| 国产日韩欧美在线精品| 天天添夜夜摸| 极品人妻少妇av视频| 无限看片的www在线观看| 精品少妇内射三级| 一二三四在线观看免费中文在| 日本91视频免费播放| 美女高潮到喷水免费观看| 亚洲天堂av无毛| www.熟女人妻精品国产| 亚洲va日本ⅴa欧美va伊人久久 | 中文字幕人妻丝袜一区二区| 丝袜人妻中文字幕| 狂野欧美激情性xxxx| xxxhd国产人妻xxx| 一边摸一边做爽爽视频免费| 国产精品影院久久| 久久久国产精品麻豆| 男人操女人黄网站| 久久国产精品人妻蜜桃| 在线 av 中文字幕| 精品人妻一区二区三区麻豆| 亚洲精品粉嫩美女一区| 国产熟女午夜一区二区三区| 久久影院123| 男人爽女人下面视频在线观看| 99国产极品粉嫩在线观看| 国产亚洲欧美在线一区二区| 久久久国产成人免费| 亚洲av电影在线观看一区二区三区| 欧美日韩精品网址| 欧美97在线视频| 久久99一区二区三区| 久久久久视频综合| 一本久久精品| 韩国精品一区二区三区| 日韩精品免费视频一区二区三区| 久久精品国产亚洲av高清一级| 国产区一区二久久| 99久久综合免费| 丝瓜视频免费看黄片| videosex国产| 欧美日韩一级在线毛片| 久久影院123| 婷婷成人精品国产| 人人妻人人爽人人添夜夜欢视频| 桃红色精品国产亚洲av| 欧美大码av| 久久国产精品影院| 80岁老熟妇乱子伦牲交| 最近最新中文字幕大全免费视频| 熟女少妇亚洲综合色aaa.| 又黄又粗又硬又大视频| 亚洲少妇的诱惑av| 国产男女超爽视频在线观看| 国产欧美日韩综合在线一区二区| 黄色a级毛片大全视频| 老汉色∧v一级毛片| 91精品三级在线观看| 最黄视频免费看| 大片免费播放器 马上看| 韩国精品一区二区三区| 久久精品人人爽人人爽视色| 久久精品久久久久久噜噜老黄| 黄片播放在线免费| 婷婷成人精品国产| 久久久水蜜桃国产精品网| 国产熟女午夜一区二区三区| 啦啦啦在线免费观看视频4| 男人舔女人的私密视频| 亚洲av电影在线进入| 色婷婷av一区二区三区视频| 免费在线观看视频国产中文字幕亚洲 | 9色porny在线观看| 不卡一级毛片| 日韩,欧美,国产一区二区三区| 最黄视频免费看| 国产免费一区二区三区四区乱码| 嫩草影视91久久| 成年人免费黄色播放视频| 后天国语完整版免费观看| 丁香六月欧美| 嫩草影视91久久| 国产免费av片在线观看野外av| 99精品欧美一区二区三区四区| 超色免费av| 99久久国产精品久久久| 亚洲国产欧美网| 亚洲国产毛片av蜜桃av| cao死你这个sao货| 欧美久久黑人一区二区| 妹子高潮喷水视频| 精品少妇黑人巨大在线播放| 一进一出抽搐动态| 99九九在线精品视频| tube8黄色片| 欧美精品一区二区大全| www.精华液| 秋霞在线观看毛片| 美女高潮喷水抽搐中文字幕| 国产成人av教育| 久久精品亚洲熟妇少妇任你| 欧美变态另类bdsm刘玥| 亚洲一卡2卡3卡4卡5卡精品中文| 精品高清国产在线一区| 国产亚洲av片在线观看秒播厂| 99国产精品一区二区蜜桃av | 99re6热这里在线精品视频| 国产成人一区二区三区免费视频网站| 一本久久精品| h视频一区二区三区| 丝袜人妻中文字幕| 久久99热这里只频精品6学生| 丝瓜视频免费看黄片| 男女下面插进去视频免费观看| 亚洲国产精品一区三区| 成年av动漫网址| 久久中文字幕一级| 法律面前人人平等表现在哪些方面 | 精品亚洲成a人片在线观看| 99国产综合亚洲精品| 夫妻午夜视频| 欧美 日韩 精品 国产| 俄罗斯特黄特色一大片| 久久久国产一区二区| 国产精品国产三级国产专区5o| 亚洲久久久国产精品| 亚洲国产av新网站| 一区二区三区乱码不卡18| 国产精品久久久人人做人人爽| 亚洲欧美精品综合一区二区三区| 亚洲国产看品久久| 国产伦人伦偷精品视频| 操美女的视频在线观看| 下体分泌物呈黄色| 午夜免费观看性视频| 久久久国产成人免费| 欧美乱码精品一区二区三区| 一个人免费看片子| 亚洲精品日韩在线中文字幕| 午夜福利视频精品| 97在线人人人人妻| av网站免费在线观看视频| 久久久久久久精品精品| videos熟女内射| 大型av网站在线播放| 久久久精品94久久精品| 亚洲中文字幕日韩| 十八禁网站网址无遮挡| 色婷婷av一区二区三区视频| 精品福利观看| 欧美黑人欧美精品刺激| 十分钟在线观看高清视频www| 久久人人爽人人片av| 精品乱码久久久久久99久播| 菩萨蛮人人尽说江南好唐韦庄| 欧美日韩av久久| 日韩三级视频一区二区三区| 少妇的丰满在线观看| 久热爱精品视频在线9| 亚洲av日韩精品久久久久久密| 老司机福利观看| 老司机靠b影院| 亚洲av日韩在线播放| 天堂8中文在线网| 伊人久久大香线蕉亚洲五| 久久精品国产亚洲av高清一级| 国产成人av激情在线播放| 亚洲av成人不卡在线观看播放网 | 真人做人爱边吃奶动态| av在线老鸭窝| 这个男人来自地球电影免费观看| 国产99久久九九免费精品| 纯流量卡能插随身wifi吗| 亚洲黑人精品在线| 99九九在线精品视频| 另类亚洲欧美激情| 久久久精品国产亚洲av高清涩受| 欧美性长视频在线观看| 啦啦啦 在线观看视频| svipshipincom国产片| 少妇粗大呻吟视频| 日韩三级视频一区二区三区| 成人国产av品久久久| 在线观看免费高清a一片| 久久精品熟女亚洲av麻豆精品| 亚洲色图 男人天堂 中文字幕| 一区二区三区乱码不卡18| 欧美午夜高清在线| 午夜老司机福利片| 黄色视频在线播放观看不卡| 日本精品一区二区三区蜜桃| 免费不卡黄色视频| 午夜激情av网站| 一个人免费在线观看的高清视频 | 亚洲精品成人av观看孕妇| 亚洲av电影在线进入| 欧美少妇被猛烈插入视频| 纯流量卡能插随身wifi吗| 国产成人精品无人区| 一级毛片精品| 国产免费av片在线观看野外av| 极品人妻少妇av视频| 亚洲熟女毛片儿| 久久人妻福利社区极品人妻图片| 欧美大码av| 91成人精品电影| 亚洲美女黄色视频免费看| 久久精品亚洲熟妇少妇任你| 十八禁网站免费在线| 亚洲第一欧美日韩一区二区三区 | 欧美乱码精品一区二区三区| 亚洲第一青青草原| 人人妻人人澡人人看| 水蜜桃什么品种好| 久久久久精品人妻al黑| 免费久久久久久久精品成人欧美视频| www.自偷自拍.com| 大香蕉久久成人网| 无限看片的www在线观看| 深夜精品福利| 日韩中文字幕视频在线看片| 久久久久久人人人人人| cao死你这个sao货| 91精品国产国语对白视频| 国产91精品成人一区二区三区 | 狂野欧美激情性xxxx| 97精品久久久久久久久久精品| 人人妻人人澡人人爽人人夜夜| 亚洲avbb在线观看| 欧美日韩亚洲高清精品| 中国美女看黄片| 91精品国产国语对白视频| 久久精品国产亚洲av香蕉五月 | 两个人免费观看高清视频| 日韩一卡2卡3卡4卡2021年| 午夜福利影视在线免费观看| 美女高潮喷水抽搐中文字幕| 亚洲avbb在线观看| av天堂在线播放| 精品熟女少妇八av免费久了| 亚洲精品一卡2卡三卡4卡5卡 | 亚洲av电影在线进入| 成人亚洲精品一区在线观看| 亚洲精品日韩在线中文字幕| 欧美人与性动交α欧美软件| 日本精品一区二区三区蜜桃| 午夜久久久在线观看| 精品人妻一区二区三区麻豆| 男人操女人黄网站| 精品少妇一区二区三区视频日本电影| 免费女性裸体啪啪无遮挡网站| 大片免费播放器 马上看| 国产一区二区三区综合在线观看| 国产高清videossex| 巨乳人妻的诱惑在线观看| 日本wwww免费看| 欧美亚洲日本最大视频资源| av超薄肉色丝袜交足视频| 高清欧美精品videossex| 一二三四社区在线视频社区8| www.精华液| 国产精品麻豆人妻色哟哟久久| 国产极品粉嫩免费观看在线| 日韩一区二区三区影片| 纵有疾风起免费观看全集完整版| av电影中文网址| 亚洲精品国产精品久久久不卡| 国产日韩一区二区三区精品不卡| 男女无遮挡免费网站观看| 久久久久久亚洲精品国产蜜桃av| 婷婷成人精品国产| 欧美日韩成人在线一区二区| 亚洲国产精品成人久久小说| 中文字幕另类日韩欧美亚洲嫩草| 91精品国产国语对白视频| 一二三四社区在线视频社区8| 久久久久精品国产欧美久久久 | 亚洲七黄色美女视频| 午夜福利影视在线免费观看| 五月开心婷婷网| 久久午夜综合久久蜜桃| 久热这里只有精品99| 成人三级做爰电影| 丝袜喷水一区| 操美女的视频在线观看| 丝瓜视频免费看黄片| 久久精品久久久久久噜噜老黄| 国产一级毛片在线| 91国产中文字幕| 亚洲精品成人av观看孕妇| 午夜福利在线观看吧| 在线 av 中文字幕| 最近中文字幕2019免费版| 国产成人精品久久二区二区免费| 亚洲精品美女久久av网站| 国产高清国产精品国产三级| 日韩精品免费视频一区二区三区| 爱豆传媒免费全集在线观看| av线在线观看网站| 精品一区二区三区四区五区乱码| 亚洲专区中文字幕在线| 国产熟女午夜一区二区三区| 亚洲精品在线美女| 亚洲精品久久久久久婷婷小说| 日韩 亚洲 欧美在线| 99国产精品一区二区蜜桃av | 免费人妻精品一区二区三区视频| av免费在线观看网站| 欧美在线一区亚洲| 男人爽女人下面视频在线观看| 狠狠狠狠99中文字幕| 欧美 亚洲 国产 日韩一| 少妇裸体淫交视频免费看高清 | av超薄肉色丝袜交足视频| 成年人午夜在线观看视频| 久久人妻熟女aⅴ| 两性午夜刺激爽爽歪歪视频在线观看 | 啦啦啦中文免费视频观看日本| 亚洲第一青青草原| 久久久久久久久免费视频了| 高清视频免费观看一区二区| 看免费av毛片| 成年人免费黄色播放视频| 亚洲精品av麻豆狂野| 人妻一区二区av| 国产高清videossex| 黄片小视频在线播放| 制服人妻中文乱码| 人妻一区二区av| 日本av手机在线免费观看| www.999成人在线观看| 黄色a级毛片大全视频| 亚洲精品粉嫩美女一区| 国产免费一区二区三区四区乱码| 成人亚洲精品一区在线观看| 天堂8中文在线网| 亚洲精品美女久久久久99蜜臀| 成人18禁高潮啪啪吃奶动态图| 18禁国产床啪视频网站| 两个人看的免费小视频| 亚洲成人国产一区在线观看| 青春草亚洲视频在线观看| 午夜福利,免费看| 欧美久久黑人一区二区| 欧美老熟妇乱子伦牲交| 亚洲情色 制服丝袜| 人人妻人人爽人人添夜夜欢视频| 91老司机精品| 日本黄色日本黄色录像| 日日摸夜夜添夜夜添小说| 欧美激情极品国产一区二区三区| 一区二区三区四区激情视频| 国产精品.久久久| 丝袜美足系列| 国产av精品麻豆| 久久精品亚洲熟妇少妇任你| 啪啪无遮挡十八禁网站| 欧美日韩av久久| 中亚洲国语对白在线视频| 日韩一区二区三区影片| 久久精品亚洲av国产电影网| 日韩中文字幕欧美一区二区| 久久精品亚洲av国产电影网| 女人被躁到高潮嗷嗷叫费观| 最黄视频免费看| 午夜免费成人在线视频| 国产成人欧美在线观看 | 9色porny在线观看| 十八禁网站网址无遮挡| 一级毛片精品| 国产男女超爽视频在线观看| 又黄又粗又硬又大视频| 99国产精品一区二区蜜桃av | 欧美日韩福利视频一区二区| 韩国精品一区二区三区| 成人国产av品久久久| 亚洲全国av大片| 少妇人妻久久综合中文| 国产精品久久久av美女十八| 色综合欧美亚洲国产小说| 少妇被粗大的猛进出69影院| 国产高清国产精品国产三级| 91精品三级在线观看| 丝袜美足系列| 性色av乱码一区二区三区2| 一本大道久久a久久精品| 美女大奶头黄色视频| 激情视频va一区二区三区| 国产国语露脸激情在线看| 免费黄频网站在线观看国产| 亚洲精品国产一区二区精华液| 国产又色又爽无遮挡免| 啦啦啦中文免费视频观看日本| 黄色视频,在线免费观看| 亚洲午夜精品一区,二区,三区| 亚洲黑人精品在线| 韩国高清视频一区二区三区| 日韩有码中文字幕| 男女午夜视频在线观看| 久久久国产一区二区| 久久av网站| 中文字幕高清在线视频| 中文字幕人妻丝袜制服| 国产精品99久久99久久久不卡| 桃红色精品国产亚洲av| 天天躁夜夜躁狠狠躁躁|