• <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ì)算方法
    色5月婷婷丁香| 亚洲色图 男人天堂 中文字幕 | 亚洲精品乱码久久久久久按摩| 成人国产麻豆网| 精品一区在线观看国产| 一本一本综合久久| 精品亚洲成a人片在线观看| 欧美老熟妇乱子伦牲交| 最新的欧美精品一区二区| 久久人人爽人人爽人人片va| 中文字幕精品免费在线观看视频 | 国产精品.久久久| 久久久亚洲精品成人影院| 99久久中文字幕三级久久日本| 中文字幕av电影在线播放| 韩国av在线不卡| 大又大粗又爽又黄少妇毛片口| 人体艺术视频欧美日本| 亚洲一级一片aⅴ在线观看| 91久久精品国产一区二区成人| 久久人人爽人人片av| 国产淫语在线视频| 日韩av不卡免费在线播放| 啦啦啦在线观看免费高清www| 国模一区二区三区四区视频| 丁香六月天网| 精品熟女少妇av免费看| 在线观看三级黄色| 人妻人人澡人人爽人人| 男人添女人高潮全过程视频| 伊人亚洲综合成人网| 七月丁香在线播放| 国产国拍精品亚洲av在线观看| 777米奇影视久久| 久久精品久久久久久噜噜老黄| 欧美激情 高清一区二区三区| 亚洲成色77777| 欧美最新免费一区二区三区| 五月开心婷婷网| 亚洲三级黄色毛片| 久久精品国产亚洲av天美| 久久精品国产亚洲网站| 婷婷成人精品国产| 婷婷色av中文字幕| 18禁观看日本| 国产亚洲精品久久久com| 国产老妇伦熟女老妇高清| 水蜜桃什么品种好| 人妻少妇偷人精品九色| 少妇熟女欧美另类| 日韩一区二区三区影片| 国产日韩欧美视频二区| 涩涩av久久男人的天堂| 国产伦理片在线播放av一区| 在线观看免费日韩欧美大片 | 亚洲精品av麻豆狂野| 丝袜在线中文字幕| 波野结衣二区三区在线| 日日撸夜夜添| 国语对白做爰xxxⅹ性视频网站| 午夜精品国产一区二区电影| 久久鲁丝午夜福利片| 制服诱惑二区| av卡一久久| 精品亚洲乱码少妇综合久久| 满18在线观看网站| 国产高清国产精品国产三级| 午夜精品国产一区二区电影| 99精国产麻豆久久婷婷| 国产亚洲欧美精品永久| 久久久久精品性色| 亚洲性久久影院| 国产精品无大码| 国产精品女同一区二区软件| 成年av动漫网址| 成人免费观看视频高清| 桃花免费在线播放| 亚洲欧美一区二区三区国产| 亚洲久久久国产精品| 国产精品无大码| 国产av国产精品国产| 狠狠精品人妻久久久久久综合| 日韩,欧美,国产一区二区三区| 国产一区二区在线观看日韩| 99久久综合免费| 国产乱人偷精品视频| 亚洲成人av在线免费| 少妇人妻 视频| 大片免费播放器 马上看| 黑人巨大精品欧美一区二区蜜桃 | 亚洲国产av影院在线观看| 好男人视频免费观看在线| 久久久久久久久久久丰满| 亚洲成人av在线免费| 最近最新中文字幕免费大全7| 亚洲综合精品二区| 亚洲欧美色中文字幕在线| 亚洲欧美一区二区三区黑人 | 七月丁香在线播放| av女优亚洲男人天堂| 久久这里有精品视频免费| 最新的欧美精品一区二区| 青青草视频在线视频观看| 美女视频免费永久观看网站| 99国产综合亚洲精品| 一本—道久久a久久精品蜜桃钙片| 高清欧美精品videossex| 国产成人精品一,二区| 满18在线观看网站| 亚洲精品久久午夜乱码| 国产亚洲精品第一综合不卡 | √禁漫天堂资源中文www| 久久精品夜色国产| 伦理电影大哥的女人| 99热网站在线观看| 97精品久久久久久久久久精品| 蜜桃在线观看..| 黄片无遮挡物在线观看| 国产免费一区二区三区四区乱码| 久久这里有精品视频免费| 只有这里有精品99| av国产久精品久网站免费入址| 午夜福利视频在线观看免费| 男女边摸边吃奶| 国产日韩一区二区三区精品不卡 | 只有这里有精品99| 最近手机中文字幕大全| 日韩在线高清观看一区二区三区| 亚洲av成人精品一二三区| 中文字幕精品免费在线观看视频 | 婷婷色麻豆天堂久久| 黄片播放在线免费| 老女人水多毛片| 国产成人精品久久久久久| 色哟哟·www| 国产精品国产三级国产专区5o| 亚洲在久久综合| 人人妻人人澡人人看| 国语对白做爰xxxⅹ性视频网站| 精品视频人人做人人爽| 五月伊人婷婷丁香| 美女国产高潮福利片在线看| 另类亚洲欧美激情| 婷婷色综合大香蕉| 亚洲av免费高清在线观看| 亚洲欧美色中文字幕在线| 亚洲精品视频女| 亚洲色图综合在线观看| 欧美一级a爱片免费观看看| 国产精品 国内视频| 日韩成人av中文字幕在线观看| 内地一区二区视频在线| 熟女电影av网| 午夜免费鲁丝| 中文字幕人妻熟人妻熟丝袜美| 满18在线观看网站| 国产成人免费无遮挡视频| av黄色大香蕉| 免费高清在线观看视频在线观看| av电影中文网址| 日韩电影二区| 国产黄频视频在线观看| 2021少妇久久久久久久久久久| 国产成人午夜福利电影在线观看| 七月丁香在线播放| 啦啦啦啦在线视频资源| 国产熟女欧美一区二区| 国产精品欧美亚洲77777| 久久精品久久久久久久性| 免费av不卡在线播放| 中文精品一卡2卡3卡4更新| 少妇熟女欧美另类| 成人亚洲精品一区在线观看| 男女高潮啪啪啪动态图| 国产在线一区二区三区精| 22中文网久久字幕| 毛片一级片免费看久久久久| 精品一区二区三卡| 五月伊人婷婷丁香| 少妇 在线观看| 国产欧美亚洲国产| 国产亚洲欧美精品永久| 亚洲图色成人| 久久免费观看电影| 精品卡一卡二卡四卡免费| 黄色欧美视频在线观看| 最黄视频免费看| 婷婷色综合www| 久久久欧美国产精品| 26uuu在线亚洲综合色| 99久久精品国产国产毛片| 国产成人午夜福利电影在线观看| 桃花免费在线播放| 亚洲国产欧美日韩在线播放| 亚洲av欧美aⅴ国产| 亚洲精品国产色婷婷电影| 女的被弄到高潮叫床怎么办| 国产永久视频网站| 国产免费一区二区三区四区乱码| 在线精品无人区一区二区三| 国产亚洲欧美精品永久| 乱人伦中国视频| 中文天堂在线官网| 国产成人a∨麻豆精品| 男女无遮挡免费网站观看| 午夜av观看不卡| 亚洲精品国产av成人精品| 热99久久久久精品小说推荐| 国产精品国产三级国产av玫瑰| 久久久精品区二区三区| 高清黄色对白视频在线免费看| 久久青草综合色| 亚洲欧美色中文字幕在线| 亚洲色图 男人天堂 中文字幕 | 99热这里只有是精品在线观看| 久久女婷五月综合色啪小说| 欧美性感艳星| 亚洲第一av免费看| 欧美bdsm另类| 国产精品国产三级国产专区5o| 最近中文字幕2019免费版| 91久久精品国产一区二区三区| 久久精品国产鲁丝片午夜精品| 久久久亚洲精品成人影院| 99久久综合免费| 97超碰精品成人国产| 我的女老师完整版在线观看| 精品99又大又爽又粗少妇毛片| 两个人的视频大全免费| 一区二区三区乱码不卡18| 欧美老熟妇乱子伦牲交| 久久久国产一区二区| 欧美 日韩 精品 国产| 男的添女的下面高潮视频| 国产av国产精品国产| 一区二区日韩欧美中文字幕 | 大片电影免费在线观看免费| 久久久久久久大尺度免费视频| 伊人亚洲综合成人网| 亚洲成人手机| 欧美老熟妇乱子伦牲交| 亚洲综合色网址| 美女视频免费永久观看网站| 男女边摸边吃奶| 一本久久精品| 人妻一区二区av| 极品人妻少妇av视频| 韩国av在线不卡| 日本-黄色视频高清免费观看| 天堂8中文在线网| 国产成人精品一,二区| 成人手机av| 99久国产av精品国产电影| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 热99久久久久精品小说推荐| 免费播放大片免费观看视频在线观看| 少妇人妻精品综合一区二区| 国产女主播在线喷水免费视频网站| 在线观看免费高清a一片| 两个人免费观看高清视频| 国产精品秋霞免费鲁丝片| 久久99热6这里只有精品| 狠狠婷婷综合久久久久久88av| 欧美日韩视频高清一区二区三区二| 国产成人精品一,二区| 国产极品粉嫩免费观看在线 | 寂寞人妻少妇视频99o| 卡戴珊不雅视频在线播放| 国产免费又黄又爽又色| 成人国语在线视频| 亚洲精品国产av蜜桃| 伊人亚洲综合成人网| 插阴视频在线观看视频| 看十八女毛片水多多多| 日日爽夜夜爽网站| 欧美人与善性xxx| 777米奇影视久久| av线在线观看网站| 欧美日韩精品成人综合77777| 亚洲av免费高清在线观看| 亚洲经典国产精华液单| 青春草亚洲视频在线观看| 欧美精品高潮呻吟av久久| 少妇被粗大猛烈的视频| 欧美3d第一页| 国产精品国产三级专区第一集| av不卡在线播放| 免费大片黄手机在线观看| av天堂久久9| 成年女人在线观看亚洲视频| 国产欧美另类精品又又久久亚洲欧美| 内地一区二区视频在线| 又大又黄又爽视频免费| 18禁观看日本| 26uuu在线亚洲综合色| 丰满饥渴人妻一区二区三| 久久久久久久久久久久大奶| 一本一本综合久久| 寂寞人妻少妇视频99o| 99九九在线精品视频| 亚洲在久久综合| 搡老乐熟女国产| 亚洲第一av免费看| xxx大片免费视频| 国产日韩欧美亚洲二区| 国产伦精品一区二区三区视频9| 嘟嘟电影网在线观看| 亚洲激情五月婷婷啪啪| 简卡轻食公司| 成人毛片60女人毛片免费| 99九九在线精品视频| 免费大片黄手机在线观看| 黑人欧美特级aaaaaa片| 午夜免费观看性视频| av又黄又爽大尺度在线免费看| 九九爱精品视频在线观看| 哪个播放器可以免费观看大片| 少妇人妻精品综合一区二区| 成人毛片60女人毛片免费| 美女cb高潮喷水在线观看| av有码第一页| av卡一久久| 大香蕉久久网| 热re99久久国产66热| 黄色配什么色好看| 国产黄片视频在线免费观看| 99re6热这里在线精品视频| 久久久久人妻精品一区果冻| 热99久久久久精品小说推荐| 丰满乱子伦码专区| 日韩精品免费视频一区二区三区 | 激情五月婷婷亚洲| 成人午夜精彩视频在线观看| 精品一区在线观看国产| 天天影视国产精品| 久久久国产一区二区| 国产日韩欧美在线精品| 一级片'在线观看视频| 丁香六月天网| 国产欧美另类精品又又久久亚洲欧美| 免费少妇av软件| 七月丁香在线播放| 在线观看人妻少妇| 色婷婷av一区二区三区视频| 久久久国产欧美日韩av| 最近最新中文字幕免费大全7| 乱人伦中国视频| 搡女人真爽免费视频火全软件| 在线观看www视频免费| 欧美精品一区二区大全| 久久久久精品性色| 亚洲精品自拍成人| 国产亚洲一区二区精品| 中文欧美无线码| 国产日韩欧美视频二区| 国产精品无大码| 国产视频内射| 超色免费av| 日韩一区二区三区影片| 如日韩欧美国产精品一区二区三区 | av电影中文网址| 欧美日韩综合久久久久久| 久久韩国三级中文字幕| av有码第一页| 男女高潮啪啪啪动态图| 日日摸夜夜添夜夜添av毛片| 久久综合国产亚洲精品| 在线精品无人区一区二区三| 人人澡人人妻人| 国产精品不卡视频一区二区| 精品一区在线观看国产| 午夜福利,免费看| 在线观看国产h片| 国产色婷婷99| 少妇人妻 视频| 丝瓜视频免费看黄片| 亚洲怡红院男人天堂| av线在线观看网站| 亚洲精品乱久久久久久| 久久精品熟女亚洲av麻豆精品| 一本—道久久a久久精品蜜桃钙片| 国产在线一区二区三区精| 少妇人妻久久综合中文| 久久久久久久亚洲中文字幕| 中文精品一卡2卡3卡4更新| 亚洲三级黄色毛片| 精品久久国产蜜桃| 少妇的逼水好多| 久久久久人妻精品一区果冻| 精品人妻熟女毛片av久久网站| 黄片播放在线免费| 亚洲中文av在线| 日日摸夜夜添夜夜爱| 女性生殖器流出的白浆| xxxhd国产人妻xxx| 18+在线观看网站| 亚洲av.av天堂| 亚洲国产精品成人久久小说| 亚洲欧美清纯卡通| av有码第一页| 国产视频首页在线观看| 午夜免费观看性视频| 亚洲av.av天堂| 亚洲人成网站在线播| 街头女战士在线观看网站| a级毛片黄视频| 国产成人精品在线电影| 欧美另类一区| 亚洲第一av免费看| 久久久久久久久大av| 18禁观看日本| 午夜av观看不卡| 免费观看的影片在线观看| 国产高清不卡午夜福利| 丰满乱子伦码专区| av专区在线播放| 久久久久久伊人网av| 高清毛片免费看| 免费观看a级毛片全部| 最近2019中文字幕mv第一页| 精品熟女少妇av免费看| 亚洲成人手机| 国产成人freesex在线| 在线观看美女被高潮喷水网站| 国产精品不卡视频一区二区| 日韩强制内射视频| 日韩电影二区| 久久久久久久久大av| 狠狠精品人妻久久久久久综合| 超色免费av| 少妇人妻精品综合一区二区| 国产成人精品无人区| 婷婷色综合www| 色94色欧美一区二区| 国精品久久久久久国模美| 18+在线观看网站| 日韩av免费高清视频| 成人综合一区亚洲| 日韩制服骚丝袜av| 亚洲精品456在线播放app| 精品人妻偷拍中文字幕| 欧美日韩国产mv在线观看视频| 国产精品.久久久| 日本av手机在线免费观看| 精品人妻在线不人妻| av又黄又爽大尺度在线免费看| 亚洲国产欧美日韩在线播放| 亚洲av国产av综合av卡| 老司机影院成人| 国产精品一区二区在线观看99| 亚洲精品国产色婷婷电影| 久久婷婷青草| 少妇人妻 视频| 在线观看www视频免费| 两个人的视频大全免费| 丰满饥渴人妻一区二区三| 最新中文字幕久久久久| 国产日韩欧美在线精品| 成人漫画全彩无遮挡| 少妇高潮的动态图| 国产 一区精品| 亚洲精品456在线播放app| 九九久久精品国产亚洲av麻豆| 午夜福利视频在线观看免费| a级毛片免费高清观看在线播放| 尾随美女入室| 精品卡一卡二卡四卡免费| 满18在线观看网站| 欧美丝袜亚洲另类| 国产成人91sexporn| 免费少妇av软件| 久久ye,这里只有精品| 精品少妇黑人巨大在线播放| 王馨瑶露胸无遮挡在线观看| 满18在线观看网站| 超碰97精品在线观看| 亚洲伊人久久精品综合| 日韩精品免费视频一区二区三区 | 欧美国产精品一级二级三级| 亚洲国产日韩一区二区| 国产精品久久久久成人av| 亚洲av中文av极速乱| 三上悠亚av全集在线观看| 亚州av有码| 久久久精品免费免费高清| 国产综合精华液| 色网站视频免费| 国产精品99久久99久久久不卡 | 夜夜爽夜夜爽视频| 中国国产av一级| 国产精品一区二区三区四区免费观看| 热99国产精品久久久久久7| 亚洲精品亚洲一区二区| 自线自在国产av| 久久人人爽人人爽人人片va| 不卡视频在线观看欧美| 国产精品久久久久成人av| 亚洲精品久久久久久婷婷小说| 成年av动漫网址| 久久人人爽人人爽人人片va| 国产探花极品一区二区| 大香蕉97超碰在线| 国产精品无大码| 日本91视频免费播放| 多毛熟女@视频| 国产成人精品福利久久| 亚洲av免费高清在线观看| 久久99蜜桃精品久久| 精品人妻一区二区三区麻豆| 国产精品三级大全| 久久久久久久久久久免费av| 精品人妻熟女av久视频| 久久久久精品久久久久真实原创| 一级毛片aaaaaa免费看小| av在线播放精品| 三级国产精品欧美在线观看| 一级黄片播放器| 看非洲黑人一级黄片| 大又大粗又爽又黄少妇毛片口| 校园人妻丝袜中文字幕| 国产免费福利视频在线观看| 日本黄色片子视频| 免费不卡的大黄色大毛片视频在线观看| 亚洲欧美成人精品一区二区| 成人无遮挡网站| 久久久亚洲精品成人影院| 熟妇人妻不卡中文字幕| 97在线视频观看| videosex国产| 久久精品夜色国产| 亚洲伊人久久精品综合| 国产精品不卡视频一区二区| 午夜福利影视在线免费观看| 久久99蜜桃精品久久| 久久久a久久爽久久v久久| 美女xxoo啪啪120秒动态图| 国产有黄有色有爽视频| 人妻系列 视频| 国产一区二区在线观看av| 热99久久久久精品小说推荐| 日韩欧美精品免费久久| 91国产中文字幕| 少妇精品久久久久久久| 久久久久久久精品精品| 人妻夜夜爽99麻豆av| 国产精品蜜桃在线观看| 午夜久久久在线观看| 久久久久国产网址| 国产精品99久久99久久久不卡 | 两个人免费观看高清视频| 毛片一级片免费看久久久久| 视频在线观看一区二区三区| 极品少妇高潮喷水抽搐| 国产精品女同一区二区软件| 欧美xxⅹ黑人| 亚洲国产成人一精品久久久| 男女边吃奶边做爰视频| 免费高清在线观看日韩| 日韩一区二区三区影片| 亚州av有码| 亚洲精品美女久久av网站| 欧美bdsm另类| 99久久精品一区二区三区| 久久 成人 亚洲| 男女边摸边吃奶| 极品少妇高潮喷水抽搐| 欧美日韩视频高清一区二区三区二| 寂寞人妻少妇视频99o| 久久久久久久大尺度免费视频| 国产亚洲午夜精品一区二区久久| .国产精品久久| 亚洲精品第二区| 中国三级夫妇交换| 欧美+日韩+精品| 国内精品宾馆在线| 哪个播放器可以免费观看大片| 日韩电影二区| av国产久精品久网站免费入址| 国产精品偷伦视频观看了| 国产日韩欧美在线精品| 91精品国产国语对白视频| 天天躁夜夜躁狠狠久久av| 老熟女久久久| 黑人欧美特级aaaaaa片| 两个人的视频大全免费| 91精品国产九色| 一边亲一边摸免费视频| 夫妻性生交免费视频一级片| 欧美丝袜亚洲另类| 国产男女内射视频| 免费高清在线观看日韩| 国产69精品久久久久777片| 黑人高潮一二区| 伦理电影大哥的女人| 少妇熟女欧美另类| 国产精品人妻久久久影院| 伦理电影大哥的女人| 国产乱来视频区| 中文精品一卡2卡3卡4更新| 国产免费现黄频在线看| 一级黄片播放器| 亚洲国产欧美在线一区| 亚洲成色77777| 免费看av在线观看网站| 国产精品三级大全| 熟女电影av网| 纯流量卡能插随身wifi吗| 国产精品一区二区在线不卡| 欧美三级亚洲精品| 中文字幕制服av| 欧美国产精品一级二级三级| 全区人妻精品视频| 18禁观看日本| 黑人欧美特级aaaaaa片| 天天操日日干夜夜撸|