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

    基于魯棒H∞濾波的鋰離子電池SOC估計(jì)

    2021-05-19 09:48:32潘鳳文弓棟梁寇亞林
    工程科學(xué)學(xué)報(bào) 2021年5期
    關(guān)鍵詞:方根值魯棒等效電路

    潘鳳文,弓棟梁,高 瑩?,寇亞林

    1) 吉林大學(xué)汽車仿真與控制國家重點(diǎn)實(shí)驗(yàn)室,長春 130025 2) 吉林大學(xué)汽車工程學(xué)院,長春 130025

    鋰離子電池能夠提供較高的能量密度和比功率,沒有記憶效應(yīng),具有較低的自放電率和較長的可循環(huán)使用壽命,且后續(xù)回收過程中對環(huán)境污染相對較少,這些特性使其成為目前新能源汽車的重要?jiǎng)恿υ碵1?3]. 但隨著鋰離子電池在新能源汽車中進(jìn)一步廣泛而深入的使用,其可靠性和安全性正在受越來越嚴(yán)苛的考驗(yàn). 電池管理系統(tǒng)(Battery management system,BMS)的核心目標(biāo)就是獲得表征鋰電池可用容量的SOC,保證電動(dòng)車輛安全有效的運(yùn)行,并通過降低電壓極限的不確定性來實(shí)現(xiàn)動(dòng)力分配的優(yōu)化,從而延長電池的使用壽命[4?6].然而,目前鋰離子電池可用的測量是電流、電壓和溫度,SOC不能直接測量[7?8]. 因此,準(zhǔn)確估計(jì)SOC對保證鋰電池的安全性和可靠性至關(guān)重要.

    對于鋰離子電池SOC的估計(jì)方法,目前可大致分為兩類:無模型法(安時(shí)積分法、開路電壓法)和基于模型的方法[9?13]. 基于模型的SOC估計(jì)中,常用的建模方式有黑箱模型、電化學(xué)模型、等效電路模型等. 黑箱模型通常使用傳遞函數(shù)來描述電池的行為,無需重建其基本的物理化學(xué)過程,但是該模型應(yīng)用場景受限,且模型精度不高. 電化學(xué)模型采用一組偏微分方程來構(gòu)建整個(gè)物理化學(xué)過程的“時(shí)?空”動(dòng)態(tài)變化,該模型中大量復(fù)雜的電化學(xué)方程計(jì)算耗時(shí)長,主要用于優(yōu)化電池的物理設(shè)計(jì)和電池相關(guān)的設(shè)計(jì)參數(shù),表征電能產(chǎn)生的基本機(jī)理等. 等效電路模型一般可分為Thevenin模型和阻抗模型兩大類. 這些基于電壓源,電阻和電容組合的模型具有直觀、實(shí)用和易于處理的特點(diǎn),通常用于控制設(shè)計(jì),且特別適合與其它電路或系統(tǒng)進(jìn)行協(xié)同設(shè)計(jì)和協(xié)同仿真. Thevenin模型又可分為一階RC等效電路模型,考慮電壓滯后的增強(qiáng)自校正模型,以及二階和高階RC模型等.

    在有關(guān)實(shí)時(shí)獲取電池SOC的文獻(xiàn)中,廣泛采用基于模型的SOC估計(jì)方法進(jìn)行研究. 文獻(xiàn)[14]采用二階RC等效電路模型,結(jié)合安時(shí)積分法與基于模型的方法對鋰電池的SOC進(jìn)行估計(jì),其結(jié)果表明,該方法增強(qiáng)了對測量噪聲的魯棒性. 文獻(xiàn)[15]~[16]針對離散化的鋰離子電池二階RC等效電路模型,采用線性參變技術(shù),設(shè)計(jì)了計(jì)算資源需求較少的SOC估計(jì)器,而且其所設(shè)計(jì)的估計(jì)器的穩(wěn)定性可由線性矩陣不等式(Linear matrix inequality,LMI)技術(shù)來分析驗(yàn)證. 文獻(xiàn)[17]采用大電容二階RC電池模型,該模型可用狀態(tài)空間函數(shù)來表示,由大容量電容的電壓值和電容值來確定SOC. 并在此模型基礎(chǔ)上設(shè)計(jì)了對未知擾動(dòng)和模型不確定性具有一定魯棒性的H∞濾波器;其所設(shè)計(jì)的H∞濾波器對SOC具有很好的估計(jì)性能. 文獻(xiàn)[18]采用與文獻(xiàn)[17]中相同的等效電路模型,相較于文獻(xiàn)[17],文獻(xiàn)[18]中設(shè)計(jì)了能夠改善觀測器響應(yīng)速度的H∞切換觀測器,進(jìn)一步提升了對SOC的估計(jì)性能. 文獻(xiàn)[9]和[19]均采用電化學(xué)鋰離子電池模型,并基于該模型設(shè)計(jì)了用于估計(jì)SOC的具有常值增益的龍伯格非線性觀測器,以及一個(gè)改進(jìn)的輸出雅可比矩陣加權(quán)增益的非線性觀測器,其中觀測器增益的設(shè)計(jì)均被轉(zhuǎn)化為了求解LMI問題.文獻(xiàn)[20]采用二階RC等效電路模型,把電池老化建模為加性故障,并對鋰離子電池模型進(jìn)行了線性化,設(shè)計(jì)了基于H∞觀測器的故障估計(jì)器,從其仿真結(jié)果來看,所設(shè)計(jì)的故障估計(jì)器對干擾具有很好的魯棒性,同時(shí)對老化故障具有很強(qiáng)的追蹤能力. 文獻(xiàn)[21]采用與文獻(xiàn)[20]中相同的鋰離子電池模型,針對SOC的估計(jì)問題,采用LMI技術(shù)設(shè)計(jì)了基于H∞方法的非線性觀測器. 文獻(xiàn)[22]針對鋰離子電池SOC的估計(jì)誤差收斂問題,設(shè)計(jì)了三種基于LMI技術(shù)的魯棒龍伯格SOC觀測器,并分別與傳統(tǒng)的基于非線性等效電路模型的擴(kuò)展卡爾曼濾波器進(jìn)行對比,對比結(jié)果表明,其設(shè)計(jì)的兩種基于LMI的SOC觀測器與廣泛使用的擴(kuò)展卡爾曼濾波器具有相近的性能,并且設(shè)計(jì)方法更為簡單. 文獻(xiàn)[23]~[24]均采用非線性分?jǐn)?shù)階鋰離子電池模型對SOC進(jìn)行估計(jì),不同之處是前者采用二階RC等效電路模型,并基于SOC估計(jì)誤差動(dòng)力學(xué)方程,運(yùn)用LMI技術(shù)來設(shè)計(jì)SOC觀測器;后者采用基于阻抗模型和分?jǐn)?shù)階卡爾曼濾波方法對SOC進(jìn)行估計(jì).

    綜上可知,大多數(shù)文獻(xiàn)在鋰離子電池SOC估計(jì)中選擇與LMI技術(shù)相結(jié)合,這樣不僅可以簡化設(shè)計(jì)過程,而且也保證了SOC的估計(jì)準(zhǔn)確性和一定的魯棒性. 然而,針對鋰離子電池模型參數(shù)隨溫度變動(dòng)而產(chǎn)生攝動(dòng),進(jìn)而影響SOC估計(jì)準(zhǔn)確性的問題討論的并不多. 因此,本文在前述工作基礎(chǔ)之上,考慮電池溫度波動(dòng)對模型參數(shù)的影響,把溫度對模型參數(shù)的影響方式建模為標(biāo)稱電阻值和電池總?cè)萘康募有宰兞浚⒁暅囟茸兓癁橄到y(tǒng)的外部擾動(dòng). 然后,采用切線法對模型進(jìn)行線性化. 繼而采用LMI技術(shù)設(shè)計(jì)了對SOC進(jìn)行估計(jì)的魯棒H∞濾波器. 最后,采用四種不同類型的動(dòng)態(tài)激勵(lì)電流,并與kalman濾波對SOC的估計(jì)進(jìn)行對比,對所設(shè)計(jì)的SOC觀測器進(jìn)行驗(yàn)證.

    1 鋰離子電池建模

    鋰離子電池的動(dòng)態(tài)特性高度依賴于它的工作條件,具有時(shí)變參數(shù)的二階RC等效電路模型可以足夠精確的來表征鋰電池系統(tǒng)特性,同時(shí)又相對簡單且計(jì)算量不大. 因此,本文主要采用二階RC等效電路模型.

    1.1 電氣模型

    表征鋰離子電池動(dòng)態(tài)特性的二階RC等效電路模型如圖1所示.

    圖1 等效電路模型Fig.1 Equivalent circuit model

    根據(jù)基爾霍夫電壓定律,鋰離子電池的二階RC等效電路模型的動(dòng)力學(xué)方程可表述為

    式中:t代表時(shí)間,s;Uocv代表與SOC相關(guān)的開路電壓(其為SOC的函數(shù)),V;Ic為進(jìn)入或離開電池正極的電流(放電時(shí)為正,充電時(shí)為負(fù)),A;R0為內(nèi)阻,Ω;R1和C1分別為電化學(xué)極化電阻(Ω)和電容(F);R2和C2分別為濃差極化電阻(Ω)和電容(F);U1和U2分別代表電化學(xué)極化電容C1和濃差極化電容C2的電壓,V;UT為端電壓,V.

    1.2 電池荷電狀態(tài)模型

    電池SOC的定義為可用容量(A·h)與電池的總?cè)萘康谋戎? 簡化的SOC估計(jì)模型如下:

    式中,Q為電池總?cè)萘?,SOC0為電池的初始SOC.需要指出的是電池容量也會(huì)隨著溫度的升高而變化,從而影響SOC估計(jì)的準(zhǔn)確性[13].

    大量試驗(yàn)表明,開路電壓Uocv與SOC在較大的工作區(qū)間內(nèi)均存在一定的線性關(guān)系(除了SOC非常高或非常低的情況,小于20%或大于80%)[14?15].因此,可以在特定的工作范圍采用SOC的線性函數(shù)來近似逼近開路電壓. 它們之間的函數(shù)關(guān)系具有以下形式:

    式中:a和b為給定的常量.

    那么根據(jù)公式(4),終端電壓UT與SOC的關(guān)系可表示為:

    所選用的電池模型參數(shù)如下表1所示[20]. 該模型參數(shù)的獲取需要滿足以下的條件:(1)T=40 ℃,Ic=1 A,SOC=0.5;(2)電池模型的參數(shù)與溫度、SOC、以及電流的大小和方向均無關(guān).

    表1 電池模型參數(shù)Table 1 Battery model parameters

    上述條件假設(shè)非常嚴(yán)格,大量的實(shí)際實(shí)驗(yàn)已表明:電池模型參數(shù)對于諸如溫度、SOC和充放電狀態(tài)以及電流等環(huán)境狀況的變化是很敏感的;其中,溫度對電池模型參數(shù)的影響最為顯著[16].

    2 模型線性化

    由于溫度對鋰離子電池的模型參數(shù)具有較大的影響,因此,本文將溫度對模型參數(shù)的影響方式建模為標(biāo)稱電阻值和電池總?cè)萘康募有宰兞? 其具有如下形式:

    式中:Rni(i=0,1,2)分別表示,,標(biāo)稱阻值;ΔRi(i=0,1,2)表示由溫度變化而引起的阻值變化;R?i(i=0,1,2)表示由溫度變化而產(chǎn)生的最新阻值;表示由溫度變化引起的電池容量變化,表示由溫度變化產(chǎn)生的最新容量.

    3 魯棒 SOC 估計(jì)器設(shè)計(jì)

    針對線性化系統(tǒng)(15),設(shè)計(jì)一個(gè)漸近穩(wěn)定的H∞濾波器,實(shí)現(xiàn)對SOC的估計(jì),使得傳遞函數(shù)(從外部擾動(dòng)到SOC估計(jì)誤差)的H∞范數(shù)小于給定值,該濾波器具有以下狀態(tài)空間形式:

    4 仿真實(shí)驗(yàn)

    通過在Matlab中求解優(yōu)化問題(29),獲得干擾抑制性能指標(biāo) γ =5.7274×10?5. 最終獲得所設(shè)計(jì)的魯棒H∞濾波器狀態(tài)空間矩陣分別為

    仿真實(shí)驗(yàn)中給電流輸入信號注入了能量為0.001的白噪聲,溫度變動(dòng)用幅值為40的正弦函數(shù)進(jìn)行模擬,仿真結(jié)果與kalman濾波對SOC的估計(jì)結(jié)果進(jìn)行對比分析. 采用四種類型的動(dòng)態(tài)測試進(jìn)行仿真實(shí)驗(yàn)驗(yàn)證,其分別為 BJDST (Beijing Dynamic Stress Test)、FUDS (Federal Urban Driving Schedule)、US06 (US06 Highway Driving Schedule)以及 BJDST?FUDS?US06聯(lián)合駕駛時(shí)間表動(dòng)態(tài)測試. 一個(gè)完整BJDST 循環(huán)是 916 s,一個(gè)完整 FUDS 循環(huán)是 1372 s,一個(gè)完整US06循環(huán)是600 s. 其完整循環(huán)時(shí)間歷程如下圖2,圖3所示.

    圖2 BJDST,F(xiàn)UDS,US06動(dòng)態(tài)測試時(shí)間歷程Fig.2 BJDST, FUDS, US06 dynamic test time history

    圖3 BJDST?FUDS?US06聯(lián)合動(dòng)態(tài)測試時(shí)間歷程Fig.3 BJDST?FUDS?US06 joint dynamic test time history

    BJDST,F(xiàn)UDS和US06動(dòng)態(tài)激勵(lì)下的基于H∞濾波的SOC估計(jì)和基于kalman濾波的SOC估計(jì)以及它們的SOC估計(jì)誤差時(shí)間歷程曲線如圖4~圖9所示. 其中藍(lán)色實(shí)線表示模型輸出,紅色實(shí)線表示H∞估計(jì),綠色實(shí)線表示kalman估計(jì).

    圖4 BJDST 激勵(lì)的 SOC 估計(jì)Fig.4 SOC estimation excited by BJDST

    圖5 BJDST 激勵(lì)的 SOC 估計(jì)誤差Fig.5 SOC estimation error excited by BJDST

    圖6 FUDS 激勵(lì)的 SOC 估計(jì)Fig.6 SOC estimation excited by FUDS

    圖7 FUDS 激勵(lì)的 SOC 估計(jì)誤差Fig.7 SOC estimation error excited by FUDS

    圖8 US06激勵(lì)的 SOC 估計(jì)Fig.8 SOC estimation excited by US06

    圖9 US06激勵(lì)的 SOC 估計(jì)誤差Fig.9 SOC estimation error excited by US06

    從圖4中可以看出,在BJDST動(dòng)態(tài)激勵(lì)下,基于魯棒H∞濾波器的SOC估計(jì)時(shí)間歷程曲線比基于kalman濾波的SOC估計(jì)時(shí)間歷程曲線更為接近模型輸出值. 圖5中BJDST動(dòng)態(tài)激勵(lì)下的SOC估計(jì)誤差表明,基于魯棒H∞濾波器的SOC估計(jì)誤差波動(dòng)幅值小于基于kalman濾波的SOC估計(jì)誤差波動(dòng)幅值. 從表2中對應(yīng)的均方根值對比結(jié)果可以看出,基于魯棒H∞濾波器的SOC估計(jì)比模型輸出大0.1%,而基于kalman濾波的SOC估計(jì)相較于模型輸出小0.27%. 基于魯棒H∞濾波器的SOC估計(jì)誤差均方根值是0.0019,明顯小于基于kalman濾波的SOC估計(jì)誤差均方根值0.0041. 圖4,圖5以及對應(yīng)表2中的均方根值分析結(jié)果表明,即使在電流信號存在干擾,電池溫度持續(xù)波動(dòng)并對模型參數(shù)產(chǎn)生一定影響的條件下,所設(shè)計(jì)的魯棒H∞濾波器依然能夠?qū)崿F(xiàn)對SOC的準(zhǔn)確估計(jì).

    表2 SOC 估計(jì)及其估計(jì)誤差的均方根值Table 2 RMS value of SOC estimation and estimation error

    四種不同類型動(dòng)態(tài)激勵(lì)下的SOC估計(jì),以及SOC 估計(jì)誤差的均方根值(Root mean square,RMS)如下表2所示.

    BJDST?FUDS?US06聯(lián)合動(dòng)態(tài)測試下的SOC估計(jì)及其估計(jì)誤差時(shí)間歷程曲線如下圖10~11所示.

    圖10 BJDST?FUDS?US06激勵(lì)的 SOC 估計(jì)Fig.10 SOC estimation excited by BJDST?FUDS?US06

    圖11 BJDST?FUDS?US06 US06激勵(lì)的 SOC 估計(jì)誤差Fig.11 SOC estimation error excited by BJDST?FUDS?US06

    在FUDS動(dòng)態(tài)激勵(lì)下,圖6中的SOC估計(jì)時(shí)間歷程曲線表明相較于基于kalman濾波的SOC估計(jì),基于魯棒H∞濾波器的SOC估計(jì)與模型輸出較為接近. 圖7中的SOC估計(jì)誤差時(shí)間歷程曲線表明了基于魯棒H∞濾波器的SOC估計(jì)誤差波動(dòng)幅值比基于kalman濾波的SOC估計(jì)誤差波動(dòng)幅值要小. 從表2中FUDS動(dòng)態(tài)測試下的均方根值可以看出,基于魯棒H∞濾波器的SOC估計(jì)的均方根值與模型輸出的均方根值近乎一致,而基于kalman濾波的SOC估計(jì)的均方根值相較于模型輸出的均方根值小0.21%. 基于kalman濾波的SOC估計(jì)誤差均方根值0.0036明顯大于基于魯棒H∞濾波器的SOC估計(jì)誤差均方根值0.0018. 圖6,圖7以及表2中相對應(yīng)的FUDS激勵(lì)下的均方根值分析對比結(jié)果表明,所設(shè)計(jì)的H∞濾波器不僅具有很好的魯棒性,而且具有優(yōu)異的SOC估計(jì)性能.

    在US06動(dòng)態(tài)激勵(lì)下,圖8中的SOC時(shí)間歷程曲線和圖9中的SOC估計(jì)誤差時(shí)間歷程曲線表明,與基于kalman濾波的SOC估計(jì)相比,基于魯棒H∞濾波器的SOC估計(jì)更為接近模型輸出值.表2中對應(yīng)的US06動(dòng)態(tài)激勵(lì)下的均方根值對比結(jié)果表明,基于魯棒H∞濾波器的SOC估計(jì)的均方根值與模型輸出的均方根值幾近一致,而基于kalman濾波的SOC估計(jì)的均方根值比模型輸出的均方根值小0.1%. 基于魯棒H∞濾波器的SOC估計(jì)誤差均方根值是0.0019,明顯小于基于kalman濾波的SOC估計(jì)誤差均方根值0.0043. 仿真實(shí)驗(yàn)與均方根值的對比結(jié)果表明,所設(shè)計(jì)的H∞濾波器具有更為優(yōu)異的SOC估計(jì)性能,同時(shí)具有很好的干擾抑制效果,即具有很強(qiáng)的魯棒性.

    在BDST?FUDS?US06聯(lián)合動(dòng)態(tài)測試下,圖10中接近10 h的時(shí)間歷程表明,H∞濾波和kalman濾波的時(shí)間歷程均比較接近于模型輸出結(jié)果. 從圖11中的SOC估計(jì)誤差時(shí)間歷程曲線中可以看出,H∞估計(jì)誤差比kalman濾波估計(jì)誤差具有更小的波動(dòng)幅值. 表2中對應(yīng)的 BDST?FUDS?US06聯(lián)合動(dòng)態(tài)激勵(lì)下的均方根值表明,H∞濾波與模型輸出的均方根值均為0.3186,而kalman濾波比模型輸出的均方根值小0.35%. 基于kalman濾波的SOC估計(jì)誤差均方根值0.0033明顯大于基于魯棒H∞濾波器的SOC估計(jì)誤差均方根值0.0022. 仿真實(shí)驗(yàn)與均方根值的對比結(jié)果表明,所設(shè)計(jì)的魯棒H∞濾波器不僅對外部擾動(dòng)表現(xiàn)出了很好的魯棒性,而且能夠?qū)崿F(xiàn)對SOC的強(qiáng)跟蹤能力.

    5 結(jié)論

    本文采用LMI技術(shù)設(shè)計(jì)了用于SOC估計(jì)的魯棒H∞濾波器. 該濾波器的設(shè)計(jì)考慮了鋰離子電池模型參數(shù)隨溫度變化產(chǎn)生攝動(dòng)影響SOC估計(jì)準(zhǔn)確性的問題. 把電池溫度變化對模型參數(shù)的影響方式建模為標(biāo)稱電阻值和電池總?cè)萘康募有宰兞?,并將電池溫度變化看作系統(tǒng)的外部擾動(dòng). 通過切線法對模型進(jìn)行了線性化,并基于該線性化模型設(shè)計(jì)了估計(jì)SOC的觀測器. 最后,采用四種不同類型的動(dòng)態(tài)電流激勵(lì),并與kalman濾波對SOC的估計(jì)結(jié)果進(jìn)行仿真對比,其結(jié)果表明所設(shè)計(jì)的魯棒H∞濾波器不僅具有優(yōu)異的SOC估計(jì)性能,而且對外部擾動(dòng)表現(xiàn)出了很好的魯棒性.

    猜你喜歡
    方根值魯棒等效電路
    磁致伸縮振動(dòng)能量收集器的全耦合非線性等效電路模型
    基于撕裂法的變壓器繞組集總參數(shù)等效電路頻率響應(yīng)計(jì)算方法
    磁流變彈性體減振單元?jiǎng)恿W(xué)分析
    客車平順性仿真及優(yōu)化
    基于學(xué)習(xí)的魯棒自適應(yīng)評判控制研究進(jìn)展
    目標(biāo)魯棒識別的抗旋轉(zhuǎn)HDO 局部特征描述
    基于Cauchy魯棒函數(shù)的UKF改進(jìn)算法
    目標(biāo)軌跡更新的點(diǎn)到點(diǎn)魯棒迭代學(xué)習(xí)控制
    Studyof velocityfluctuations in the plenum of a 3/4open jet automotive wind tunnel
    變頻調(diào)速異步電動(dòng)機(jī)電磁轉(zhuǎn)矩計(jì)算方法
    久99久视频精品免费| 看非洲黑人一级黄片| 在线天堂最新版资源| 国产v大片淫在线免费观看| 日韩欧美精品免费久久| 亚洲av成人av| 亚洲欧美成人综合另类久久久 | 精品久久久久久久末码| 久久婷婷人人爽人人干人人爱| 日韩av不卡免费在线播放| av福利片在线观看| 亚洲成av人片在线播放无| 国产精品伦人一区二区| 午夜免费激情av| 激情 狠狠 欧美| 亚洲国产欧美在线一区| 亚洲精品成人久久久久久| 国产精品综合久久久久久久免费| 亚洲中文字幕日韩| 在线观看免费视频日本深夜| 精品久久久久久久久久免费视频| 午夜福利在线观看吧| 毛片一级片免费看久久久久| 我要搜黄色片| 有码 亚洲区| av免费观看日本| 男的添女的下面高潮视频| 在线a可以看的网站| 国产伦理片在线播放av一区 | 亚洲人成网站在线播| 午夜激情欧美在线| 国产人妻一区二区三区在| 国产亚洲av嫩草精品影院| 亚洲婷婷狠狠爱综合网| 中文字幕熟女人妻在线| 综合色av麻豆| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 老师上课跳d突然被开到最大视频| 干丝袜人妻中文字幕| a级毛色黄片| 免费在线观看成人毛片| 狂野欧美白嫩少妇大欣赏| 午夜老司机福利剧场| 嫩草影院新地址| 日韩强制内射视频| 精品国产三级普通话版| 亚洲自拍偷在线| 男女啪啪激烈高潮av片| 天天躁夜夜躁狠狠久久av| 熟女电影av网| 国产精品国产高清国产av| 人妻制服诱惑在线中文字幕| 亚洲精华国产精华液的使用体验 | eeuss影院久久| 又粗又爽又猛毛片免费看| 亚洲第一区二区三区不卡| 波野结衣二区三区在线| 黄色一级大片看看| 欧美日韩国产亚洲二区| 欧美激情国产日韩精品一区| 人人妻人人看人人澡| 啦啦啦韩国在线观看视频| 国产综合懂色| 久久久a久久爽久久v久久| 免费黄网站久久成人精品| 中文字幕人妻熟人妻熟丝袜美| 亚洲av不卡在线观看| 久久九九热精品免费| 成人性生交大片免费视频hd| 久久久久网色| 青春草亚洲视频在线观看| 日韩欧美精品免费久久| 99久久成人亚洲精品观看| 婷婷色av中文字幕| 欧美成人精品欧美一级黄| 欧美丝袜亚洲另类| 啦啦啦啦在线视频资源| 久久精品国产亚洲av天美| 精品人妻一区二区三区麻豆| 99热这里只有是精品50| 国产在线精品亚洲第一网站| www.av在线官网国产| 成人鲁丝片一二三区免费| 中文字幕人妻熟人妻熟丝袜美| 中文在线观看免费www的网站| 精品久久久久久久久久免费视频| 国产成人午夜福利电影在线观看| 成人永久免费在线观看视频| 免费av不卡在线播放| 国产女主播在线喷水免费视频网站 | 亚洲精华国产精华液的使用体验 | 亚洲综合色惰| 久久中文看片网| 国产精品人妻久久久影院| АⅤ资源中文在线天堂| 美女xxoo啪啪120秒动态图| 亚洲aⅴ乱码一区二区在线播放| 国产精品女同一区二区软件| 国模一区二区三区四区视频| 能在线免费看毛片的网站| 亚洲精品久久国产高清桃花| 一本一本综合久久| 免费在线观看成人毛片| 国产一级毛片在线| 国产久久久一区二区三区| 99热这里只有精品一区| 女的被弄到高潮叫床怎么办| 国产高清三级在线| 久久久久久久久中文| 成年版毛片免费区| 国产高潮美女av| 亚洲成人中文字幕在线播放| 自拍偷自拍亚洲精品老妇| 婷婷色综合大香蕉| 免费电影在线观看免费观看| 一区福利在线观看| 五月玫瑰六月丁香| 嫩草影院精品99| 国产黄片视频在线免费观看| 国产亚洲精品久久久com| 成年免费大片在线观看| 亚洲天堂国产精品一区在线| 欧美高清成人免费视频www| 搞女人的毛片| 波野结衣二区三区在线| 精品一区二区三区人妻视频| 精品人妻视频免费看| 99久久中文字幕三级久久日本| 一级毛片电影观看 | 久久久成人免费电影| 一级二级三级毛片免费看| 欧美日韩一区二区视频在线观看视频在线 | 美女 人体艺术 gogo| 午夜精品国产一区二区电影 | 成年女人永久免费观看视频| 嫩草影院新地址| 99在线视频只有这里精品首页| 久久人人爽人人片av| 在线播放国产精品三级| 欧美变态另类bdsm刘玥| 免费av不卡在线播放| 欧美又色又爽又黄视频| 亚洲综合色惰| 偷拍熟女少妇极品色| 嫩草影院入口| 国产亚洲欧美98| 直男gayav资源| 久久这里有精品视频免费| 国产免费男女视频| 不卡视频在线观看欧美| 欧美区成人在线视频| 1000部很黄的大片| 久久精品人妻少妇| 黄片无遮挡物在线观看| 国产精品一及| 欧美性感艳星| 午夜亚洲福利在线播放| 久久久久网色| 精品久久久久久久久久久久久| 人人妻人人看人人澡| 国产高清不卡午夜福利| 国产真实伦视频高清在线观看| 国产男人的电影天堂91| 嫩草影院新地址| 老司机福利观看| 大型黄色视频在线免费观看| 22中文网久久字幕| 亚洲av.av天堂| 一级二级三级毛片免费看| 美女cb高潮喷水在线观看| 免费在线观看成人毛片| 久久久a久久爽久久v久久| 亚洲精华国产精华液的使用体验 | 亚洲精品国产成人久久av| 两个人视频免费观看高清| 国产精品一区www在线观看| 国产v大片淫在线免费观看| 老女人水多毛片| 国产成人一区二区在线| 久久精品91蜜桃| av天堂在线播放| 长腿黑丝高跟| 亚洲婷婷狠狠爱综合网| 国产美女午夜福利| 国产日本99.免费观看| 国产又黄又爽又无遮挡在线| 国产午夜精品一二区理论片| 亚洲国产精品合色在线| 亚洲真实伦在线观看| 我的女老师完整版在线观看| 一级二级三级毛片免费看| 禁无遮挡网站| 欧洲精品卡2卡3卡4卡5卡区| 99热精品在线国产| 91狼人影院| 天天躁夜夜躁狠狠久久av| 男女下面进入的视频免费午夜| 久久精品综合一区二区三区| 色播亚洲综合网| 婷婷色av中文字幕| 国产激情偷乱视频一区二区| 亚洲最大成人中文| 日韩精品青青久久久久久| 最新中文字幕久久久久| 国产一区二区激情短视频| 午夜免费激情av| 校园春色视频在线观看| 在线国产一区二区在线| 两性午夜刺激爽爽歪歪视频在线观看| 在线a可以看的网站| 成人三级黄色视频| 国产精品一区二区三区四区久久| 成人午夜高清在线视频| 国产亚洲欧美98| 亚洲人成网站在线观看播放| 久久精品91蜜桃| 男人舔女人下体高潮全视频| 舔av片在线| 国产亚洲av嫩草精品影院| 成人亚洲精品av一区二区| 午夜免费男女啪啪视频观看| 久久精品91蜜桃| 最好的美女福利视频网| 日韩制服骚丝袜av| 国产高清有码在线观看视频| 午夜激情欧美在线| 精品人妻一区二区三区麻豆| 日韩精品有码人妻一区| 亚洲国产精品成人久久小说 | 国产伦在线观看视频一区| 亚洲精品久久国产高清桃花| 男人和女人高潮做爰伦理| 秋霞在线观看毛片| 亚洲av成人精品一区久久| 人妻制服诱惑在线中文字幕| 亚洲中文字幕一区二区三区有码在线看| 99久久精品热视频| 91久久精品国产一区二区三区| 欧美成人a在线观看| 国产亚洲欧美98| 99久久人妻综合| 在线国产一区二区在线| 国产亚洲精品久久久久久毛片| 最好的美女福利视频网| 男女边吃奶边做爰视频| 69人妻影院| 亚洲av免费高清在线观看| 亚洲人与动物交配视频| 午夜福利成人在线免费观看| 国产乱人偷精品视频| 老熟妇乱子伦视频在线观看| 99国产极品粉嫩在线观看| 久久久久久久久大av| 中文字幕熟女人妻在线| 性插视频无遮挡在线免费观看| 中文字幕精品亚洲无线码一区| 亚洲国产精品久久男人天堂| 男女那种视频在线观看| 黄色配什么色好看| av在线蜜桃| 亚洲欧洲国产日韩| 亚洲精品粉嫩美女一区| 国产精品一区二区三区四区久久| 国产一区二区激情短视频| 久久久精品94久久精品| 中国国产av一级| 国产成人影院久久av| 午夜福利成人在线免费观看| 久久精品国产亚洲av涩爱 | 日韩,欧美,国产一区二区三区 | av卡一久久| 韩国av在线不卡| 又黄又爽又刺激的免费视频.| h日本视频在线播放| 黄片无遮挡物在线观看| 舔av片在线| 人妻系列 视频| 亚洲丝袜综合中文字幕| 99久久久亚洲精品蜜臀av| 欧美bdsm另类| 亚洲最大成人av| 春色校园在线视频观看| 日韩欧美国产在线观看| 亚洲经典国产精华液单| 亚洲精品成人久久久久久| 中国国产av一级| 麻豆精品久久久久久蜜桃| av专区在线播放| 悠悠久久av| 国产精品一区二区三区四区久久| 午夜精品在线福利| 国产成人a区在线观看| 不卡一级毛片| 国产成人91sexporn| 亚洲五月天丁香| 亚洲av电影不卡..在线观看| 亚洲人成网站在线观看播放| 91午夜精品亚洲一区二区三区| 在线免费观看的www视频| 一本一本综合久久| 亚洲在线自拍视频| 夜夜看夜夜爽夜夜摸| 少妇被粗大猛烈的视频| а√天堂www在线а√下载| 色视频www国产| 免费看光身美女| 国产午夜精品久久久久久一区二区三区| 综合色av麻豆| 在现免费观看毛片| 十八禁国产超污无遮挡网站| 欧美日韩一区二区视频在线观看视频在线 | 国产成人影院久久av| 亚洲人成网站在线播| 变态另类成人亚洲欧美熟女| 久久亚洲精品不卡| 久久久久九九精品影院| 久久久久久久久久久丰满| 91狼人影院| 中文欧美无线码| 日本av手机在线免费观看| 亚洲av电影不卡..在线观看| www.av在线官网国产| 婷婷六月久久综合丁香| 在线观看午夜福利视频| 给我免费播放毛片高清在线观看| 久久久精品大字幕| 三级国产精品欧美在线观看| 亚洲国产欧美人成| 三级毛片av免费| 尾随美女入室| 国产精品综合久久久久久久免费| 亚洲熟妇中文字幕五十中出| 国产精品嫩草影院av在线观看| 成人欧美大片| 亚洲第一区二区三区不卡| 国产精品一区二区三区四区久久| 午夜免费激情av| 天堂av国产一区二区熟女人妻| 久久精品人妻少妇| 国产精品免费一区二区三区在线| 亚洲精品久久国产高清桃花| 亚洲av电影不卡..在线观看| 国产成人一区二区在线| 国产精品无大码| 欧美成人a在线观看| 观看免费一级毛片| 午夜福利在线观看免费完整高清在 | 国产精品一区二区三区四区久久| 久久久久九九精品影院| 久久精品国产鲁丝片午夜精品| 变态另类成人亚洲欧美熟女| 精品午夜福利在线看| 国产黄色视频一区二区在线观看 | 欧美高清性xxxxhd video| 九草在线视频观看| 日日撸夜夜添| 亚洲成人久久爱视频| 午夜免费男女啪啪视频观看| 久久欧美精品欧美久久欧美| 少妇裸体淫交视频免费看高清| 三级毛片av免费| 日韩国内少妇激情av| av免费观看日本| 97热精品久久久久久| 亚洲第一电影网av| 成人毛片60女人毛片免费| 亚洲av成人av| 国产老妇伦熟女老妇高清| 日韩视频在线欧美| 久久精品人妻少妇| 中国国产av一级| 亚洲欧美精品综合久久99| 欧美区成人在线视频| 久久久久久久午夜电影| 欧美成人免费av一区二区三区| 色综合色国产| 美女黄网站色视频| 午夜精品在线福利| 久久久久免费精品人妻一区二区| 97人妻精品一区二区三区麻豆| 青青草视频在线视频观看| 桃色一区二区三区在线观看| 国产精品不卡视频一区二区| 欧美日韩精品成人综合77777| 国产黄片视频在线免费观看| 级片在线观看| 两个人的视频大全免费| 淫秽高清视频在线观看| 亚洲av不卡在线观看| 欧美日韩综合久久久久久| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | av福利片在线观看| kizo精华| 国产乱人视频| 18禁在线播放成人免费| 中文字幕免费在线视频6| 国产熟女欧美一区二区| 亚洲国产欧美在线一区| 久久精品国产亚洲网站| 欧美人与善性xxx| 91麻豆精品激情在线观看国产| 国产亚洲欧美98| 在线天堂最新版资源| 国产成人福利小说| 啦啦啦啦在线视频资源| 国产老妇女一区| 色综合站精品国产| 国产精品久久电影中文字幕| 九九在线视频观看精品| 99在线人妻在线中文字幕| 国产精品久久久久久久久免| 国产精品一区二区性色av| АⅤ资源中文在线天堂| 久久99热6这里只有精品| 亚洲,欧美,日韩| 亚洲精品日韩av片在线观看| 男女视频在线观看网站免费| 国产精品电影一区二区三区| 亚洲精品成人久久久久久| 最后的刺客免费高清国语| 日韩视频在线欧美| 精品人妻视频免费看| 日日干狠狠操夜夜爽| 日韩大尺度精品在线看网址| 亚洲真实伦在线观看| 亚洲高清免费不卡视频| 亚洲av中文字字幕乱码综合| 深夜a级毛片| 99riav亚洲国产免费| 成人一区二区视频在线观看| 2022亚洲国产成人精品| 亚洲高清免费不卡视频| 日韩在线高清观看一区二区三区| 午夜福利高清视频| 国产一区二区亚洲精品在线观看| 美女高潮的动态| 悠悠久久av| 国产黄色小视频在线观看| 久久久久久久亚洲中文字幕| 亚洲激情五月婷婷啪啪| 99久久无色码亚洲精品果冻| 国产 一区 欧美 日韩| 国产真实伦视频高清在线观看| 久久午夜福利片| 性欧美人与动物交配| 人妻制服诱惑在线中文字幕| 欧美精品国产亚洲| 欧洲精品卡2卡3卡4卡5卡区| 午夜福利高清视频| 日韩强制内射视频| 色尼玛亚洲综合影院| av在线蜜桃| 91在线精品国自产拍蜜月| 一边摸一边抽搐一进一小说| 国产黄色视频一区二区在线观看 | 成人欧美大片| 26uuu在线亚洲综合色| 性色avwww在线观看| 国产伦理片在线播放av一区 | 久久韩国三级中文字幕| 真实男女啪啪啪动态图| av天堂在线播放| 国模一区二区三区四区视频| 天堂av国产一区二区熟女人妻| 婷婷六月久久综合丁香| 乱系列少妇在线播放| 伦精品一区二区三区| 国产精品一及| 麻豆精品久久久久久蜜桃| 女同久久另类99精品国产91| 久久久精品94久久精品| 人妻久久中文字幕网| 国产精品一区二区性色av| 国产精品久久视频播放| 成人欧美大片| www.色视频.com| 熟女人妻精品中文字幕| 黄色日韩在线| 日韩强制内射视频| 亚洲国产精品国产精品| 亚洲乱码一区二区免费版| 在现免费观看毛片| 国产精品乱码一区二三区的特点| 亚洲一级一片aⅴ在线观看| 久99久视频精品免费| 久久久久久久久久久免费av| 国产精品日韩av在线免费观看| 久久精品国产亚洲av天美| 国产三级中文精品| 一本久久精品| 天天一区二区日本电影三级| 欧美性猛交黑人性爽| 老女人水多毛片| 美女大奶头视频| 成人性生交大片免费视频hd| 国产精品.久久久| 嫩草影院精品99| 中文在线观看免费www的网站| 国产精品无大码| 毛片女人毛片| 免费看光身美女| 哪里可以看免费的av片| 久久精品国产鲁丝片午夜精品| 最好的美女福利视频网| 国产精华一区二区三区| 欧美日韩国产亚洲二区| 草草在线视频免费看| av卡一久久| 小蜜桃在线观看免费完整版高清| 只有这里有精品99| 91在线精品国自产拍蜜月| 天堂av国产一区二区熟女人妻| 九九久久精品国产亚洲av麻豆| 国产日本99.免费观看| 日韩欧美 国产精品| 久久午夜福利片| 日日摸夜夜添夜夜添av毛片| 国产精品国产三级国产av玫瑰| 成人二区视频| 久久精品人妻少妇| 性色avwww在线观看| 在线a可以看的网站| 69人妻影院| 黄色欧美视频在线观看| 精品国产三级普通话版| 最近最新中文字幕大全电影3| 亚洲中文字幕日韩| 亚洲七黄色美女视频| 久久久久久久久中文| 亚洲精品粉嫩美女一区| 国产色爽女视频免费观看| 久久精品夜夜夜夜夜久久蜜豆| 少妇熟女欧美另类| 男的添女的下面高潮视频| 国产精品嫩草影院av在线观看| 成人亚洲欧美一区二区av| 在线观看美女被高潮喷水网站| 日韩一区二区视频免费看| 国产成人91sexporn| 99久国产av精品| 国产三级在线视频| 最近2019中文字幕mv第一页| 你懂的网址亚洲精品在线观看 | av又黄又爽大尺度在线免费看 | 国产精品精品国产色婷婷| 亚洲av中文字字幕乱码综合| 97超视频在线观看视频| 亚洲美女搞黄在线观看| 久久欧美精品欧美久久欧美| 禁无遮挡网站| 黄色配什么色好看| 日日干狠狠操夜夜爽| 欧美潮喷喷水| 亚洲自偷自拍三级| 欧美另类亚洲清纯唯美| 久久久午夜欧美精品| 国产精品一区www在线观看| 大型黄色视频在线免费观看| 午夜精品在线福利| 看免费成人av毛片| 一级二级三级毛片免费看| 欧美最黄视频在线播放免费| 色噜噜av男人的天堂激情| 蜜桃久久精品国产亚洲av| 亚洲成人中文字幕在线播放| 婷婷色av中文字幕| 此物有八面人人有两片| 精品一区二区三区视频在线| 国产高清不卡午夜福利| 欧美成人a在线观看| 国产精品一区二区三区四区免费观看| 国产黄a三级三级三级人| 乱码一卡2卡4卡精品| 最新中文字幕久久久久| 久久久国产成人精品二区| 2022亚洲国产成人精品| 99热全是精品| 亚洲精品乱码久久久久久按摩| 国产片特级美女逼逼视频| 中文亚洲av片在线观看爽| 性色avwww在线观看| 国产高清不卡午夜福利| 一区二区三区四区激情视频 | 久久久久久国产a免费观看| 99热这里只有是精品50| 爱豆传媒免费全集在线观看| 精品久久久噜噜| 亚洲av男天堂| 欧洲精品卡2卡3卡4卡5卡区| eeuss影院久久| 欧美丝袜亚洲另类| 久久人人精品亚洲av| 久久精品国产亚洲网站| 此物有八面人人有两片| 国产精品美女特级片免费视频播放器| 日本黄色视频三级网站网址| 一区二区三区免费毛片| 在线观看免费视频日本深夜| 一个人观看的视频www高清免费观看| 亚洲,欧美,日韩| 一级二级三级毛片免费看| 午夜免费男女啪啪视频观看| 国产一区二区三区在线臀色熟女| 97超碰精品成人国产| 国产精品爽爽va在线观看网站| 亚洲丝袜综合中文字幕| 国产av不卡久久| 日韩视频在线欧美| 一边亲一边摸免费视频| 少妇的逼水好多| 国产一区二区激情短视频| 啦啦啦韩国在线观看视频| 色哟哟·www| 尾随美女入室| 亚洲av一区综合|