• <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在线| 国产成人精品久久久久久| freevideosex欧美| 777米奇影视久久| 在线观看美女被高潮喷水网站| 9191精品国产免费久久| 免费日韩欧美在线观看| 中文字幕制服av| 在线观看免费视频网站a站| 黄色配什么色好看| 香蕉丝袜av| 超色免费av| 成人国产av品久久久| 精品久久蜜臀av无| 亚洲欧美日韩卡通动漫| 热re99久久精品国产66热6| 久久精品国产a三级三级三级| 成年女人在线观看亚洲视频| 国产熟女欧美一区二区| 亚洲熟女精品中文字幕| 成人国语在线视频| 国内精品宾馆在线| 亚洲精品一区蜜桃| 黄色一级大片看看| 免费人妻精品一区二区三区视频| 久久久国产一区二区| 精品少妇内射三级| 久久久久精品性色| 人妻 亚洲 视频| 成人亚洲精品一区在线观看| 免费高清在线观看视频在线观看| 国产 一区精品| 97在线人人人人妻| 国产片特级美女逼逼视频| 美女福利国产在线| 精品一区在线观看国产| 777米奇影视久久| 一区二区av电影网| 亚洲伊人久久精品综合| 日本wwww免费看| 色吧在线观看| freevideosex欧美| 天堂中文最新版在线下载| 国产一区二区激情短视频 | 一级片免费观看大全| 精品国产露脸久久av麻豆| 国产一区二区激情短视频 | 国产精品一区二区在线不卡| 久久ye,这里只有精品| 精品久久久久久电影网| 高清av免费在线| 2018国产大陆天天弄谢| 桃花免费在线播放| 我的女老师完整版在线观看| 日韩三级伦理在线观看| 精品午夜福利在线看| 久久99蜜桃精品久久| 国产精品一国产av| 夫妻性生交免费视频一级片| 中国美白少妇内射xxxbb| 一边摸一边做爽爽视频免费| 一本久久精品| 91精品国产国语对白视频| 国产精品一区二区在线不卡| 夫妻午夜视频| 在线观看三级黄色| 色吧在线观看| 国产一区亚洲一区在线观看| 一区二区av电影网| 欧美bdsm另类| av福利片在线| 人人妻人人澡人人看| 亚洲成国产人片在线观看| 国国产精品蜜臀av免费| av.在线天堂| 久久久久久久国产电影| 国产成人欧美| 99热这里只有是精品在线观看| 国产69精品久久久久777片| 精品人妻一区二区三区麻豆| 午夜免费观看性视频| 热99久久久久精品小说推荐| 午夜影院在线不卡| 天美传媒精品一区二区| 色婷婷av一区二区三区视频| 日日撸夜夜添| av播播在线观看一区| 色94色欧美一区二区| 成年人免费黄色播放视频| 国产又爽黄色视频| 久久精品国产鲁丝片午夜精品| 99热这里只有是精品在线观看| 熟女人妻精品中文字幕| 久久99热6这里只有精品| 国产乱人偷精品视频| 乱人伦中国视频| 精品人妻一区二区三区麻豆| 日韩视频在线欧美| 国产深夜福利视频在线观看| 精品久久久久久电影网| 大码成人一级视频| 国产亚洲最大av| 欧美另类一区| 91在线精品国自产拍蜜月| www.熟女人妻精品国产 | 久久久久久久久久久免费av| 人妻 亚洲 视频| 国产成人免费观看mmmm| xxx大片免费视频| 97在线视频观看| √禁漫天堂资源中文www| 欧美另类一区| 国产1区2区3区精品| 国产精品免费大片| 亚洲色图综合在线观看| 91精品伊人久久大香线蕉| 日本色播在线视频| 五月开心婷婷网| 99精国产麻豆久久婷婷| 国产成人精品久久久久久| 亚洲av电影在线进入| 久久人妻熟女aⅴ| 免费人妻精品一区二区三区视频| 少妇人妻久久综合中文| 一边亲一边摸免费视频| 美女内射精品一级片tv| 香蕉精品网在线| 日韩,欧美,国产一区二区三区| 黄片无遮挡物在线观看| 国产片内射在线| 中文字幕精品免费在线观看视频 | 在线精品无人区一区二区三| 欧美精品人与动牲交sv欧美| 男人操女人黄网站| 看非洲黑人一级黄片| 免费av中文字幕在线| 精品一区在线观看国产| 最近中文字幕高清免费大全6| 美女内射精品一级片tv| a级毛片黄视频| 久久韩国三级中文字幕| 国产精品一区二区在线观看99| 成人亚洲欧美一区二区av| 欧美日韩综合久久久久久| 日韩熟女老妇一区二区性免费视频| av有码第一页| 老司机影院毛片| 1024视频免费在线观看| 精品福利永久在线观看| 成人毛片60女人毛片免费| 2022亚洲国产成人精品| 欧美日韩成人在线一区二区| 人人妻人人澡人人看| 如日韩欧美国产精品一区二区三区| 18在线观看网站| 少妇人妻 视频| 午夜久久久在线观看| 久久国产亚洲av麻豆专区| 成年人午夜在线观看视频| 青春草视频在线免费观看| 成人毛片60女人毛片免费| 2018国产大陆天天弄谢| 黄色视频在线播放观看不卡| 人人妻人人爽人人添夜夜欢视频| 国产欧美日韩综合在线一区二区| 毛片一级片免费看久久久久| 成人18禁高潮啪啪吃奶动态图| 国产亚洲最大av| 少妇熟女欧美另类| 欧美日本中文国产一区发布| 久久这里只有精品19| 午夜91福利影院| 精品视频人人做人人爽| 菩萨蛮人人尽说江南好唐韦庄| videossex国产| 亚洲国产精品国产精品| 熟女电影av网| 男女边吃奶边做爰视频| 久久精品夜色国产| 免费在线观看完整版高清| 国产一区二区三区av在线| 校园人妻丝袜中文字幕| 99久久人妻综合| 老司机影院毛片| 国产精品国产三级国产av玫瑰| 一区二区三区乱码不卡18| 久久久久久久亚洲中文字幕| 日韩熟女老妇一区二区性免费视频| 免费在线观看黄色视频的| 亚洲高清免费不卡视频| 国产伦理片在线播放av一区| 午夜激情久久久久久久| 国产深夜福利视频在线观看| 亚洲美女黄色视频免费看| 精品久久久精品久久久| 日本黄色日本黄色录像| 成人影院久久| 久久久久久久大尺度免费视频| 全区人妻精品视频| 精品一区在线观看国产| 欧美亚洲日本最大视频资源| 看十八女毛片水多多多| 女人精品久久久久毛片| 下体分泌物呈黄色| 赤兔流量卡办理| 欧美精品一区二区免费开放| av在线播放精品| 日韩制服丝袜自拍偷拍| 你懂的网址亚洲精品在线观看| 国产亚洲欧美精品永久| 少妇人妻精品综合一区二区| 久久精品国产亚洲av涩爱| 99热国产这里只有精品6| 大片免费播放器 马上看| 精品亚洲成a人片在线观看| 寂寞人妻少妇视频99o| 多毛熟女@视频| 热re99久久精品国产66热6| 国产一区二区激情短视频 | 中国美白少妇内射xxxbb| 欧美精品一区二区免费开放| 午夜福利影视在线免费观看| 日韩在线高清观看一区二区三区| 欧美日韩精品成人综合77777| 亚洲高清免费不卡视频| 午夜视频国产福利| 欧美精品一区二区大全| 看免费av毛片| 极品少妇高潮喷水抽搐| 人成视频在线观看免费观看| 久久99热这里只频精品6学生| 亚洲图色成人| 制服丝袜香蕉在线| 久久久久人妻精品一区果冻| 高清毛片免费看| 日本黄大片高清| 久久午夜综合久久蜜桃| 亚洲少妇的诱惑av| 大香蕉久久成人网| videossex国产| 免费人成在线观看视频色| 久久99精品国语久久久| 日日摸夜夜添夜夜爱| 日本91视频免费播放| 亚洲激情五月婷婷啪啪| 久久国产亚洲av麻豆专区| 欧美精品高潮呻吟av久久| 亚洲精品乱码久久久久久按摩| 亚洲精品久久午夜乱码| 天天影视国产精品| 18在线观看网站| 久久 成人 亚洲| 日本欧美视频一区| 日产精品乱码卡一卡2卡三| 亚洲国产精品一区三区| 乱码一卡2卡4卡精品| 人成视频在线观看免费观看| 男女国产视频网站| 久久久久久久国产电影| 成人毛片a级毛片在线播放| 日韩制服骚丝袜av| 毛片一级片免费看久久久久| 韩国高清视频一区二区三区| 老女人水多毛片| 免费av不卡在线播放| 国产精品国产三级专区第一集| 久久久久视频综合| 91精品国产国语对白视频| 99热6这里只有精品| 久久国产精品男人的天堂亚洲 | 日韩av免费高清视频| 好男人视频免费观看在线| 午夜激情av网站| 在线观看一区二区三区激情| 国产一区二区激情短视频 | videossex国产| 九九爱精品视频在线观看| 亚洲久久久国产精品| 欧美日韩成人在线一区二区| 国产成人av激情在线播放| 在线观看人妻少妇| 中文字幕另类日韩欧美亚洲嫩草| 大码成人一级视频| 五月天丁香电影| 欧美人与善性xxx| 精品一区二区免费观看| 99热6这里只有精品| 日韩欧美一区视频在线观看| 欧美精品人与动牲交sv欧美| 久久青草综合色| 啦啦啦啦在线视频资源| 亚洲经典国产精华液单| 久久 成人 亚洲| 亚洲性久久影院| 日韩,欧美,国产一区二区三区| 全区人妻精品视频| 亚洲国产欧美在线一区| 日韩 亚洲 欧美在线| h视频一区二区三区| 欧美精品一区二区大全| 国产精品欧美亚洲77777| 一边摸一边做爽爽视频免费| 97超碰精品成人国产| 国产国语露脸激情在线看| 嫩草影院入口| 精品亚洲成a人片在线观看| 大陆偷拍与自拍| 欧美激情国产日韩精品一区| 女人精品久久久久毛片| 久久久久国产精品人妻一区二区| 日韩成人伦理影院| 国产国语露脸激情在线看| 最近手机中文字幕大全| 国产精品不卡视频一区二区| 大香蕉久久成人网| 亚洲欧美清纯卡通| 国产xxxxx性猛交| 亚洲色图 男人天堂 中文字幕 | 成人影院久久| 免费观看a级毛片全部| 成人国语在线视频| 久久精品夜色国产| 久久久国产精品麻豆| 视频区图区小说| 久久人人爽人人片av| 亚洲成人av在线免费| 99热网站在线观看| 51国产日韩欧美| 精品酒店卫生间| 99热国产这里只有精品6| 亚洲高清免费不卡视频| 免费少妇av软件| 久久久久久久久久久久大奶| 国产精品女同一区二区软件| 乱人伦中国视频| 亚洲欧美精品自产自拍| 国产亚洲最大av| 欧美 日韩 精品 国产| 黑丝袜美女国产一区| 国产成人精品婷婷| 欧美精品亚洲一区二区| 黄色配什么色好看| 亚洲综合色网址| 国产一级毛片在线| 亚洲av电影在线进入| 欧美少妇被猛烈插入视频| 熟妇人妻不卡中文字幕| 精品国产一区二区三区久久久樱花| 亚洲综合色惰| 人妻一区二区av| 十分钟在线观看高清视频www| 久久99热6这里只有精品| 久久国产精品大桥未久av| 男女高潮啪啪啪动态图| 97在线人人人人妻| 国产精品成人在线| 欧美 亚洲 国产 日韩一| 国产成人一区二区在线| 精品久久蜜臀av无| 国产乱来视频区| 免费大片黄手机在线观看| 亚洲成国产人片在线观看| 三级国产精品片| 少妇高潮的动态图| 日韩三级伦理在线观看| 欧美 日韩 精品 国产| 成年美女黄网站色视频大全免费| 午夜久久久在线观看| 久热久热在线精品观看| 国产熟女欧美一区二区| 国产片内射在线| 国产亚洲精品第一综合不卡 | 毛片一级片免费看久久久久| 久久久精品免费免费高清| 亚洲国产av影院在线观看| 97精品久久久久久久久久精品| 亚洲av电影在线进入| 国产日韩一区二区三区精品不卡| 精品人妻一区二区三区麻豆| 亚洲精品,欧美精品| 99久久精品国产国产毛片| 国产探花极品一区二区| 一级a做视频免费观看| 国产在线一区二区三区精| 欧美亚洲 丝袜 人妻 在线| 欧美国产精品va在线观看不卡| 日韩在线高清观看一区二区三区| 欧美亚洲 丝袜 人妻 在线| 亚洲av中文av极速乱| 欧美精品av麻豆av| 亚洲av在线观看美女高潮| 22中文网久久字幕| 一二三四中文在线观看免费高清| 两个人免费观看高清视频| 菩萨蛮人人尽说江南好唐韦庄| 性色av一级| 天天躁夜夜躁狠狠久久av| 成年人午夜在线观看视频| 人人妻人人澡人人爽人人夜夜| 青春草亚洲视频在线观看| 最近手机中文字幕大全| 少妇人妻精品综合一区二区| 国产av码专区亚洲av| 亚洲情色 制服丝袜| 亚洲国产看品久久| 菩萨蛮人人尽说江南好唐韦庄| 一边摸一边做爽爽视频免费| 国产熟女欧美一区二区| 欧美另类一区| 欧美xxxx性猛交bbbb| 亚洲精品美女久久久久99蜜臀 | 97人妻天天添夜夜摸| 伊人久久国产一区二区| 制服诱惑二区| 春色校园在线视频观看| 十八禁网站网址无遮挡| av免费观看日本| 下体分泌物呈黄色| 插逼视频在线观看| 亚洲色图综合在线观看| 欧美变态另类bdsm刘玥| av免费在线看不卡| 欧美亚洲日本最大视频资源| 日韩一本色道免费dvd| 一本大道久久a久久精品| 五月开心婷婷网| 黄片播放在线免费| 久久这里有精品视频免费| 中国三级夫妇交换| 欧美 日韩 精品 国产| 汤姆久久久久久久影院中文字幕| 一本大道久久a久久精品| 久久女婷五月综合色啪小说| 国产综合精华液| 男女无遮挡免费网站观看| 日韩精品有码人妻一区| 国产精品久久久久久精品古装| 最近中文字幕2019免费版| 成人无遮挡网站| 日本-黄色视频高清免费观看| 男人爽女人下面视频在线观看| 日韩av在线免费看完整版不卡| 欧美精品亚洲一区二区| 香蕉国产在线看| 久久久久久人人人人人| 我的女老师完整版在线观看| 国产av一区二区精品久久| 青青草视频在线视频观看| 国产免费现黄频在线看| 99久国产av精品国产电影| a级毛片在线看网站| 国产亚洲一区二区精品| 精品一区二区免费观看| 秋霞伦理黄片| 看免费av毛片| 全区人妻精品视频| 97精品久久久久久久久久精品| 日韩精品免费视频一区二区三区 | xxx大片免费视频| 成人午夜精彩视频在线观看| 人人澡人人妻人| 欧美日本中文国产一区发布| 日韩免费高清中文字幕av| 国产女主播在线喷水免费视频网站| 秋霞在线观看毛片| 男人操女人黄网站| 91成人精品电影| 男女边吃奶边做爰视频| 久久鲁丝午夜福利片| videossex国产| av网站免费在线观看视频| 少妇 在线观看| 亚洲精品乱码久久久久久按摩| 卡戴珊不雅视频在线播放| 最近手机中文字幕大全| 成人亚洲精品一区在线观看| 亚洲高清免费不卡视频| 美女脱内裤让男人舔精品视频| 国产在视频线精品| 国产免费一级a男人的天堂| 黄色视频在线播放观看不卡| 视频在线观看一区二区三区| 亚洲国产色片| 在线看a的网站| 午夜福利视频精品| 黄色配什么色好看| 亚洲情色 制服丝袜| 热99久久久久精品小说推荐| 久久人人爽人人爽人人片va| 97在线人人人人妻| 日韩精品有码人妻一区| 久热久热在线精品观看| 91aial.com中文字幕在线观看| 一二三四在线观看免费中文在 | 两性夫妻黄色片 | 18+在线观看网站| 久久久国产一区二区| 国产精品麻豆人妻色哟哟久久| 亚洲av成人精品一二三区| 免费看av在线观看网站| 亚洲一区二区三区欧美精品| 亚洲四区av| 中文字幕另类日韩欧美亚洲嫩草| 高清不卡的av网站| 91精品三级在线观看| 亚洲精品色激情综合| 夜夜爽夜夜爽视频| 中文字幕最新亚洲高清| 人人妻人人爽人人添夜夜欢视频| 亚洲情色 制服丝袜| 久久久久久久久久成人| 精品少妇久久久久久888优播| 欧美日本中文国产一区发布| 精品国产乱码久久久久久小说| 高清欧美精品videossex| 国产69精品久久久久777片| 国产黄频视频在线观看| 久久久久久久国产电影| 黄片无遮挡物在线观看| 黄色 视频免费看| 中文字幕人妻熟女乱码| 久久久久久人妻| 女人精品久久久久毛片| 91成人精品电影| 777米奇影视久久| 久久毛片免费看一区二区三区| 欧美少妇被猛烈插入视频| 国产午夜精品一二区理论片| 免费av中文字幕在线| 宅男免费午夜| 日韩中文字幕视频在线看片| 久久综合国产亚洲精品| 日韩精品免费视频一区二区三区 | 丰满乱子伦码专区| 日韩欧美一区视频在线观看| 国产视频首页在线观看| 赤兔流量卡办理| 国产亚洲午夜精品一区二区久久| 成年女人在线观看亚洲视频| 七月丁香在线播放| 晚上一个人看的免费电影| 国产精品无大码| 韩国av在线不卡| 边亲边吃奶的免费视频| 9热在线视频观看99| 全区人妻精品视频| 卡戴珊不雅视频在线播放| 国产免费现黄频在线看| 卡戴珊不雅视频在线播放| 亚洲经典国产精华液单| 欧美日韩国产mv在线观看视频| 亚洲人成77777在线视频| 咕卡用的链子| 99视频精品全部免费 在线| 最近手机中文字幕大全| 九九在线视频观看精品| 亚洲综合精品二区| 伦理电影大哥的女人| 亚洲 欧美一区二区三区| 亚洲人成77777在线视频| 国产精品嫩草影院av在线观看| 国产一区二区在线观看日韩| 亚洲丝袜综合中文字幕| 少妇猛男粗大的猛烈进出视频| 久久久亚洲精品成人影院| 亚洲伊人色综图| 成人18禁高潮啪啪吃奶动态图| 香蕉丝袜av| 国产无遮挡羞羞视频在线观看| 五月玫瑰六月丁香| 伦理电影大哥的女人| 99热6这里只有精品| 一本色道久久久久久精品综合| 亚洲精品美女久久av网站| 欧美精品亚洲一区二区| av片东京热男人的天堂| 亚洲精品美女久久久久99蜜臀 | 男人操女人黄网站| h视频一区二区三区| videossex国产| 两性夫妻黄色片 | 中国三级夫妇交换| 亚洲国产日韩一区二区| 少妇 在线观看| 爱豆传媒免费全集在线观看| 亚洲精品美女久久av网站| 狂野欧美激情性bbbbbb| 亚洲一级一片aⅴ在线观看| 久久鲁丝午夜福利片| 久久久久国产精品人妻一区二区| 国产麻豆69| 91精品三级在线观看| 欧美精品亚洲一区二区| 最近最新中文字幕免费大全7| 看免费成人av毛片| 99热全是精品| 中文精品一卡2卡3卡4更新| 国产成人精品在线电影| 热re99久久精品国产66热6| 国产男女内射视频| 夫妻性生交免费视频一级片| 热re99久久精品国产66热6| 午夜久久久在线观看| 在线天堂最新版资源| 午夜精品国产一区二区电影|