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

    多堆燃料電池系統(tǒng)溫度模型預(yù)測控制

    2022-10-08 09:50:38陳春光裴馮來
    關(guān)鍵詞:效果模型系統(tǒng)

    沈 偉,石 霖,陳春光,周 蘇,2,王 寧,裴馮來

    (1. 同濟(jì)大學(xué)汽車學(xué)院,上海 201804;2. 同濟(jì)大學(xué)中德學(xué)院,上海 201804;3. 上海機(jī)動(dòng)車檢測認(rèn)證技術(shù)研究中心有限公司,上海 201805)

    質(zhì)子交換膜燃料電池(proton exchange membrane fuel cell,PEMFC)以氫氣為燃料,除了對(duì)環(huán)境友好之外,還具有工作溫度低、比能量高、啟動(dòng)速度快和壽命較長等優(yōu)點(diǎn),是包括汽車動(dòng)力電源在內(nèi)的多個(gè)應(yīng)用領(lǐng)域的潛在零排放發(fā)電裝置之一[1-2]。但是,傳統(tǒng)的單堆PEMFC 系統(tǒng)很難滿足諸如商用車、船舶和重型機(jī)械等高功率需求應(yīng)用場景;并且單堆PEMFC 是由多個(gè)單池串聯(lián)組成的電結(jié)構(gòu),若某個(gè)單池出現(xiàn)故障,會(huì)導(dǎo)致整個(gè)單堆系統(tǒng)失效,這種低容錯(cuò)率也制約了PEMFC系統(tǒng)向大功率應(yīng)用場景的拓展[3-4]。為解決上述問題,可將多個(gè)單堆PEMFC進(jìn)行組合構(gòu)成多堆燃料電池系統(tǒng)(multi-stack fuel cell system,MFCS)。這樣不僅可以提高系統(tǒng)的輸出功率和能量轉(zhuǎn)換效率、拓展應(yīng)用場景,還可以增加系統(tǒng)運(yùn)行的冗余度,提高系統(tǒng)工作的可靠性。

    當(dāng)前,對(duì)MFCS 的研究主要集中于系統(tǒng)層面。Palma 等[5]利用DC/DC 將單個(gè)PEMFC 電堆分成多個(gè)模塊,提高了系統(tǒng)的容錯(cuò)率。Marx 等[6]對(duì)MFCS結(jié)構(gòu)的研究結(jié)果表明,MFCS 的輸出性能和壽命等與MFCS 結(jié)構(gòu)具有一定的關(guān)聯(lián)。Long 等[7]研究發(fā)現(xiàn),MFCS 并聯(lián)結(jié)構(gòu)有利于溫度控制。戴麗君[8]和馬天才等[9]針對(duì)MFCS熱管理問題設(shè)計(jì)了系統(tǒng)的管理方案以及零部件應(yīng)用方案。不論是單堆PEMFC系統(tǒng)還是MFCS,其涉及的電化學(xué)反應(yīng)區(qū)域都需要具備合適的溫度條件[10-11],這對(duì)燃料電池的催化劑活性、質(zhì)子交換膜含水量、電化學(xué)反應(yīng)速率、反應(yīng)氣體壓力和電堆壽命等都有重要影響[12-13]。因此,采用合適的熱管理子系統(tǒng)以及相應(yīng)的控制策略來維持MFCS各電堆的工作溫度是必不可少的。

    本文針對(duì)MFCS 各電堆溫度控制問題,借鑒單堆PEMFC 系統(tǒng)模型和熱管理控制方法,設(shè)計(jì)了面向MFCS 的并聯(lián)式熱管理子系統(tǒng)并建立了相關(guān)模型;基于熱平衡動(dòng)態(tài)分析,采用模型預(yù)測控制(model predictive control,MPC)算法對(duì)時(shí)滯性較大的各電堆溫度過程進(jìn)行實(shí)時(shí)控制;引入電堆氣體出口溫度作為修正項(xiàng)重新定義了電堆溫度;應(yīng)用模型參數(shù)辨識(shí)方法,建立對(duì)應(yīng)于典型工況點(diǎn)的多個(gè)預(yù)測模型,在不同工作點(diǎn)處通過模型切換,控制MFCS各電堆溫度;基于設(shè)計(jì)的測試工況,比較并分析了多個(gè)預(yù)測模型對(duì)各電堆溫度的控制效果。

    1 MFCS熱管理子系統(tǒng)模型搭建

    1.1 熱管理子系統(tǒng)

    適用于大功率應(yīng)用場景的MFCS 采用水冷方式,由于電堆數(shù)量的增加,需要重新設(shè)計(jì)熱管理子系統(tǒng)。如圖1所示,針對(duì)MFCS設(shè)計(jì)了并聯(lián)式熱管理子系統(tǒng)。旁通閥和分流器的引入可以調(diào)節(jié)各支路的冷卻液流量以滿足各電堆的冷卻需求,混合器用于混合各電堆出口冷卻液,節(jié)溫器用于控制大小循環(huán)。

    圖1 MFCS熱管理子系統(tǒng)Fig.1 MFCS thermal management subsystem

    1.2 熱平衡分析

    MFCS 的傳熱過程模型有如下假設(shè):①燃料電池外部環(huán)境溫度恒定;②氣體和冷卻液在管路內(nèi)的流動(dòng)為層流;③冷卻液與電堆充分均勻換熱;管路為絕熱系統(tǒng),忽略冷卻液在管路中流動(dòng)的散熱;④冷卻液在流經(jīng)中冷器組時(shí)與流進(jìn)電堆的空氣進(jìn)行換熱,換熱后的冷卻液溫度變化由查表所得,后續(xù)建模不再給出相關(guān)計(jì)算公式;⑤忽略水泵運(yùn)轉(zhuǎn)對(duì)系統(tǒng)冷卻液溫度的影響。

    涉及的MFCS 由3 個(gè)70kW 的PEMFC 單堆組成。假設(shè)3 個(gè)單堆的冷卻液與外界的換熱過程一致,故以單堆為例對(duì)其換熱過程進(jìn)行分析。單堆的熱平衡關(guān)系包括電堆與外界的熱交換和自身內(nèi)部電化學(xué)反應(yīng)產(chǎn)熱[12],可由式(1)表示,式(1)中的各項(xiàng)分別由式(2)-(5)表示:

    式中:Q·為熱功率值;I為電堆電流;n為電堆所包含的單池片數(shù);EHHV為電堆等效電壓;V為電堆對(duì)外輸出電壓;ε為斯特藩-玻爾茲曼常量;φ為電堆黑度;A為電堆表面積;Tst為電堆溫度;T0為環(huán)境溫度;h-為電堆表面對(duì)流換熱系數(shù);下標(biāo)an為陽極,ca為陰極,cool為冷卻液,in、out代表流入、流出電堆的物質(zhì),st為電堆,cell 為單池,react 為電化學(xué)反應(yīng),rad 為熱輻射,conv為熱對(duì)流。

    式(2)中各項(xiàng)包含了PEMFC 陰陽極氣體和電堆冷卻液流入流出所發(fā)生的熱交換[12],分別由式(6)-(11)表示:

    式中:m·為質(zhì)量流量;Cp為定壓比熱容;r為汽化潛熱;上標(biāo)vap為水蒸氣,liq為液態(tài)水。

    通過式(12)、(13)可以計(jì)算陰極氧氣和陽極氫氣的消耗量,然后通過維持陰陽極相應(yīng)的氣體過量系數(shù)可以實(shí)時(shí)計(jì)算得到陰極氣體的入堆流量。

    式中:M為分子質(zhì)量;F為法拉第常數(shù)。

    熱管理子系統(tǒng)的熱平衡過程可由圖2表示。散熱器出口為低溫冷卻液,混合器中為高溫冷卻液,二者之間的溫差是由散熱器散熱造成的。

    圖2 MFCS熱管理子系統(tǒng)熱平衡圖Fig.2 Heat balance diagram of MFCS thermal management subsystem

    混合器將流經(jīng)中冷器組和不同電堆的高溫冷卻液混合均勻,其熱平衡過程可由式(14)描述:

    散熱器將高溫冷卻液中的熱量釋放,其熱平衡過程可由式(15)表示:

    模型中各參數(shù)如表1所示。

    表1 模型參數(shù)Tab.1 Model parameter

    1.3 被控溫度定義

    在常規(guī)模型中,通常將PEMFC 冷卻液出口溫度視為電堆溫度,并且將此溫度作為被控溫度[14]。實(shí)際的PEMFC 內(nèi)部溫度分布并不均勻,熱管理本質(zhì)上要控制的是各單池電化學(xué)反應(yīng)區(qū)域(膜電極)的溫度,并且相關(guān)研究[15]表明,PEMFC 氣體出口溫度也與膜電極溫度相關(guān)。為了更好地表征電池內(nèi)部膜電極的溫度狀況,將氣體出口溫度作為修正項(xiàng),重新定義電堆溫度為

    式中:λ為相關(guān)物質(zhì)的溫度系數(shù),可以根據(jù)相關(guān)流體的比熱容和質(zhì)量流量確定。本文中,冷卻液、陰極和陽極氣體的溫度系數(shù)分別取0.894、0.100和0.006。

    2 MFCS溫度模型預(yù)測控制

    2.1 MPC控制器設(shè)計(jì)

    MPC算法的原理如圖3所示。MPC控制器利用每一個(gè)當(dāng)前采樣時(shí)刻獲取的測量信息,在線求解一個(gè)有限時(shí)域開環(huán)優(yōu)化問題,并將所得到的控制序列的第1個(gè)元素作為控制輸出,作用于被控對(duì)象。通過不斷重復(fù)該過程,得到每一時(shí)刻的最優(yōu)控制量[16]。

    圖3 MPC基本原理Fig.3 MPC fundamentals

    MFCS的熱管理子系統(tǒng)是一個(gè)內(nèi)部結(jié)構(gòu)復(fù)雜的多輸入多輸出系統(tǒng),其中各電堆的溫度過程時(shí)滯性較大。相較于其他算法而言,MPC算法能夠更好地處理多變量、大時(shí)滯的控制問題。因此,選用MPC算法對(duì)被控溫度進(jìn)行控制,以期獲得更好的控制效果。選取風(fēng)扇、水泵、分流器1 和2 作為執(zhí)行器??刂朴成潢P(guān)系由式(17)表示:

    控制器的輸出,即控制變量為風(fēng)扇轉(zhuǎn)速nfan、水泵轉(zhuǎn)速npump和分流器1、2的開度θdis1、θdis2,系統(tǒng)輸出為各電堆的溫度Tst1、Tst2、Tst3以及入堆冷卻液溫度Tst,in。MFCS各電堆負(fù)載電流Ii為可測的系統(tǒng)擾動(dòng)輸入。模型選擇器根據(jù)負(fù)載電流的大小選擇預(yù)測模型對(duì)應(yīng)的參數(shù)集。建立的MPC控制器結(jié)構(gòu)如圖4所示。

    圖4 MPC控制器結(jié)構(gòu)Fig.4 MPC controller architecture

    2.2 預(yù)測模型

    PEMFC電堆涉及的溫度過程非線性強(qiáng)、時(shí)滯性大,單一的預(yù)測模型無法反映不同工作條件下的溫度過程特性。因此,分別選取活化極化區(qū)后段、歐姆極化區(qū)全段和濃差極化區(qū)前段的典型工作點(diǎn),利用系統(tǒng)辨識(shí)方法建立相應(yīng)工作點(diǎn)的預(yù)測模型集。溫度MPC運(yùn)行時(shí),根據(jù)實(shí)際系統(tǒng)的工作狀態(tài)選用與之相匹配的預(yù)測模型。

    利用8 階M 序列生成辨識(shí)信號(hào),分別在單堆負(fù)載電流為80A、170A和250A這3個(gè)工作點(diǎn)處進(jìn)行系統(tǒng)辨識(shí)。以風(fēng)扇轉(zhuǎn)速、水泵轉(zhuǎn)速、分流器1 與2 的開度和3 個(gè)電堆的實(shí)際電流為被控系統(tǒng)的7 個(gè)輸入變量,各電堆溫度和入堆冷卻液溫度作為4 個(gè)輸出變量,這樣一個(gè)多輸入、多輸出熱過程的狀態(tài)空間模型可由式(18)表示,其中A、B和C分別代表系統(tǒng)矩陣、輸入矩陣和觀測矩陣。

    其中,狀態(tài)空間的維數(shù)(或模型的階數(shù))與MFCS熱管理子系統(tǒng)(被控對(duì)象)的熱儲(chǔ)能單元數(shù)量有關(guān)。如圖1所示,MFCS熱管理子系統(tǒng)中每個(gè)電堆具有電堆質(zhì)量體、冷卻液流道、空氣流道和氫氣流道4個(gè)熱儲(chǔ)能單元,堆外具有混合器、水箱、中冷器空氣流道與冷卻液流道和散熱器空氣流道與冷卻液流道6個(gè)熱儲(chǔ)能單元,因此,3個(gè)電堆的MFCS熱管理子系統(tǒng)(被控對(duì)象)的熱儲(chǔ)能單元數(shù)量為18??紤]到小電流工況時(shí),電堆內(nèi)空氣流道與冷卻液流道的熱儲(chǔ)能功能較弱,模型的階數(shù)根據(jù)模型辨識(shí)的結(jié)果可以小于18。

    以80A 工作點(diǎn)處為例,模型辨識(shí)的結(jié)果分別由式(19)-(21)表示,相應(yīng)的狀態(tài)空間維數(shù)為17。

    MFCS中各堆的電功率輸出相近時(shí)才能夠達(dá)到系統(tǒng)效率最優(yōu)[17-18]。因此,設(shè)定MFCS 中各堆工作功率相差不超過15kW,并且按照3個(gè)電堆負(fù)載電流的平均值大小進(jìn)行模型切換。模型切換條件如表2所示。

    表2 MPC模式選擇Tab.2 MPC mode selection

    2.3 控制器參數(shù)設(shè)定

    MPC控制器的采樣時(shí)間Ts根據(jù)系統(tǒng)特性確定,由于被控的溫度過程時(shí)滯性大,取采樣時(shí)間Ts=24s。預(yù)測時(shí)域P對(duì)控制算法耗時(shí)影響較大,控制時(shí)域M則對(duì)控制效果影響較大,不同的控制時(shí)域和預(yù)測時(shí)域?qū)刂扑惴ê臅r(shí)和控制效果的組合影響分別由圖5和圖6所示。

    圖5 不同預(yù)測和控制時(shí)域下的控制算法耗時(shí)Fig.5 Time consuming in different prediction and control time domains

    圖6 不同預(yù)測和控制時(shí)域下Tst1、Tst2、Tst3和Tstin的控制效果Fig.6 Control effects in different prediction and control time domains of Tst1,Tst2,Tst3,and Tstin

    綜合控制算法耗時(shí)和控制效果,最終確定控制器相應(yīng)的參數(shù)為Ts=24s,P=28,M=2。

    3 MFCS溫度控制仿真

    Simulink 是MATLAB 軟件中的一種可視化仿真工具,在Simulink 平臺(tái)上可以對(duì)動(dòng)態(tài)系統(tǒng)進(jìn)行建模、仿真和分析。本文在Simulink 平臺(tái)上對(duì)MFCS以及相應(yīng)的控制器進(jìn)行建模與仿真,分析系統(tǒng)的動(dòng)態(tài)特性和控制器的控制效果,仿真結(jié)果可以為MFCS熱管理控制提供一定的指導(dǎo)與借鑒。

    3.1 測試工況設(shè)計(jì)

    根據(jù)MFCS中相關(guān)部件工作特性的約束條件以及工作環(huán)境的約束條件,仿真過程中的各個(gè)約束與邊界條件如表3所示。

    表3 MPC控制器約束Tab.3 MPC controller constraints

    為了驗(yàn)證MPC控制效果,所設(shè)計(jì)的測試工況應(yīng)該能夠反映MFCS的大部分運(yùn)行狀況。其運(yùn)行工況有以下6種:①1、3個(gè)電堆中1個(gè)電堆的負(fù)載上升,其余2 個(gè)維持不變;②2、3 個(gè)電堆中2 個(gè)電堆的負(fù)載上升,其余1個(gè)維持不變;③3個(gè)電堆的負(fù)載同時(shí)上升;④3 個(gè)電堆中2 個(gè)電堆的負(fù)載下降,其余1 個(gè)維持不變;⑤3個(gè)電堆的負(fù)載同時(shí)下降;⑥3個(gè)電堆中1個(gè)電堆的負(fù)載下降,其余2個(gè)維持不變。

    根據(jù)上述運(yùn)行工況,分別設(shè)計(jì)MFCS 中3 個(gè)電堆的測試工況,組合后形成的熱管理控制器測試工況如圖7 所示。圖中分別標(biāo)出了對(duì)應(yīng)于6 個(gè)測試工況的典型,每個(gè)工況的開始時(shí)刻分別為2 500s、4 000s、5 000s、6 000s、6 500s、7 000s,持續(xù)時(shí)間均為500s。

    圖7 MFCS測試工況Fig.7 MFCS test working conditions

    3.2 仿真結(jié)果分析

    基于測試工況,利用MPC控制器對(duì)MFCS熱管理系統(tǒng)進(jìn)行控制仿真,仿真結(jié)果如圖8 所示。部分典型時(shí)間點(diǎn)的控制效果如表4 所示??梢钥闯?,在MPC 控制器的作用下,3 個(gè)電堆的溫度在整個(gè)測試工況中都維持在設(shè)定值(75℃)附近,超調(diào)量最大為1.10℃,穩(wěn)定時(shí)間最小為100s,最大為450s;入堆冷卻液溫度在整個(gè)測試工況下也維持在設(shè)定值(65℃)附近,超調(diào)量最大為1.30℃,穩(wěn)定時(shí)間最大為450s。由此可知,整體上MPC控制器對(duì)MFCS的溫度控制具有較好的效果。但是,在2 500s、4 000s、4 500s 和6 500s 附近控制效果不佳,超調(diào)量較大,穩(wěn)定時(shí)間較長。

    圖8 MPC控制仿真結(jié)果Fig.8 MPC control simulation results

    測試工況對(duì)應(yīng)的電堆平均電流、辨識(shí)電流和切換電流與時(shí)間的關(guān)系如圖9 所示。對(duì)比圖8 和9 可知,在2 500s、4 000s、4 500s 和6 500s 時(shí)間節(jié)點(diǎn)處控制效果較差是因?yàn)檫@3個(gè)點(diǎn)均處于預(yù)測模型切換邊界。模型切換會(huì)導(dǎo)致燃料電池系統(tǒng)的工作點(diǎn)距離模型辨識(shí)點(diǎn)較遠(yuǎn),預(yù)測模型與實(shí)際系統(tǒng)間差異較大,從而使得控制效果變差。另外,4 500s和6 500s處還存在整體電流變化較大的情況,較大的系統(tǒng)干擾會(huì)造成較大的超調(diào)量和較長的穩(wěn)定時(shí)間。

    圖9 測試工況電流信息Fig.9 Current information under test condition

    針對(duì)上述問題,對(duì)MPC控制器中的預(yù)測模型數(shù)目加以改進(jìn),增設(shè)120A和220A這2個(gè)系統(tǒng)辨識(shí)點(diǎn)。當(dāng)辨識(shí)點(diǎn)增多后,MPC控制器根據(jù)燃料電池的不同工況點(diǎn)可以更好地就近選擇預(yù)測模型,在一定程度上改善預(yù)測模型不精準(zhǔn)的缺點(diǎn)。

    對(duì)改進(jìn)后的系統(tǒng)采用同樣的測試工況進(jìn)行仿真,結(jié)果如圖10 所示,預(yù)測模型數(shù)目改進(jìn)后部分時(shí)間點(diǎn)的控制效果如表5 所示。對(duì)比圖8 和圖10 可知,在2 500s、4 000s 和6 500s 這3 個(gè)時(shí)間點(diǎn)處的控制效果發(fā)生了不同的變化。2 500s處電堆和冷卻液溫度的穩(wěn)定時(shí)間均減小了50s,電堆1 和電堆3 的超調(diào)量也分別減小了0.15℃和0.20℃,說明在改進(jìn)后的MPC控制器的作用下,電堆2的工況變化對(duì)電堆1和電堆3 的影響減??;4 000s處的電堆溫度和冷卻液溫度的穩(wěn)定時(shí)間都減小了50s,超調(diào)量最高可減小0.80℃;6 500s 點(diǎn)處冷卻液溫度的超調(diào)量減小了0.50℃。對(duì)控制器的改進(jìn)使得其對(duì)MFCS中各電堆的溫度控制效果得到了提升。

    圖10 改進(jìn)的MPC控制器仿真結(jié)果Fig.10 Simulation results of improved MPC controller

    表5 模型改進(jìn)后部分時(shí)間點(diǎn)控制效果Tab.5 Partial time point control effects after model improvement

    4 結(jié)語

    搭建了用于MFCS 各電堆溫度控制的并聯(lián)式熱管理子系統(tǒng)模型,引入PEMFC 陰極和陽極的出口氣體溫度作為修正項(xiàng)重新定義了電堆溫度,設(shè)計(jì)了MPC 控制器對(duì)MFCS 中各電堆溫度進(jìn)行控制,并利用測試工況進(jìn)行了仿真研究。研究結(jié)果顯示,MPC 控制器對(duì)MFCS 各電堆溫度的控制具有較好的效果,可為水冷型MFCS 的電堆溫度控制提供一定的參考。在MPC 控制器的作用下,MFCS中電堆溫度的波動(dòng)能夠維持在1.10℃以內(nèi),電堆溫度調(diào)節(jié)時(shí)間能夠維持在450s 以內(nèi)。對(duì)MPC 控制器進(jìn)行改進(jìn),將預(yù)測模型數(shù)目由3 個(gè)增加至5 個(gè),能夠緩解模型切換過程中預(yù)測模型與PEMFC 差異性較大的問題,進(jìn)一步改善控制效果,4 000s 處的電堆溫度波動(dòng)量最高能夠減少0.80℃,電堆溫度調(diào)節(jié)時(shí)間能夠減少50s。

    作者貢獻(xiàn)聲明:

    沈 偉:論文撰寫、數(shù)據(jù)分析。

    石 霖:論文撰寫、仿真研究。

    陳春光:論文撰寫、仿真研究。

    周 蘇:論文思路指導(dǎo)。

    王 寧:論文思路指導(dǎo)。

    裴馮來:論文思路指導(dǎo)及論文數(shù)據(jù)分析。

    猜你喜歡
    效果模型系統(tǒng)
    一半模型
    Smartflower POP 一體式光伏系統(tǒng)
    按摩效果確有理論依據(jù)
    WJ-700無人機(jī)系統(tǒng)
    ZC系列無人機(jī)遙感系統(tǒng)
    北京測繪(2020年12期)2020-12-29 01:33:58
    重要模型『一線三等角』
    重尾非線性自回歸模型自加權(quán)M-估計(jì)的漸近分布
    迅速制造慢門虛化效果
    抓住“瞬間性”效果
    中華詩詞(2018年11期)2018-03-26 06:41:34
    連通與提升系統(tǒng)的最后一塊拼圖 Audiolab 傲立 M-DAC mini
    国产视频首页在线观看| 国产精品人妻久久久久久| 国产av码专区亚洲av| 成人手机av| 国产国拍精品亚洲av在线观看| 99九九在线精品视频| 天堂俺去俺来也www色官网| 亚洲丝袜综合中文字幕| 久久免费观看电影| 九色成人免费人妻av| 人妻人人澡人人爽人人| 蜜桃国产av成人99| 亚洲无线观看免费| 亚洲色图 男人天堂 中文字幕 | 国产男女内射视频| a 毛片基地| 在线观看人妻少妇| 中文字幕人妻丝袜制服| 爱豆传媒免费全集在线观看| av不卡在线播放| 高清在线视频一区二区三区| 日本免费在线观看一区| 亚州av有码| 国产免费视频播放在线视频| 久久久久久久国产电影| 人人妻人人澡人人爽人人夜夜| 午夜av观看不卡| 国产精品一区二区在线不卡| 黄色视频在线播放观看不卡| 国产精品国产av在线观看| 亚洲精品aⅴ在线观看| 一本—道久久a久久精品蜜桃钙片| 亚洲成人手机| 精品亚洲成国产av| 人体艺术视频欧美日本| av免费在线看不卡| 一本一本综合久久| 日韩一区二区视频免费看| 满18在线观看网站| 精品久久国产蜜桃| 久久人人爽人人片av| 国产淫语在线视频| 丝瓜视频免费看黄片| 看免费成人av毛片| 亚洲国产av影院在线观看| 日日爽夜夜爽网站| 国产乱人偷精品视频| 女人久久www免费人成看片| 精品久久久久久久久亚洲| 在线观看美女被高潮喷水网站| 麻豆精品久久久久久蜜桃| 少妇被粗大的猛进出69影院 | 视频在线观看一区二区三区| 国产精品女同一区二区软件| 国产精品蜜桃在线观看| 国产在线免费精品| av在线观看视频网站免费| 国产成人精品婷婷| 日韩 亚洲 欧美在线| 黄色欧美视频在线观看| 两个人的视频大全免费| av天堂久久9| 久久精品国产亚洲网站| 日本黄色日本黄色录像| 亚洲怡红院男人天堂| 男男h啪啪无遮挡| 中文字幕亚洲精品专区| 99精国产麻豆久久婷婷| av又黄又爽大尺度在线免费看| 精品久久久精品久久久| 国产色婷婷99| 日韩一本色道免费dvd| 午夜久久久在线观看| 国产精品嫩草影院av在线观看| 亚洲欧美色中文字幕在线| 欧美三级亚洲精品| 丰满乱子伦码专区| 超碰97精品在线观看| 久久久久久久久久成人| 草草在线视频免费看| 丝袜美足系列| 黄色欧美视频在线观看| 久久久久久久久久久久大奶| 一级二级三级毛片免费看| 国产69精品久久久久777片| 中国美白少妇内射xxxbb| 日本与韩国留学比较| 亚洲av电影在线观看一区二区三区| 麻豆乱淫一区二区| 九草在线视频观看| 各种免费的搞黄视频| 欧美 亚洲 国产 日韩一| 亚洲av成人精品一区久久| 日本-黄色视频高清免费观看| 色婷婷av一区二区三区视频| 亚洲国产色片| 久久毛片免费看一区二区三区| 午夜日本视频在线| 高清不卡的av网站| 十八禁高潮呻吟视频| 成人亚洲欧美一区二区av| 亚洲欧美成人综合另类久久久| 王馨瑶露胸无遮挡在线观看| 秋霞伦理黄片| 日韩中字成人| 日韩亚洲欧美综合| 精品一品国产午夜福利视频| 日本午夜av视频| 免费观看无遮挡的男女| 国产亚洲av片在线观看秒播厂| 寂寞人妻少妇视频99o| 午夜激情久久久久久久| 青春草国产在线视频| 在线精品无人区一区二区三| 亚洲av成人精品一二三区| 一本久久精品| 热99久久久久精品小说推荐| 国产精品秋霞免费鲁丝片| 最后的刺客免费高清国语| av在线老鸭窝| 亚洲一区二区三区欧美精品| 中文字幕av电影在线播放| 久久精品国产a三级三级三级| 国产av码专区亚洲av| 免费久久久久久久精品成人欧美视频 | kizo精华| 日本wwww免费看| 插阴视频在线观看视频| 两个人免费观看高清视频| 免费高清在线观看视频在线观看| 91在线精品国自产拍蜜月| 蜜桃久久精品国产亚洲av| 国产熟女欧美一区二区| 精品一区在线观看国产| 纯流量卡能插随身wifi吗| 国产精品一区二区在线观看99| 少妇人妻 视频| 亚洲欧美清纯卡通| 欧美老熟妇乱子伦牲交| 国产成人免费观看mmmm| 久久99热6这里只有精品| 亚洲国产日韩一区二区| xxx大片免费视频| 在线免费观看不下载黄p国产| 精品少妇黑人巨大在线播放| 日韩中文字幕视频在线看片| 久久久精品区二区三区| 岛国毛片在线播放| 欧美精品国产亚洲| 人妻系列 视频| 黄色配什么色好看| 国产男人的电影天堂91| 亚洲国产欧美日韩在线播放| 满18在线观看网站| 久久国产亚洲av麻豆专区| 乱码一卡2卡4卡精品| 91久久精品电影网| 美女福利国产在线| videos熟女内射| 久久免费观看电影| 成人无遮挡网站| 男女无遮挡免费网站观看| 国产精品国产三级专区第一集| 啦啦啦视频在线资源免费观看| 久久久国产欧美日韩av| 亚洲色图 男人天堂 中文字幕 | 精品人妻熟女av久视频| 精品人妻熟女毛片av久久网站| 女人精品久久久久毛片| 亚洲欧美中文字幕日韩二区| 夫妻性生交免费视频一级片| 国产高清国产精品国产三级| 久久精品国产亚洲av天美| 久久99一区二区三区| 18禁在线播放成人免费| 国产成人精品一,二区| 亚洲精品亚洲一区二区| 99九九在线精品视频| 久久久精品免费免费高清| 69精品国产乱码久久久| 亚洲国产精品成人久久小说| 性高湖久久久久久久久免费观看| 成人亚洲欧美一区二区av| 国产亚洲欧美精品永久| 午夜精品国产一区二区电影| 满18在线观看网站| 毛片一级片免费看久久久久| 久久毛片免费看一区二区三区| 日韩精品有码人妻一区| 少妇 在线观看| 伦理电影大哥的女人| 五月玫瑰六月丁香| 狠狠精品人妻久久久久久综合| 中文字幕人妻丝袜制服| 成人漫画全彩无遮挡| 亚洲精品乱码久久久v下载方式| 在线看a的网站| 亚洲情色 制服丝袜| 精品亚洲成a人片在线观看| 久久久午夜欧美精品| 黄色毛片三级朝国网站| 蜜桃在线观看..| 亚洲国产精品专区欧美| 丝袜在线中文字幕| 国产永久视频网站| 精品久久久精品久久久| 岛国毛片在线播放| 亚洲av电影在线观看一区二区三区| 伊人久久国产一区二区| 亚洲精品久久久久久婷婷小说| 91精品国产国语对白视频| 成年女人在线观看亚洲视频| 日韩制服骚丝袜av| 欧美3d第一页| 精品人妻偷拍中文字幕| 亚洲高清免费不卡视频| 国产毛片在线视频| 欧美国产精品一级二级三级| 亚洲欧美色中文字幕在线| 成人黄色视频免费在线看| 亚洲av免费高清在线观看| 国产成人一区二区在线| 少妇的逼好多水| 18禁观看日本| 91午夜精品亚洲一区二区三区| 丝袜美足系列| 尾随美女入室| 免费av中文字幕在线| 亚洲精品乱码久久久v下载方式| 国产永久视频网站| 2021少妇久久久久久久久久久| 伦精品一区二区三区| 国产欧美另类精品又又久久亚洲欧美| 熟妇人妻不卡中文字幕| 亚洲精品久久久久久婷婷小说| av电影中文网址| 亚洲久久久国产精品| 欧美精品国产亚洲| 涩涩av久久男人的天堂| 欧美日韩在线观看h| 国产精品久久久久成人av| 亚洲,欧美,日韩| 亚洲综合色网址| 国产爽快片一区二区三区| 母亲3免费完整高清在线观看 | 亚洲国产精品专区欧美| 亚洲av不卡在线观看| 边亲边吃奶的免费视频| 亚洲精品国产av成人精品| 十八禁网站网址无遮挡| 青春草亚洲视频在线观看| 久久久久久久亚洲中文字幕| 一区二区av电影网| 999精品在线视频| 交换朋友夫妻互换小说| 王馨瑶露胸无遮挡在线观看| 九草在线视频观看| 观看美女的网站| 亚洲第一av免费看| 日本wwww免费看| 性色av一级| 大话2 男鬼变身卡| videossex国产| 国产毛片在线视频| av天堂久久9| 免费少妇av软件| 日韩不卡一区二区三区视频在线| 香蕉精品网在线| 久久热精品热| 日韩亚洲欧美综合| 久久人妻熟女aⅴ| 国产黄片视频在线免费观看| 中文字幕久久专区| 国产精品国产三级国产av玫瑰| 999精品在线视频| 日韩av不卡免费在线播放| 欧美日韩国产mv在线观看视频| 人妻一区二区av| 一级a做视频免费观看| 91成人精品电影| 26uuu在线亚洲综合色| 人妻系列 视频| 我的老师免费观看完整版| 国产精品久久久久久av不卡| 不卡视频在线观看欧美| 欧美国产精品一级二级三级| 在线观看免费视频网站a站| 久久久久久久久大av| 超碰97精品在线观看| 天天躁夜夜躁狠狠久久av| 精品久久久精品久久久| 大陆偷拍与自拍| 黄色怎么调成土黄色| 日韩欧美一区视频在线观看| 18禁裸乳无遮挡动漫免费视频| 熟妇人妻不卡中文字幕| 国产极品粉嫩免费观看在线 | 我的老师免费观看完整版| 国产成人a∨麻豆精品| 国产精品不卡视频一区二区| 一级毛片电影观看| 国产av一区二区精品久久| 亚洲av.av天堂| 精品99又大又爽又粗少妇毛片| 大香蕉久久成人网| 水蜜桃什么品种好| 国产成人精品无人区| 精品久久久精品久久久| 最近中文字幕2019免费版| 日日爽夜夜爽网站| 久久鲁丝午夜福利片| 黑人欧美特级aaaaaa片| 在线播放无遮挡| 一级片'在线观看视频| 十八禁高潮呻吟视频| 精品国产乱码久久久久久小说| 天堂俺去俺来也www色官网| 久久免费观看电影| 97精品久久久久久久久久精品| 久久99一区二区三区| 日韩免费高清中文字幕av| 伦理电影大哥的女人| 日本av手机在线免费观看| 欧美变态另类bdsm刘玥| 日本猛色少妇xxxxx猛交久久| tube8黄色片| 九九久久精品国产亚洲av麻豆| 亚洲精品456在线播放app| 欧美日韩成人在线一区二区| 久久久精品区二区三区| 国产精品无大码| 丝袜脚勾引网站| 免费观看在线日韩| 美女cb高潮喷水在线观看| .国产精品久久| 日日爽夜夜爽网站| videosex国产| av卡一久久| 免费大片黄手机在线观看| 18禁在线播放成人免费| 欧美激情国产日韩精品一区| 免费观看性生交大片5| 日韩成人伦理影院| 成人影院久久| 热99久久久久精品小说推荐| 搡老乐熟女国产| 亚洲人成网站在线观看播放| av天堂久久9| 国产精品国产av在线观看| 2018国产大陆天天弄谢| 久久精品熟女亚洲av麻豆精品| 爱豆传媒免费全集在线观看| 黄色欧美视频在线观看| 午夜激情久久久久久久| 一个人免费看片子| 自拍欧美九色日韩亚洲蝌蚪91| 国产日韩欧美亚洲二区| 国产成人免费观看mmmm| 国产日韩一区二区三区精品不卡 | 亚洲国产av影院在线观看| 99热这里只有精品一区| www.色视频.com| 内地一区二区视频在线| 免费高清在线观看视频在线观看| 国产69精品久久久久777片| av在线观看视频网站免费| 少妇熟女欧美另类| 国产一区二区三区av在线| 国产无遮挡羞羞视频在线观看| 国产精品成人在线| 美女视频免费永久观看网站| 一区在线观看完整版| 九色亚洲精品在线播放| 男女免费视频国产| 国产精品久久久久久精品电影小说| 国产免费福利视频在线观看| 黄色视频在线播放观看不卡| 日韩av不卡免费在线播放| 亚洲欧美中文字幕日韩二区| 久久免费观看电影| 韩国高清视频一区二区三区| 日本av免费视频播放| 国产精品一二三区在线看| 久久99蜜桃精品久久| 热re99久久精品国产66热6| 汤姆久久久久久久影院中文字幕| 麻豆精品久久久久久蜜桃| 大香蕉久久成人网| 特大巨黑吊av在线直播| www.色视频.com| 国产亚洲精品第一综合不卡 | 免费观看的影片在线观看| 欧美日韩一区二区视频在线观看视频在线| 在线 av 中文字幕| 日韩av在线免费看完整版不卡| 九色亚洲精品在线播放| 在线观看免费日韩欧美大片 | 一本—道久久a久久精品蜜桃钙片| 久久久国产欧美日韩av| 高清av免费在线| 国产淫语在线视频| 搡老乐熟女国产| 色5月婷婷丁香| 国产成人av激情在线播放 | 欧美少妇被猛烈插入视频| 少妇的逼水好多| 一级二级三级毛片免费看| 99国产精品免费福利视频| 久久精品久久精品一区二区三区| 亚洲人成77777在线视频| 精品久久久精品久久久| 交换朋友夫妻互换小说| 青春草国产在线视频| 免费高清在线观看视频在线观看| 国产精品熟女久久久久浪| 性高湖久久久久久久久免费观看| 黑人欧美特级aaaaaa片| 国产成人a∨麻豆精品| 99九九线精品视频在线观看视频| av在线app专区| 亚洲美女视频黄频| 成人二区视频| 一区二区三区精品91| 国产精品99久久99久久久不卡 | 狠狠婷婷综合久久久久久88av| 一区二区三区四区激情视频| a级片在线免费高清观看视频| 国国产精品蜜臀av免费| 日韩精品免费视频一区二区三区 | 九九久久精品国产亚洲av麻豆| 亚洲人成网站在线播| 成人黄色视频免费在线看| 国产成人aa在线观看| 伊人亚洲综合成人网| 国产精品一区www在线观看| 欧美三级亚洲精品| 精品久久久精品久久久| 成人亚洲精品一区在线观看| 日日摸夜夜添夜夜爱| 国产成人一区二区在线| 久久综合国产亚洲精品| 考比视频在线观看| 多毛熟女@视频| 精品视频人人做人人爽| 国产男人的电影天堂91| tube8黄色片| 久久久久久久久大av| 99久久人妻综合| 五月伊人婷婷丁香| 亚洲五月色婷婷综合| 亚洲精品视频女| 国产成人aa在线观看| 交换朋友夫妻互换小说| 日日爽夜夜爽网站| 91精品三级在线观看| 中国国产av一级| 午夜影院在线不卡| 亚洲av男天堂| 丝袜在线中文字幕| 18在线观看网站| 国产成人aa在线观看| 在线看a的网站| 亚洲四区av| .国产精品久久| 美女福利国产在线| 91aial.com中文字幕在线观看| 一本一本综合久久| 高清黄色对白视频在线免费看| 亚洲精品成人av观看孕妇| 精品久久久久久久久亚洲| 精品亚洲成国产av| 亚洲精品国产色婷婷电影| 欧美人与善性xxx| 精品卡一卡二卡四卡免费| 亚洲国产精品一区二区三区在线| 久久精品国产亚洲网站| 日本wwww免费看| 青春草视频在线免费观看| 99九九线精品视频在线观看视频| 波野结衣二区三区在线| 熟女人妻精品中文字幕| 日本色播在线视频| 少妇精品久久久久久久| 国产视频首页在线观看| 亚洲成色77777| 简卡轻食公司| 欧美人与性动交α欧美精品济南到 | 永久免费av网站大全| 丰满饥渴人妻一区二区三| 亚洲欧美成人综合另类久久久| 久久久久久久亚洲中文字幕| 9色porny在线观看| 国产av码专区亚洲av| 免费黄色在线免费观看| 国产精品秋霞免费鲁丝片| 精品视频人人做人人爽| 91在线精品国自产拍蜜月| 精品熟女少妇av免费看| 三级国产精品片| 最新中文字幕久久久久| 免费人成在线观看视频色| 日韩成人av中文字幕在线观看| 在线天堂最新版资源| 一级黄片播放器| 少妇丰满av| 美女中出高潮动态图| 18禁在线播放成人免费| 在线天堂最新版资源| 久久精品国产亚洲网站| 日日撸夜夜添| 熟女av电影| 国产免费一区二区三区四区乱码| 亚洲精品亚洲一区二区| 99久久人妻综合| 国国产精品蜜臀av免费| videosex国产| 日韩 亚洲 欧美在线| 最新中文字幕久久久久| 午夜免费观看性视频| 中文字幕av电影在线播放| 两个人的视频大全免费| 亚洲av不卡在线观看| 午夜激情久久久久久久| 九草在线视频观看| 欧美日韩av久久| 精品久久久久久久久亚洲| 中文字幕久久专区| 高清不卡的av网站| 母亲3免费完整高清在线观看 | 精品少妇内射三级| 久久免费观看电影| 亚洲av在线观看美女高潮| 一级毛片电影观看| 国产高清不卡午夜福利| 51国产日韩欧美| 搡女人真爽免费视频火全软件| 成人午夜精彩视频在线观看| 少妇猛男粗大的猛烈进出视频| 两个人的视频大全免费| 亚洲国产欧美日韩在线播放| 少妇的逼好多水| 伦理电影大哥的女人| 日韩 亚洲 欧美在线| 草草在线视频免费看| 久久久欧美国产精品| 2018国产大陆天天弄谢| 国产成人精品福利久久| 男人爽女人下面视频在线观看| 少妇被粗大猛烈的视频| h视频一区二区三区| 卡戴珊不雅视频在线播放| 亚洲av欧美aⅴ国产| 女的被弄到高潮叫床怎么办| 色哟哟·www| 少妇高潮的动态图| a 毛片基地| 亚洲怡红院男人天堂| 亚洲欧美成人精品一区二区| 婷婷色综合www| 日本av免费视频播放| 夜夜爽夜夜爽视频| 久久国产亚洲av麻豆专区| 男女边摸边吃奶| 极品少妇高潮喷水抽搐| 久久人人爽av亚洲精品天堂| 亚洲第一区二区三区不卡| 欧美成人精品欧美一级黄| 一边亲一边摸免费视频| 色94色欧美一区二区| 麻豆乱淫一区二区| 日韩三级伦理在线观看| 黄色一级大片看看| 午夜激情福利司机影院| 99热网站在线观看| 欧美激情国产日韩精品一区| a级毛片黄视频| 草草在线视频免费看| 国产精品不卡视频一区二区| av不卡在线播放| 少妇丰满av| 免费看不卡的av| 中文字幕精品免费在线观看视频 | 啦啦啦中文免费视频观看日本| av网站免费在线观看视频| 国产在线一区二区三区精| 久久人人爽av亚洲精品天堂| 在线看a的网站| 蜜桃久久精品国产亚洲av| av天堂久久9| 国产成人精品久久久久久| 晚上一个人看的免费电影| 久久久精品94久久精品| 欧美日韩一区二区视频在线观看视频在线| 久热久热在线精品观看| 欧美另类一区| 欧美日韩一区二区视频在线观看视频在线| 飞空精品影院首页| a级毛片免费高清观看在线播放| 国产熟女欧美一区二区| 亚洲国产欧美日韩在线播放| 激情五月婷婷亚洲| 美女脱内裤让男人舔精品视频| 欧美亚洲日本最大视频资源| 多毛熟女@视频| 日本免费在线观看一区| 三上悠亚av全集在线观看| 在线观看国产h片| 人妻夜夜爽99麻豆av| 亚洲精品一区蜜桃| 人人妻人人添人人爽欧美一区卜| 天天操日日干夜夜撸| 久久精品国产a三级三级三级| 在线观看美女被高潮喷水网站| 22中文网久久字幕| 国产日韩欧美亚洲二区|