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

    船用跨臨界CO2引射制冷系統(tǒng)動(dòng)態(tài)性能研究

    2016-11-04 12:35:52周愛(ài)民朱求源張?jiān)缧?/span>
    艦船科學(xué)技術(shù) 2016年9期
    關(guān)鍵詞:引射器制冷系統(tǒng)氣液

    周愛(ài)民,朱求源,余 濤,段 晨,張?jiān)缧?/p>

    ( 1. 武漢第二船舶設(shè)計(jì)研究所,湖北 武漢 430064;2. 西安交通大學(xué) 化工學(xué)院,陜西 西安 710049)

    船用跨臨界CO2引射制冷系統(tǒng)動(dòng)態(tài)性能研究

    周愛(ài)民1,朱求源1,余濤1,段晨1,張?jiān)缧?

    ( 1. 武漢第二船舶設(shè)計(jì)研究所,湖北武漢 430064;2. 西安交通大學(xué) 化工學(xué)院,陜西西安 710049)

    針對(duì)船用跨臨界CO2兩相流引射制冷系統(tǒng)建立動(dòng)態(tài)模型,模型中蒸發(fā)器采用移動(dòng)邊界模型,氣冷器采用集總參數(shù)模型,引射器使用 1-D 恒壓混合模型。通過(guò) Matlab 計(jì)算動(dòng)態(tài)模型,分析系統(tǒng)各參數(shù)在壓縮機(jī)轉(zhuǎn)速、膨脹閥開(kāi)度及引射器尺寸階躍擾動(dòng)下的動(dòng)態(tài)響應(yīng)。模擬結(jié)果顯示系統(tǒng)對(duì)壓縮機(jī)轉(zhuǎn)速變化更加敏感,系統(tǒng)參數(shù)耦合作用強(qiáng),并且動(dòng)態(tài)響應(yīng)中存在逆向響應(yīng)。本文研究為控制器設(shè)計(jì)提供了依據(jù)和指導(dǎo)。

    跨臨界;引射制冷系統(tǒng);動(dòng)態(tài)性能;二氧化碳

    0 引 言

    制冷系統(tǒng)對(duì)溫室效應(yīng)的影響包括兩方面:制冷劑泄漏和系統(tǒng)能耗所產(chǎn)生的 CO2排放。為了降低制冷劑對(duì)溫室效應(yīng)的作用,Lorrentzen 提出使用自然工質(zhì) CO2作為替代制冷劑,并考慮 CO2高臨界壓力、低臨界溫度的特點(diǎn),設(shè)計(jì)了跨臨界 CO2制冷循環(huán)。自 20 世紀(jì)90 年代以來(lái),國(guó)內(nèi)科研院所對(duì) CO2在船用制冷系統(tǒng)中應(yīng)用的研究方興未艾,由于傳統(tǒng)的跨臨界 CO2制冷循環(huán)中,存在高低壓差大,節(jié)流過(guò)程能量損失嚴(yán)重,系統(tǒng)性能較差。為了降低節(jié)流損失,有學(xué)者提出跨臨界CO2引射制冷循環(huán),使用引射器代替節(jié)流閥回收部分高壓能,從而提高系統(tǒng)性能。目前已有很多學(xué)者就跨臨界 CO2引射制冷循環(huán)進(jìn)行理論和實(shí)驗(yàn)研究[1-2],研究顯示引射器的使用對(duì)系統(tǒng)性能有很大提高。此外,動(dòng)態(tài)性能是制冷系統(tǒng)研究的另一個(gè)重要方面,是系統(tǒng)控制系統(tǒng)設(shè)計(jì)的基礎(chǔ),然而目前鮮有跨臨界 CO2引射制冷系統(tǒng)動(dòng)態(tài)性能的研究成果發(fā)表。

    本文針對(duì)跨臨界 CO2兩相流引射制冷系統(tǒng),通過(guò)Matlab 軟件建立動(dòng)態(tài)模型,研究在各設(shè)備運(yùn)行參數(shù)階躍擾動(dòng)下系統(tǒng)的動(dòng)態(tài)性能,為系統(tǒng)控制策略設(shè)計(jì)提供依據(jù)。

    1 系統(tǒng)流程及動(dòng)態(tài)模型

    1.1跨臨界 CO2引射制冷循環(huán)流程

    典型跨臨界 CO2引射制冷循環(huán)流程如圖1所示,引射器的結(jié)構(gòu)如圖2所示。高壓超臨界 CO2在氣冷器中放熱降溫,作為主動(dòng)流在引射器噴嘴中產(chǎn)生膨脹,在噴嘴出口產(chǎn)生低壓超音速兩相流,在引射段引射來(lái)自蒸發(fā)器的過(guò)熱 CO2,兩流體在混合段產(chǎn)生混合動(dòng)量交換,壓力升高,通過(guò)擴(kuò)壓段再次提升壓力,在引射器出口獲得中間壓力的兩相 CO2;流體經(jīng)過(guò)氣液分離器分離,氣相 CO2通過(guò)壓縮機(jī)升壓升溫進(jìn)入氣冷器,完成高壓循環(huán);液相 CO2在膨脹閥中節(jié)流產(chǎn)生低壓流體進(jìn)入蒸發(fā)器完成低壓循環(huán)。因此,引射器通過(guò)利用高壓側(cè)流體的壓力引射低壓流體,獲得中間壓力,降低了節(jié)流損失和壓縮功,從而提高系統(tǒng)的性能。

    圖1 跨臨界 CO2引射制冷系統(tǒng)流程Fig. 1 Transcritical CO2two-phase ejector expansion refrigeration system drawing

    圖2 引射器結(jié)構(gòu)圖Fig. 2 Ejector structure drawing

    本文在建立各設(shè)備動(dòng)態(tài)模型的基礎(chǔ)上,通過(guò)分析設(shè)備之間的質(zhì)量、能量守恒關(guān)系,模擬計(jì)算系統(tǒng)的動(dòng)態(tài)性能,分析在壓縮機(jī)轉(zhuǎn)速、膨脹閥開(kāi)度及噴射器噴嘴尺寸階躍擾動(dòng)下,系統(tǒng)各參數(shù)的階躍響應(yīng)。

    1.2系統(tǒng)動(dòng)態(tài)模型

    在建立系統(tǒng)各設(shè)備模型時(shí),考慮蒸發(fā)器、氣冷器、氣液分離器為工質(zhì)存儲(chǔ)設(shè)備,慣性作用較大,時(shí)間常數(shù)大,有延遲效應(yīng),選擇動(dòng)態(tài)模型;而壓縮機(jī)、膨脹閥、引射器作為流量控制設(shè)備,時(shí)間常數(shù)小,為了降低計(jì)算的復(fù)雜度,選用穩(wěn)態(tài)模型,并不會(huì)對(duì)系統(tǒng)動(dòng)態(tài)性能有很大影響。

    蒸發(fā)器采用套管式換熱器,二級(jí)流體為水,并假設(shè)出口總處于過(guò)熱狀態(tài)。在蒸發(fā)器中,換熱過(guò)程可分為兩相相變換熱區(qū)和單相過(guò)熱換熱區(qū) 2個(gè)區(qū)。為了在模型中反映不同的換熱段,蒸發(fā)器模型選擇移動(dòng)邊界模型。該模型在制冷系統(tǒng)動(dòng)態(tài)分析中經(jīng)常使用,降低了模型計(jì)算復(fù)雜度,同時(shí)能夠反應(yīng)蒸發(fā)器中由相變換熱到單相換熱的特點(diǎn)。

    氣冷器中 CO2為超臨界態(tài),換熱過(guò)程不存在相變,采用套管式換熱器,二級(jí)流體為水。模型采用集總參數(shù)模型。該方法通過(guò)建立氣冷器質(zhì)量守恒、能量守恒的偏微分方程,沿管長(zhǎng)方向積分,從而獲得設(shè)備常微分方程組。具體的換熱器模型建立過(guò)程可參考文獻(xiàn)[3],二氧化碳的換熱系數(shù)計(jì)算采用文獻(xiàn)[4-5]的經(jīng)驗(yàn)公式計(jì)算。

    壓縮機(jī)穩(wěn)態(tài)模型基于流量和能量方程,設(shè)備容積效率和等熵效率使用文獻(xiàn)[6-7]所提供的經(jīng)驗(yàn)公式;膨脹閥中的膨脹過(guò)程假設(shè)為等焓節(jié)流過(guò)程,流量由孔板節(jié)流流量計(jì)算公式計(jì)算,具體公式參考文獻(xiàn)[8]。

    如圖2所示,引射器包括噴嘴、引射段、混合段、擴(kuò)壓段 4個(gè)部分。假設(shè)流體在入口、出口處速度為0;噴嘴出口壓力等于蒸發(fā)壓力,引射流體在引射段內(nèi)沒(méi)有膨脹過(guò)程;混合段中主動(dòng)流與引射流壓力相同,是恒壓混合過(guò)程。對(duì)每個(gè)部分建立質(zhì)量、能量及動(dòng)量守恒方程,所得模型如式(1)~式(7)所示:

    式中的等熵系數(shù)采用文獻(xiàn)[2]所提供的值,假設(shè) ηn為0.7,ηm為0.95,ηd為0.8。

    氣液分離器模型如圖3所示。應(yīng)用質(zhì)量、能量守恒方程,選擇壓力和制冷劑儲(chǔ)量為狀態(tài)變量,建立狀態(tài)空間形式的動(dòng)態(tài)模型如式(8)~式(10)。

    圖3 氣液分離器模型Fig. 3 Gas-Liquid separator model drawing

    式中系數(shù):

    2 仿真結(jié)果及分析

    2.1設(shè)備參數(shù)

    根據(jù)圖1中所示的系統(tǒng)流程,本文設(shè)計(jì)的設(shè)備主要參數(shù)如表1所示。

    表1 系統(tǒng)設(shè)備主要參數(shù)表Tab. 1 The system parts'main parameters

    2.2仿真計(jì)算工況

    在模擬計(jì)算時(shí),考慮初始穩(wěn)態(tài)工況如表2所示。將本文開(kāi)發(fā)的模型應(yīng)用于上述系統(tǒng)穩(wěn)態(tài)工況,在 Matlab 中模擬計(jì)算系統(tǒng)的動(dòng)態(tài)性能。穩(wěn)態(tài)工況下,2 s 時(shí)引射器噴嘴面積階躍降低 5%;在一定時(shí)間的穩(wěn)定后,在 100 s 時(shí)膨脹閥開(kāi)度階躍降低 5%;該階段穩(wěn)定后,在 200 s 時(shí)壓縮機(jī)轉(zhuǎn)速突然階躍增加 5%。根據(jù)上述一系列擾動(dòng)下計(jì)算主要系統(tǒng)參數(shù)動(dòng)態(tài)響應(yīng)。

    表2 仿真計(jì)算初始穩(wěn)態(tài)工況Tab. 2 The original steady condition by simulation

    2.3仿真結(jié)果分析

    如圖4~圖7和表3計(jì)算結(jié)果所示,在初始噴嘴面積減少時(shí),噴嘴中膨脹作用增強(qiáng),噴嘴出口壓力降低,引射作用增強(qiáng),使得蒸發(fā)器液體儲(chǔ)量減少,壓力降低;同時(shí)主動(dòng)流流量減少,氣冷器出口流量突然減低,導(dǎo)致壓力升高;引射器混合流體變化小,氣液分離器的儲(chǔ)液量變化較小。所以在噴嘴出口面積減少時(shí),更多的 CO2將向氣冷器中集中。

    圖4 蒸發(fā)壓力動(dòng)態(tài)響應(yīng)Fig. 4 The dynamic response of evaporator pressure

    圖5 氣冷器壓力的動(dòng)態(tài)響應(yīng)Fig. 5 The dynamic response of gas cooler pressure

    圖6 氣液分離器壓力的動(dòng)態(tài)響應(yīng)Fig. 6 The dynamic response of gas-liquid separator pressure

    當(dāng)膨脹閥開(kāi)度減少時(shí),進(jìn)入蒸發(fā)器的流量減少,蒸發(fā)壓力降低,導(dǎo)致引射流量減少;在較大的噴嘴膨脹壓差作用下,主動(dòng)流流量增加,氣冷器壓力降低。在膨脹閥開(kāi)度降低時(shí),CO2將向氣液分離器集中,氣液分離器中液位上升。

    圖7 儲(chǔ)液量的動(dòng)態(tài)響應(yīng)Fig. 7 The dynamic response of supplying quantity

    表3 相同時(shí)間下系統(tǒng)參數(shù)動(dòng)態(tài)變化Tab. 3 The system parameters'dynamic trends under same time rate of change

    當(dāng)壓縮機(jī)轉(zhuǎn)速增加時(shí),更多的氣相 CO2壓縮進(jìn)入氣冷器中,氣冷器壓力升高,主動(dòng)流引射作用增強(qiáng),引射蒸發(fā)器中的過(guò)熱 CO2,使得蒸發(fā)器壓力降低;氣液分離器中氣相減少,壓力降低,液位升高。另外,從圖中分析可知,在相同時(shí)間(100s)變化的擾動(dòng)下,系統(tǒng)參數(shù)對(duì)壓縮機(jī)轉(zhuǎn)速擾動(dòng)更敏感。

    結(jié)合本文將開(kāi)展的系統(tǒng)運(yùn)行參數(shù)控制方面的研究,從圖4~圖7也可以看出,不同的控制輸入對(duì)系統(tǒng)各參數(shù)都有很明顯的作用,系統(tǒng)輸入輸出通道之間有強(qiáng)烈耦合作用,多通道 SISO 控制器的設(shè)計(jì)存在困難,而且效果可能較差;圖中顯示在膨脹閥開(kāi)度階躍變化時(shí)會(huì)產(chǎn)生逆向響應(yīng),即系統(tǒng)模型存在 RHP 零點(diǎn),這一性質(zhì)將限制系統(tǒng)響應(yīng)速度,膨脹閥開(kāi)度擾動(dòng)下系統(tǒng)穩(wěn)定時(shí)間長(zhǎng)也驗(yàn)證了這一結(jié)論,同時(shí) RHP 零點(diǎn)也會(huì)增加控制器設(shè)計(jì)的復(fù)雜度。

    3 結(jié) 語(yǔ)

    本文通過(guò)建立的船用跨臨界 CO2兩相流引射制冷系統(tǒng)的動(dòng)態(tài)模型,在一定尺寸的設(shè)備及穩(wěn)態(tài)工況下,通過(guò) Matlab 求解模型,獲得了系統(tǒng)運(yùn)行參數(shù)(蒸發(fā)壓力、氣冷器壓力、氣液分離器壓力和儲(chǔ)液量)在壓縮機(jī)轉(zhuǎn)速、膨脹閥開(kāi)度、引射器噴嘴尺寸的階躍擾動(dòng)下的動(dòng)態(tài)響應(yīng)。分析結(jié)果顯示:

    1)系統(tǒng)運(yùn)行參數(shù),尤其是各運(yùn)行壓力對(duì)壓縮機(jī)轉(zhuǎn)速更加敏感;

    2)在膨脹閥作為擾動(dòng)輸入的通道中存在 RHP 零點(diǎn),表現(xiàn)在 100 s 膨脹閥開(kāi)度出現(xiàn)擾動(dòng)時(shí),系統(tǒng)運(yùn)行參數(shù),如蒸發(fā)壓力、氣冷器壓力和氣液分離器壓力均先短時(shí)上升后再快速下降,會(huì)產(chǎn)生逆向響應(yīng);

    3)系統(tǒng)各輸入對(duì)所有系統(tǒng)參數(shù)都有明顯作用,耦合作用強(qiáng)烈。這些動(dòng)態(tài)特性會(huì)增加制冷系統(tǒng)控制器設(shè)計(jì)的復(fù)雜度,而本文關(guān)于動(dòng)態(tài)特性的研究對(duì)系統(tǒng)控制器的設(shè)計(jì)有重要的指導(dǎo)意義。

    [1]LEE J S, KIM M S, KIM M S. Experimental study on the improvement of CO2air conditioning system performance using an ejector[J]. International Journal of Refrigeration, 2011, 34(7):1614-1625.

    [2]DENG J Q, JIANG P X, LU T, et al. Particular characteristics of transcritical CO2refrigeration cycle with an ejector[J]. Applied Thermal Engineering, 2007, 27(2/3): 381-388.

    [3]ALLEYNE A G, RASMUSSEN B P. Control-oriented modeling of transcritical vapor compression systems[J]. Journal of Dynamic Systems, Measurement and Control, 2004, 126(1):54-64.

    [4]YOON S K, KIM J H, HWANG Y W, et al. Heat transfer and pressure drop characteristics during the in-tube cooling process of carbon dioxide in the supercritical region[J]. International Journal of Refrigeration, 2003, 26(8): 857-864.

    [5]KIM M H, BULLARD C W. Development of a microchannel evaporator model for a CO2air-conditioning system[J]. Energy,2001, 26(10): 931-948.

    [6]劉洪勝, 楊濤, 陳江平, 等. 跨臨界二氧化碳制冷系統(tǒng)動(dòng)態(tài)性能仿真研究[J]. 工程熱物理學(xué)報(bào), 2006, 27(S2): 69-72. LIU Hong-sheng, YANG Tao, CHEN Jiang-ping, et al. Research on dynamic simulation of transcritical CO2refrigeration system[J]. Journal of Engineering Thermophysics, 2006, 27(S2):69-72.

    [7]楊濤, 陳江平, 陳芝久. 二氧化碳汽車(chē)空調(diào)系統(tǒng)的動(dòng)態(tài)仿真[J].上海汽車(chē), 2007(6): 15-19. YANG Tao, CHEN Jiang-ping, CHEN Zhi-jiu. The dynamic simulation of the air condition system of the CO2in vehicle[J]. Shanghai Auto, 2007(6): 15-19.

    [8]張?jiān)缧#?馮霄, 郁永章. 制冷與熱泵[M]. 北京: 化學(xué)工業(yè)出版社, 2000. ZHANG Zao-xiao, FENG Xiao, YU Yong-zhang. Refrigeration and Heat Pump[M]. Beijing: Chemical Industry Press, 2000.

    Research on dynamic simulation of marine transcritical CO2ejector expansion refrigeration system

    ZHOU Ai-min1, ZHU Qiu-yuan1, YU Tao1, DUAN Chen1, ZHANG Zao-xiao2
    (1. Wuhan Second Ship Design and Research Institute, Wuhan 430064, China 2. School of Chemical Engineering and Technology, Xi'an Jiaotong University, Xi'an 710049, China)

    This paper firstly establishes the dynamic simulation model for marine transcritical CO2two-phase ejector expansion refrigeration system with Matlab softwareincluding the evaporator using the moving boundary model, the gas cooler using the lumped parameter model and the ejector using the 1-D constant pressure hybrid model. The dynamic response is analyzed under the step disturbances of system parameters including the compressor's speed, the expansion valve's opening degree and the ejector's dimension to be turbulent in by means of Matlab's calculation dynamic model. The simulation results indicate that the compressor's speed status is more sensitive than other system parameters, and the couple action of system parameters are powerful. Also the system dynamic response exists the reverse case. This paper provides the reference for controller design.

    transcritical;ejector expansion refrigeration system;dynamic simulation;carbon dioxide

    TB66

    A

    1672-7619(2016)09-0126-04

    10.3404/j.issn.1672-7619.2016.09.026

    2016-03-02;

    2016-03-30

    周愛(ài)民(1979-),男,高級(jí)工程師,主要從事船舶大氣環(huán)境控制系統(tǒng)研究。

    猜你喜歡
    引射器制冷系統(tǒng)氣液
    Φ1.2 m 高超聲速風(fēng)洞引射系統(tǒng)設(shè)計(jì)與性能試驗(yàn)
    R290/R170單級(jí)壓縮回?zé)嶂评湎到y(tǒng)模擬研究
    微重力下兩相控溫型儲(chǔ)液器內(nèi)氣液界面仿真分析
    變循環(huán)發(fā)動(dòng)機(jī)后涵道引射器調(diào)節(jié)工況的數(shù)值研究
    氣液分離罐液位計(jì)接管泄漏分析
    高超聲速風(fēng)洞兩級(jí)引射器氣動(dòng)性能試驗(yàn)研究
    PEMFC系統(tǒng)引射器設(shè)計(jì)及仿真研究
    CO2 驅(qū)低液量高氣液比井下氣錨模擬與優(yōu)化
    奧迪A6L車(chē)制冷系統(tǒng)工作異常
    氣液接觸法制備球形氫氧化鎂
    国产欧美另类精品又又久久亚洲欧美| 日韩 亚洲 欧美在线| 偷拍熟女少妇极品色| 久久精品久久精品一区二区三区| 99久国产av精品国产电影| 三级国产精品片| 极品少妇高潮喷水抽搐| 亚洲最大成人中文| 国内精品美女久久久久久| 日韩欧美精品v在线| 最近最新中文字幕免费大全7| 久久精品综合一区二区三区| 亚洲最大成人手机在线| 成年av动漫网址| 波多野结衣巨乳人妻| 成人亚洲精品一区在线观看 | 国产精品麻豆人妻色哟哟久久| 人人妻人人看人人澡| 欧美最新免费一区二区三区| 蜜臀久久99精品久久宅男| 99热网站在线观看| 久久ye,这里只有精品| videossex国产| 午夜激情福利司机影院| 午夜精品一区二区三区免费看| 一个人看视频在线观看www免费| 我的老师免费观看完整版| 亚洲天堂av无毛| 亚洲国产高清在线一区二区三| freevideosex欧美| 亚洲欧洲日产国产| 一本久久精品| 日本三级黄在线观看| 午夜免费男女啪啪视频观看| av在线蜜桃| 国产一区二区在线观看日韩| 日韩强制内射视频| 热re99久久精品国产66热6| 免费av观看视频| 中文字幕av成人在线电影| 中文精品一卡2卡3卡4更新| 性色avwww在线观看| 久久久久久久久久久丰满| 白带黄色成豆腐渣| 久久久久久伊人网av| 大片免费播放器 马上看| 九九爱精品视频在线观看| 国产成人精品福利久久| 精品一区二区免费观看| 精品久久久精品久久久| 欧美极品一区二区三区四区| 小蜜桃在线观看免费完整版高清| 亚洲欧美精品自产自拍| 国产av码专区亚洲av| 2021天堂中文幕一二区在线观| 七月丁香在线播放| 欧美激情在线99| 51国产日韩欧美| 亚洲天堂国产精品一区在线| 久久久久国产网址| 哪个播放器可以免费观看大片| 国产成人福利小说| 最后的刺客免费高清国语| 又爽又黄无遮挡网站| 国产日韩欧美在线精品| 波野结衣二区三区在线| 久久精品久久久久久噜噜老黄| 亚洲av免费在线观看| 91在线精品国自产拍蜜月| 高清毛片免费看| 老师上课跳d突然被开到最大视频| 国产伦在线观看视频一区| 肉色欧美久久久久久久蜜桃 | 欧美日韩国产mv在线观看视频 | 免费看日本二区| 欧美97在线视频| 国产白丝娇喘喷水9色精品| 色综合色国产| 亚洲av中文av极速乱| 我的女老师完整版在线观看| 99九九线精品视频在线观看视频| 狂野欧美白嫩少妇大欣赏| 国产在线一区二区三区精| 麻豆成人av视频| 日韩人妻高清精品专区| 国产成人福利小说| 国产人妻一区二区三区在| 极品教师在线视频| 国产精品久久久久久精品古装| 在线观看av片永久免费下载| 大香蕉97超碰在线| 一级二级三级毛片免费看| xxx大片免费视频| 久久热精品热| 国产成人freesex在线| 夜夜爽夜夜爽视频| 久久精品国产a三级三级三级| 国产亚洲一区二区精品| 在线免费观看不下载黄p国产| 爱豆传媒免费全集在线观看| 欧美精品国产亚洲| 另类亚洲欧美激情| 麻豆精品久久久久久蜜桃| 亚洲精品亚洲一区二区| 2021天堂中文幕一二区在线观| 国内揄拍国产精品人妻在线| 亚洲高清免费不卡视频| 国产亚洲一区二区精品| 欧美激情国产日韩精品一区| 国产精品国产三级专区第一集| 亚洲av福利一区| 亚洲国产最新在线播放| tube8黄色片| 日本熟妇午夜| 久久女婷五月综合色啪小说 | 亚洲av欧美aⅴ国产| 国产又色又爽无遮挡免| tube8黄色片| 在线观看一区二区三区| 99视频精品全部免费 在线| 男女无遮挡免费网站观看| 亚洲精品456在线播放app| 最近2019中文字幕mv第一页| 97在线视频观看| 亚洲综合色惰| 我的女老师完整版在线观看| 久久99精品国语久久久| 国产黄频视频在线观看| 国产成人精品久久久久久| 最近中文字幕2019免费版| 国产精品.久久久| 精品一区二区三卡| 夜夜看夜夜爽夜夜摸| 哪个播放器可以免费观看大片| 精品久久国产蜜桃| 男人和女人高潮做爰伦理| 久久精品国产亚洲av天美| 亚洲欧美日韩东京热| 亚洲精品,欧美精品| 国产亚洲最大av| 国产精品熟女久久久久浪| 观看免费一级毛片| 国产有黄有色有爽视频| 少妇人妻 视频| 欧美激情在线99| 国产精品久久久久久精品古装| 91精品一卡2卡3卡4卡| 久久热精品热| 国产精品爽爽va在线观看网站| 最近2019中文字幕mv第一页| 日日摸夜夜添夜夜添av毛片| 国产一区二区三区av在线| 日韩大片免费观看网站| 日本一二三区视频观看| 天堂网av新在线| 精品国产一区二区三区久久久樱花 | 亚洲精品一区蜜桃| 十八禁网站网址无遮挡 | 女的被弄到高潮叫床怎么办| 欧美人与善性xxx| 国产成人精品久久久久久| 少妇被粗大猛烈的视频| av在线播放精品| 女人被狂操c到高潮| 直男gayav资源| 成年版毛片免费区| 婷婷色麻豆天堂久久| 精品一区在线观看国产| 韩国高清视频一区二区三区| 人妻制服诱惑在线中文字幕| 精品国产乱码久久久久久小说| av国产精品久久久久影院| 如何舔出高潮| 久久鲁丝午夜福利片| 亚洲精品成人久久久久久| 亚洲精品久久午夜乱码| 免费黄网站久久成人精品| 国产伦理片在线播放av一区| 26uuu在线亚洲综合色| 久久6这里有精品| 欧美极品一区二区三区四区| 成人亚洲精品一区在线观看 | 日韩一本色道免费dvd| 一级a做视频免费观看| 男人爽女人下面视频在线观看| 亚洲av成人精品一区久久| 精品人妻偷拍中文字幕| 亚洲婷婷狠狠爱综合网| 三级经典国产精品| 国产一区二区在线观看日韩| 亚洲四区av| 日韩成人伦理影院| 精品午夜福利在线看| 国产视频内射| 肉色欧美久久久久久久蜜桃 | 欧美+日韩+精品| 少妇人妻久久综合中文| 韩国高清视频一区二区三区| 九九在线视频观看精品| 黄色怎么调成土黄色| 久久久久久久亚洲中文字幕| 国产黄频视频在线观看| 久久国内精品自在自线图片| 女人久久www免费人成看片| 五月玫瑰六月丁香| 内地一区二区视频在线| 男女边摸边吃奶| 久久韩国三级中文字幕| 国产成人精品一,二区| 视频区图区小说| 国产精品熟女久久久久浪| 亚洲在久久综合| 亚洲一级一片aⅴ在线观看| 亚洲av免费高清在线观看| 亚洲在线观看片| 国产永久视频网站| 另类亚洲欧美激情| 午夜激情福利司机影院| 精品亚洲乱码少妇综合久久| 亚洲国产日韩一区二区| 国产av码专区亚洲av| 五月玫瑰六月丁香| 简卡轻食公司| 亚洲图色成人| 天堂网av新在线| 久久精品国产亚洲av涩爱| 久久热精品热| 日韩三级伦理在线观看| 国产精品爽爽va在线观看网站| 黄色日韩在线| 国产成年人精品一区二区| 日本av手机在线免费观看| 综合色丁香网| 久久精品国产亚洲网站| 又黄又爽又刺激的免费视频.| 一个人看视频在线观看www免费| 国产免费福利视频在线观看| 26uuu在线亚洲综合色| 亚洲伊人久久精品综合| 视频中文字幕在线观看| 秋霞在线观看毛片| .国产精品久久| 18禁在线播放成人免费| 久久久精品免费免费高清| 亚洲美女视频黄频| 99久国产av精品国产电影| 中文资源天堂在线| 最新中文字幕久久久久| 日本猛色少妇xxxxx猛交久久| 日本免费在线观看一区| 丰满人妻一区二区三区视频av| 成人特级av手机在线观看| 欧美成人精品欧美一级黄| 波野结衣二区三区在线| 26uuu在线亚洲综合色| 性色avwww在线观看| 亚洲欧美成人精品一区二区| 国产午夜精品久久久久久一区二区三区| 一级毛片久久久久久久久女| 欧美成人精品欧美一级黄| 久久鲁丝午夜福利片| 国产精品秋霞免费鲁丝片| 久久久久久国产a免费观看| 国产精品久久久久久精品电影小说 | 亚洲精品国产色婷婷电影| 亚洲四区av| 女人被狂操c到高潮| 赤兔流量卡办理| 国产探花极品一区二区| 美女主播在线视频| 91久久精品国产一区二区三区| 成年人午夜在线观看视频| 男女无遮挡免费网站观看| 欧美日韩一区二区视频在线观看视频在线 | 亚洲av.av天堂| 天天躁日日操中文字幕| 午夜精品国产一区二区电影 | 久久精品综合一区二区三区| xxx大片免费视频| 国产乱人视频| 日韩av不卡免费在线播放| 97超视频在线观看视频| 久久国内精品自在自线图片| 精品久久久噜噜| 日韩,欧美,国产一区二区三区| 国产精品成人在线| 九色成人免费人妻av| 最近手机中文字幕大全| 亚洲成人久久爱视频| 国产淫片久久久久久久久| 成人午夜精彩视频在线观看| 麻豆久久精品国产亚洲av| 国内揄拍国产精品人妻在线| 青春草亚洲视频在线观看| 51国产日韩欧美| 欧美一级a爱片免费观看看| 午夜福利在线观看免费完整高清在| 日韩大片免费观看网站| av卡一久久| 丝袜喷水一区| 中文资源天堂在线| 欧美日韩在线观看h| 少妇猛男粗大的猛烈进出视频 | 精品人妻偷拍中文字幕| 天天躁夜夜躁狠狠久久av| 日韩av不卡免费在线播放| 国产精品国产三级国产av玫瑰| 亚洲最大成人av| 亚洲欧美日韩另类电影网站 | 亚洲综合精品二区| 内地一区二区视频在线| 国产精品99久久久久久久久| 看非洲黑人一级黄片| 日韩制服骚丝袜av| 日韩欧美精品免费久久| 中文字幕人妻熟人妻熟丝袜美| 中文字幕久久专区| 日韩人妻高清精品专区| 亚州av有码| 舔av片在线| 免费黄色在线免费观看| 国产精品99久久久久久久久| freevideosex欧美| 欧美国产精品一级二级三级 | 少妇丰满av| 午夜福利视频精品| 真实男女啪啪啪动态图| 人妻制服诱惑在线中文字幕| 欧美日韩在线观看h| 一区二区三区四区激情视频| 久久精品久久久久久久性| 十八禁网站网址无遮挡 | 欧美人与善性xxx| eeuss影院久久| 免费观看性生交大片5| 国产精品蜜桃在线观看| 久久鲁丝午夜福利片| 亚洲av国产av综合av卡| 在线观看国产h片| 亚洲欧美一区二区三区国产| 日韩视频在线欧美| 精品国产露脸久久av麻豆| 成人特级av手机在线观看| 又爽又黄无遮挡网站| 99久久九九国产精品国产免费| 日韩一本色道免费dvd| 欧美成人a在线观看| 亚洲天堂国产精品一区在线| 老司机影院毛片| 国产黄a三级三级三级人| 久久久久性生活片| 免费av观看视频| 搡女人真爽免费视频火全软件| 日本欧美国产在线视频| 美女被艹到高潮喷水动态| 日本与韩国留学比较| 国产成人a∨麻豆精品| 又黄又爽又刺激的免费视频.| 亚洲精品日韩av片在线观看| 午夜爱爱视频在线播放| 免费看日本二区| 国产人妻一区二区三区在| 中文资源天堂在线| 亚洲精品自拍成人| 自拍欧美九色日韩亚洲蝌蚪91 | 日韩大片免费观看网站| 国产黄色视频一区二区在线观看| 欧美潮喷喷水| 婷婷色综合www| 亚洲精品国产av蜜桃| 人人妻人人爽人人添夜夜欢视频 | 亚洲精品色激情综合| 国产免费视频播放在线视频| 国产成人a区在线观看| 永久免费av网站大全| 亚洲av一区综合| 国产爽快片一区二区三区| 中文天堂在线官网| 日本av手机在线免费观看| 尤物成人国产欧美一区二区三区| 肉色欧美久久久久久久蜜桃 | 中文资源天堂在线| 国产成人免费观看mmmm| 狂野欧美激情性bbbbbb| 久久久久国产网址| 欧美少妇被猛烈插入视频| 97在线人人人人妻| 精品国产乱码久久久久久小说| 少妇 在线观看| 一个人看的www免费观看视频| 91在线精品国自产拍蜜月| 国产成人freesex在线| 观看美女的网站| 看黄色毛片网站| 日韩人妻高清精品专区| 黄色日韩在线| 女人久久www免费人成看片| 男女边摸边吃奶| 久久午夜福利片| 免费av观看视频| 亚洲美女视频黄频| 成年人午夜在线观看视频| 99久久人妻综合| 亚洲精品日韩在线中文字幕| 精品一区二区三区视频在线| 日韩国内少妇激情av| 日产精品乱码卡一卡2卡三| 亚洲欧洲日产国产| 又爽又黄无遮挡网站| 91久久精品电影网| 一本久久精品| 国产一区二区三区综合在线观看 | h日本视频在线播放| 国产国拍精品亚洲av在线观看| 亚洲伊人久久精品综合| 九色成人免费人妻av| 色婷婷久久久亚洲欧美| 国产又色又爽无遮挡免| 看非洲黑人一级黄片| 久久午夜福利片| 色5月婷婷丁香| 在线观看三级黄色| 嘟嘟电影网在线观看| 成人毛片60女人毛片免费| 国产黄a三级三级三级人| 观看免费一级毛片| 黄片无遮挡物在线观看| 久久久色成人| 尾随美女入室| 亚州av有码| 欧美成人午夜免费资源| 国产亚洲5aaaaa淫片| 春色校园在线视频观看| 搞女人的毛片| 一边亲一边摸免费视频| 蜜桃亚洲精品一区二区三区| 成人毛片a级毛片在线播放| 男女国产视频网站| 免费黄频网站在线观看国产| 亚洲精品日韩av片在线观看| 少妇人妻精品综合一区二区| 99精国产麻豆久久婷婷| 亚洲欧美成人综合另类久久久| 久久久久国产精品人妻一区二区| 免费观看av网站的网址| 女人十人毛片免费观看3o分钟| 亚洲av中文字字幕乱码综合| 天美传媒精品一区二区| 波野结衣二区三区在线| 麻豆成人午夜福利视频| 亚洲精品,欧美精品| 亚洲精品456在线播放app| 国产精品秋霞免费鲁丝片| 一级爰片在线观看| 国产精品精品国产色婷婷| 99久久九九国产精品国产免费| 亚洲高清免费不卡视频| 国产有黄有色有爽视频| 97精品久久久久久久久久精品| 国产 一区精品| 内地一区二区视频在线| 草草在线视频免费看| 国产亚洲av片在线观看秒播厂| 高清av免费在线| 亚洲无线观看免费| 97热精品久久久久久| 午夜福利在线观看免费完整高清在| 国产亚洲精品久久久com| 全区人妻精品视频| 好男人视频免费观看在线| 亚洲国产日韩一区二区| 观看美女的网站| 人人妻人人看人人澡| 久久久精品94久久精品| 中文字幕人妻熟人妻熟丝袜美| 91aial.com中文字幕在线观看| 久久久久国产精品人妻一区二区| 欧美xxⅹ黑人| 久久6这里有精品| 亚洲精品一二三| 高清日韩中文字幕在线| 人妻系列 视频| 人妻制服诱惑在线中文字幕| 免费大片18禁| 联通29元200g的流量卡| 少妇人妻 视频| 欧美高清性xxxxhd video| 中文欧美无线码| 免费观看的影片在线观看| .国产精品久久| 国产综合精华液| 中国美白少妇内射xxxbb| 久久影院123| 看黄色毛片网站| 亚洲精品国产成人久久av| 精品视频人人做人人爽| 中国三级夫妇交换| 直男gayav资源| 久久久久精品久久久久真实原创| 欧美 日韩 精品 国产| 婷婷色麻豆天堂久久| 欧美xxxx性猛交bbbb| 91久久精品电影网| 熟女人妻精品中文字幕| 男的添女的下面高潮视频| 丝袜脚勾引网站| 男的添女的下面高潮视频| 在线免费十八禁| 大话2 男鬼变身卡| 久久久久网色| 国产精品一区二区性色av| 男人爽女人下面视频在线观看| 1000部很黄的大片| 国产人妻一区二区三区在| 日韩欧美 国产精品| 人妻少妇偷人精品九色| 精品人妻偷拍中文字幕| 国产美女午夜福利| 日韩人妻高清精品专区| 亚洲精品第二区| 男人舔奶头视频| 午夜激情福利司机影院| 日产精品乱码卡一卡2卡三| a级毛色黄片| 亚洲精品日韩av片在线观看| 午夜精品国产一区二区电影 | 舔av片在线| 亚洲在久久综合| 久久97久久精品| 两个人的视频大全免费| 国产一级毛片在线| 国产成人精品婷婷| 九草在线视频观看| 黄色视频在线播放观看不卡| 亚洲av成人精品一二三区| 成人午夜精彩视频在线观看| 亚洲精品影视一区二区三区av| 成人毛片60女人毛片免费| 欧美日韩亚洲高清精品| 大片电影免费在线观看免费| 十八禁网站网址无遮挡 | 国模一区二区三区四区视频| 亚洲精品久久久久久婷婷小说| 亚洲美女搞黄在线观看| 夫妻午夜视频| 国产白丝娇喘喷水9色精品| 麻豆成人午夜福利视频| 免费播放大片免费观看视频在线观看| 成人午夜精彩视频在线观看| 欧美xxxx黑人xx丫x性爽| 亚洲精品久久久久久婷婷小说| av黄色大香蕉| 嫩草影院精品99| 国产有黄有色有爽视频| 久久久久网色| 汤姆久久久久久久影院中文字幕| 国产69精品久久久久777片| 汤姆久久久久久久影院中文字幕| 少妇熟女欧美另类| 一区二区三区精品91| 干丝袜人妻中文字幕| 国产一区二区三区av在线| 亚洲欧美成人综合另类久久久| 欧美成人午夜免费资源| 一级毛片 在线播放| 老师上课跳d突然被开到最大视频| 亚洲av福利一区| 国产在线男女| 一级二级三级毛片免费看| 日日摸夜夜添夜夜添av毛片| 18禁在线播放成人免费| 国产精品99久久久久久久久| 偷拍熟女少妇极品色| 日韩,欧美,国产一区二区三区| 毛片一级片免费看久久久久| 色综合色国产| 在线观看免费高清a一片| 熟女电影av网| 午夜爱爱视频在线播放| 2021天堂中文幕一二区在线观| 中文字幕亚洲精品专区| 一级毛片 在线播放| 可以在线观看毛片的网站| 欧美一区二区亚洲| 人妻夜夜爽99麻豆av| 内地一区二区视频在线| 中国国产av一级| 国内精品美女久久久久久| 蜜桃久久精品国产亚洲av| 成人特级av手机在线观看| 国产精品久久久久久精品电影| 99久久人妻综合| 亚州av有码| 波多野结衣巨乳人妻| 在线观看美女被高潮喷水网站| 大香蕉久久网| 免费观看的影片在线观看| 亚洲激情五月婷婷啪啪| 亚洲精品456在线播放app| www.色视频.com| 久久这里有精品视频免费| 午夜亚洲福利在线播放| 精品国产一区二区三区久久久樱花 | 视频区图区小说| 日韩成人伦理影院| 国产淫片久久久久久久久| www.色视频.com| 伊人久久精品亚洲午夜| 狂野欧美白嫩少妇大欣赏| 亚洲av免费在线观看| 成年版毛片免费区| 狂野欧美激情性bbbbbb| 亚洲欧美成人综合另类久久久| 免费少妇av软件| 成人漫画全彩无遮挡|