亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        基于DS數(shù)據(jù)融合與SVR-PF的鋰離子電池RUL預(yù)測方法

        2015-04-21 01:27:20王常虹董漢成凌明祥李清華
        東北石油大學(xué)學(xué)報 2015年4期
        關(guān)鍵詞:充放電鋰離子容量

        王常虹,董漢成,凌明祥,李清華

        ?

        基于DS數(shù)據(jù)融合與SVR-PF的鋰離子電池RUL預(yù)測方法

        王常虹,董漢成,凌明祥,李清華

        (哈爾濱工業(yè)大學(xué) 空間控制與慣性技術(shù)研究中心,黑龍江 哈爾濱 150001 )

        為防止鋰離子電池失效導(dǎo)致的系統(tǒng)失效,提出一種基于DS數(shù)據(jù)融合與支持向量回歸機粒子濾波(Support Vector Regression-Particle Filter,SVR-PF)的鋰離子電池剩余有效工作時間(Remaining Useful Life,RUL)預(yù)測方法.結(jié)果表明:該預(yù)測方法能夠融合不同數(shù)據(jù)源對鋰離子電池RUL的預(yù)測結(jié)果,改進可用數(shù)據(jù)較少時RUL的預(yù)測準(zhǔn)確度.

        鋰離子電池; RUL; DS數(shù)據(jù)融合; SVR-PF

        0 引言

        鋰離子電池被廣泛應(yīng)用于各種電子設(shè)備,為設(shè)備的正常運行提供必要的能源[1-2].電池的剩余有效工作時間(Remaining Useful Life,RUL)定義為在一個特定的運行條件下,電池從當(dāng)前運行狀態(tài)衰退至無法正常工作狀態(tài)所需要的RUL,為電池視情維護(Condition Based Maintenance,CBM)的重要指標(biāo)[3].因此,尋找一個可靠且準(zhǔn)確的途徑預(yù)測鋰離子電池的RUL,可以提供鋰離子電池系統(tǒng)維護或更換的相關(guān)信息[4-5].

        數(shù)據(jù)驅(qū)動算法(Data-driven Algorithm)可以利用實測電池數(shù)據(jù)建立電池的回歸模型,用于鋰離子電池RUL預(yù)測[6].Juan Carlos Alvarez Anton等[7]應(yīng)用支持向量回歸機(Support Vector Regression,SVR),分析鋰離子電池荷電狀態(tài)與電池電壓、電流、溫度的關(guān)系,建立估計電池使用狀態(tài)的數(shù)據(jù)驅(qū)動模型,分析電池的健康狀態(tài)(State-of-Health,SoH),預(yù)測電池的RUL.Weng Caihao[8]與Adnan Nuhic[9]等利用支持向量回歸機,建立電池的健康狀態(tài)(State of Health,SoH)模型,即電池容量與充放電周期、電池開路電壓的關(guān)系,預(yù)測電池性能衰退與RUL.遞歸濾波算法(Recursive Filter Algorithm)基于電池狀態(tài)空間模型,可以利用實測電池數(shù)據(jù)分析、估計與預(yù)測電池的故障,用于鋰離子電池RUL預(yù)測[10].Dare Andre等[11]建立含有無跡卡爾曼濾波器(Unscented Kalman Filter,UKF)與卡爾曼濾波器的雙濾波器估計模型,利用鋰離子電池的電壓、電流與充放電周期數(shù)估計電池的剩余容量與荷電狀態(tài),預(yù)測鋰離子電池的RUL.He Wei[12]與Miao Qiang[13]等將衛(wèi)星鋰離子電池的容量作為分析參數(shù),對容量隨時間衰退的特性進行數(shù)據(jù)挖掘,建立衰退趨勢方程,以估計容量隨時間的衰退趨勢,估計鋰離子電池的當(dāng)前工作狀態(tài),預(yù)測鋰離子電池的RUL.Bhaskar Saha等[14]應(yīng)用向量機平滑后得到的阻抗數(shù)據(jù)與容量數(shù)據(jù),分析衛(wèi)星鋰離子電池性能隨時間衰退的趨勢,預(yù)測鋰離子電池的RUL.

        鋰離子電池的RUL預(yù)測需要分析可用數(shù)據(jù),以估計電池的當(dāng)前狀態(tài),預(yù)測電池性能的衰退趨勢,在可用數(shù)據(jù)量很少時,鋰離子電池RUL的預(yù)測準(zhǔn)確度降低,現(xiàn)有文獻未提出有效解決方法.筆者采用基于阻抗的鋰離子電池RUL預(yù)測方法[15]、基于容量的鋰離子電池RUL預(yù)測方法[16],以及提出的基于DS數(shù)據(jù)融合理論(Dempster-Shafer Theory)與支持向量回歸機粒子濾波(Support Vector Regression-Particle Filter,SVR-PF[15])的預(yù)測方法,融合不同數(shù)據(jù)源對鋰離子電池RUL分析結(jié)果,比較3種方法的預(yù)測性能,以在可用數(shù)據(jù)較少時得到準(zhǔn)確的預(yù)測結(jié)果.

        1 DS數(shù)據(jù)融合

        1.1 理論

        DS數(shù)據(jù)融合理論(Dempster-Shafer Theory)[17]利用基本概率賦值 (Basic Probability Assignment,BPA) 函數(shù)對來自不同源的數(shù)據(jù)進行融合,對辨識框架(Frame of Discernment)中所有可能發(fā)生的命題進行可信度(Belief)分析,利用融合規(guī)則結(jié)合所有命題發(fā)生的可信度,從而達到數(shù)據(jù)融合的目的.

        設(shè)m1與m2為在辨識框架中兩個BPA函數(shù),m1,2為融合后的BPA函數(shù),則DS數(shù)據(jù)融合可以表示為

        m1,2(?)=0,

        (1)

        (2)

        其中

        (3)

        DS數(shù)據(jù)融合理論結(jié)合對同一問題不同源的相同觀點,同時剔除所有的沖突觀點,可以得到更加可靠的融合后驗BPA函數(shù).

        1.2 應(yīng)用

        基于DS數(shù)據(jù)融合的鋰離子電池RUL預(yù)測具有兩個數(shù)據(jù)源:(1)利用阻抗數(shù)據(jù)的RUL預(yù)測;(2)利用容量數(shù)據(jù)的RUL預(yù)測.筆者提出的基于DS數(shù)據(jù)融合與SVR-PF的鋰離子電池RUL預(yù)測融合兩種方法得到的預(yù)測結(jié)果,最終得到融合預(yù)測結(jié)果.

        將全集即辨識框架Ω記為

        Ω={im,cap},

        (4)

        則冪集2Ω為

        2Ω={?,{im},{cap},{im∪cap}},

        其中冪集2Ω中所有命題的含義為

        (1){im}代表利用阻抗數(shù)據(jù)預(yù)測得到的RUL可信度;

        (2){cap}代表利用容量數(shù)據(jù)預(yù)測得到的RUL可信度;

        (3){im∪cap}代表利用阻抗或容量數(shù)據(jù)預(yù)測得到的RUL可信度.

        同時,定義在冪集2Ω上的BPA函數(shù)m1與m2代表的含義為

        (1)m1為利用阻抗數(shù)據(jù)對RUL進行預(yù)測時對冪集2Ω中命題的可信度分配;

        (2)m2為利用容量數(shù)據(jù)對RUL進行預(yù)測時對冪集2Ω中命題的可信度分配.

        基于DS數(shù)據(jù)融合的BPA函數(shù)結(jié)合方式見表 1,其中x=m1(im∪cap)m2(im∪cap).

        表1 基于DS數(shù)據(jù)融合的BPA函數(shù)結(jié)合方式Table 1 The combination of BPA function based on DS theory

        基于表 1,可以得到數(shù)據(jù)融合后的后驗BPA函數(shù)m為

        (6)

        (7)

        其中

        a=m1(im)m2(im)+m1(im∪cap)m2(im)+m1(im)m2(im∪cap),

        (8)

        b=m1(cap)m2(cap)+m1(im∪cap)m2(cap)+m1(cap)m2(im∪cap),

        (9)

        (10)

        由此得到DS數(shù)據(jù)融合在鋰離子電池RUL預(yù)測中方法,利用后驗融合BPA函數(shù)式與將兩類方法的預(yù)測結(jié)果結(jié)合,以得到更加準(zhǔn)確的預(yù)測結(jié)果.

        2 實測鋰離子電池阻抗與容量數(shù)據(jù)

        所用鋰離子電池阻抗與容量數(shù)據(jù)在美國愛達荷國家實驗室(IdahoNationalLaboratory)測得[18].實驗用鋰離子電池型號為Gen2 18650-size,實驗在室溫(24 ℃)與43 ℃溫度條件下恒溫進行,分為3個循環(huán)的測試階段:(1)充電:1.5A恒流充電至電壓4.2V,隨后恒壓充電至電流20mA;(2)放電:2.0A電流恒壓放電至電壓2.5V;(3)阻抗:進行0.1~5.0kHz掃頻,得到電化學(xué)阻抗譜,進而得到電池阻抗.

        共用4塊電池.電池5、6與7的實測數(shù)據(jù)見圖1,在室溫(24 ℃)條件下恒溫測得;電池32的實測數(shù)據(jù)見圖2,在43 ℃溫度條件下恒溫測得.由圖1-2可以看出,鋰離子電池容量隨著充放電周期的增加而減少,阻抗隨著充放電周期的增加而增加.利用測得的電池數(shù)據(jù)驗證文中提出的鋰離子電池RUL預(yù)測方法的性能.

        圖1 電池5、6、7容量與阻抗數(shù)據(jù)Fig.1 Capacity and impedance data of battery 5,6 and 7

        圖2 電池32容量與阻抗數(shù)據(jù)Fig.2 Capacity and impedance data of battery 32

        持續(xù)的充放電過程加速電池的性能衰退,同時阻抗測試提供電池內(nèi)部參數(shù)的變化方式.當(dāng)電池達到壽命終點時實驗停止,將電池達到壽命終點(EndofLife,EOL)時容量占標(biāo)稱容量的百分比稱為EOL門限.由于各個電池停止實驗時的容量值不同,因此EOL門限也不同,電池5、6、7、32的EOL門限分別為70%、65%、75%、95%.

        3 基于DS數(shù)據(jù)融合與SVR-PF的鋰離子電池RUL預(yù)測方法

        3.1 預(yù)測過程

        3.1.1 初始狀態(tài)

        (1)首先分析BPA函數(shù)m1,i(im)的初值,由中心極限定理可知,大量的阻抗數(shù)據(jù)測量誤差服從正態(tài)分布,即

        (11)

        (12)

        由文獻[15]可知:

        (13)

        因此容量估計Cim,N服從正態(tài)分布,即

        (14)

        其中

        μCim,N=αN(μRe,N+μRct,N)+βN=CN,

        (15)

        (16)

        由此得到BPA函數(shù)m1,i(im)的初值m1,N+1(im)為

        (17)

        (2)然后計算BPA函數(shù)m2,i(cap),同樣由中心極限定理可知,大量的容量數(shù)據(jù)測量誤差服從正態(tài)分布,即

        (18)

        μCcap,N=CN,

        (19)

        則BPA函數(shù)m2,i(cap)的初值m2,N+1(cap)為

        (20)

        在得到BPA函數(shù)m1,N+1(im)與m2,N+1(cap)值后,對其他的BPA函數(shù)值進行分析.由于兩類預(yù)測方法沒有關(guān)聯(lián),可得

        m1,N+1(cap)=m2,N+1(im)=0.

        (21)

        根據(jù)BPA函數(shù)的性質(zhì)可得

        m1,N+1(im∪cap)=1-m1,N+1(im)-m1,N+1(cap)=1-m1,N+1(im),

        (22)

        m2,N+1(im∪cap)=1-m2,N+1(im)-m2,N+1(cap)=1-m2,N+1(cap).

        (23)

        (3)在確定所有的BPA函數(shù)值后,計算后驗融合BPA函數(shù)值,將融合后的BPA函數(shù)表示為m,由式(6)與式(7)可知:

        (24)

        其中

        (25)

        其中

        m1,N+1(im)m2,N+1(cap)=m1,N+1(im)m2,N+1(cap).

        (26)

        (27)

        由此得到基于DS數(shù)據(jù)融合與SVR-PF鋰離子電池RUL預(yù)測的初始狀態(tài),即

        Initial=[m1,N+1(im),m1,N+1(cap),m1,N+1(im∪cap),m2,N+1(im),m2,N+1(cap),

        (28)

        3.1.2 預(yù)測過程

        (29)

        (30)

        由于兩類預(yù)測模型沒有關(guān)聯(lián),可知

        m1,k+1(cap)=m2,k+1(im)=0.

        (31)

        基于BPA函數(shù)性質(zhì),有

        m1,k+1(im∪cap)=1-m1,k+1(im),

        (32)

        m2,k+1(im∪cap)=1-m2,k+1(cap).

        (33)

        步驟2:預(yù)測第k+1個充放電周期的容量值.由式(6-7)、式(29-33)可知,后驗融合BPA函數(shù)mk+1(im)與mk+1(cap)為

        (34)

        (35)

        (36)

        (37)

        3.2 預(yù)測模型

        建立基于DS數(shù)據(jù)融合與SVR-PF的鋰離子電池RUL預(yù)測模型:

        (38)

        其中狀態(tài)Xk為

        (39)

        式中:Xim,k為利用阻抗數(shù)據(jù)分析得到的狀態(tài);Xcap,k為利用容量數(shù)據(jù)分析得到的狀態(tài);XDS,k為利用DS數(shù)據(jù)融合分析得到的狀態(tài),并且有

        (40)

        (41)

        (42)

        式(40-42)中:λR,k為第k個充放電周期辨識出的阻抗衰退參數(shù)[15];R1,k、λR1,k、R2,k和λR2,k為在第k個充放電周期的基于容量數(shù)據(jù)分析得到的容量衰退參數(shù)[16];*為對應(yīng)變量的預(yù)測值.

        由此,結(jié)合文獻[15-16]RUL預(yù)測模型,可得預(yù)測模型式(38)的狀態(tài)方程,其中利用阻抗數(shù)據(jù)分析的部分狀態(tài)方程為

        (43)

        式中:αN、βN為第N個充放電周期的基于阻抗數(shù)據(jù)分析得到的容量衰退參數(shù)估計值.

        利用容量數(shù)據(jù)分析的部分狀態(tài)方程:

        (44)

        利用DS數(shù)據(jù)融合的部分狀態(tài)方程:

        (45)

        組合式(43-45)得到預(yù)測模型式(38)的狀態(tài)方程,式(38)的測量方程表示為

        (46)

        圖3 基于DS數(shù)據(jù)融合與SVR-PF的鋰離子電池RUL預(yù)測方法流程Fig.3 Process of Lithium-ion battery RUL prediction based on DS theory and SVR-PF

        綜上所述,基于DS數(shù)據(jù)融合與SVR-PF的鋰離子電池RUL預(yù)測方法流程見圖3.

        4 仿真實例

        為了對比提出的基于DS數(shù)據(jù)融合與SVR-PF的鋰離子電池RUL預(yù)測方法,與文獻[15]利用阻抗的RUL預(yù)測方法和文獻[16]利用容量的RUL預(yù)測方法在可用數(shù)據(jù)較少時的RUL預(yù)測性能,應(yīng)用電池5、6、7、32的數(shù)據(jù).將電池5、6、7、32的RUL預(yù)測門限分別設(shè)為90%、85%、95%、95%標(biāo)稱容量.RUL預(yù)測門限值定義為低于電池最大容量(100%)的百分比,其對應(yīng)容量高于電池EOL門限對應(yīng)容量.高RUL預(yù)測門限使得算法用于分析電池健康狀態(tài)的數(shù)據(jù)變少,使用條件變得苛刻,以檢驗算法在可用數(shù)據(jù)較少時的性能.

        電池RUL與EOL的預(yù)測結(jié)果見表2.電池RUL與EOL的誤差預(yù)測結(jié)果見表3.3種方法的電池5、6、7、32的數(shù)據(jù)預(yù)測結(jié)果見圖4-7.

        由表2-3和圖4-7可知,文中提出的基于DS數(shù)據(jù)融合與SVR-PF的RUL預(yù)測方法比文獻[15]與文獻[16]方法的預(yù)測結(jié)果更準(zhǔn)確,具有更小的預(yù)測誤差.原因是提出的基于DS數(shù)據(jù)融合與SVR-PF的RUL預(yù)測方法可以融合另外2種方法的預(yù)測結(jié)果,相對于單一方法的預(yù)測結(jié)果更加可靠.

        表2 3種方法電池的RUL與EOL預(yù)測結(jié)果 充放電周期Table 2 RUL and EOL prediction results of the 3 methods %

        表3 3種方法電池的RUL與EOL預(yù)測誤差結(jié)果 充放電周期Table 3 RUL and EOL prediction errors of the 3 methods %

        圖4 電池5的3種方法預(yù)測數(shù)據(jù)Fig.4 Battery 5 prediction results of the 3 methods

        圖5 電池6的3種方法預(yù)測數(shù)據(jù)Fig.5 Battery 6 prediction results of the 3 methods

        圖6 電池7的3種方法預(yù)測數(shù)據(jù)Fig.6 Battery 7 prediction results of the 3 methods

        圖7 電池32的3種方法預(yù)測數(shù)據(jù)Fig.7 Battery 7 prediction results of the 3 methods

        5 結(jié)束語

        提出一種基于DS數(shù)據(jù)融合與SVR-PF的鋰離子電池RUL預(yù)測方法,它有效結(jié)合利用阻抗數(shù)據(jù)與利用容量數(shù)據(jù)的RUL預(yù)測方法,使得RUL預(yù)測過程更加可靠.在可用數(shù)據(jù)較少時,與利用阻抗和利用容量的RUL方法相比鋰離子電池RUL預(yù)測結(jié)果更加準(zhǔn)確.

        [1] Lu Languang, Han Xuebing, Li Jianqiu, et al.A review on the key issues for lithium-ion battery management in electric vehicles [J].Journal of Power Sources, 2013,226(1):272-288.

        [2] Xing Yinjiao, Eden W M, Tsui K L, et al.Battery management systems in electric and hybrid vehicles [J].Energies, 2011,4(12):1840-1857.

        [3] Si Xiaosheng, Wang Wenbin, Hu Changhua, et al.Remaining useful life estimation-a review on the statistical data driven approaches [J].European Journal of Operational Research, 2011,213(1):1-14.

        [4] Wang Dong, Miao Qiang, Michael P.Prognostics of lithium-ion batteries based on relevance vectors and a conditional three-parameter capacity degradation model [J].Journal of Power Sources, 2013,239(1):253-264.

        [5] Eddahech A, Briat O, Woirgard E, et al.Remaining useful life prediction of lithium batteries in calendar ageing for automotive applications [J].Microelectronics Reliability, 2012,52(1):2438-2442.

        [6] Seyed Mohammad Rezvanizaniani, Liu Zongchang, Chen Yan, et al.Review and recent advances in battery health monitoring and prognostics technologies for electric vehicle (Ev) safety and mobility [J].Journal of Power Sources, 2014,256(1):110-124.

        [7] Juan Carlos Alvarez Anton, Paulino Jose Garca Nieto, Cecilio Blanco Viejo, et al.Support vector machines used to estimate the battery state of charge [J].IEEE Transactions on Power Electronics, 2013,28(12):5919-5926.

        [8] Weng Caihao, Sun Jing, Peng Huei.A unified open-circuit-voltage model of Lithium-Ion Batteries for State-of-Charge estimation and State-of-Health monitoring [J].Journal of Power Sources, 2014,258(1):228-237.

        [9] Adnan Nuhic, Tarik Terzimehic, Thomas Soczka-Guth, et al.Health diagnosis and remaining useful Life prognostics of lithium-ion batteries using data-driven methods [J].Journal of Power Sources, 2013,239(1):680-688.

        [10] Gregory L P.Extended Kalman filtering for battery management systems of lipb-based hev battery packs [J].Journal of Power Sources, 2004,134(2):252-292.

        [11] Dave Andre, Christian Appel, Thomas Soczka-Guth, et al.Advanced mathematical methods of SOC and SOH estimation for lithium-ion batteries [J].Journal of Power Sources, 2013,224(1):20-27.

        [12] He Wei, Nicholas Williard, Michael Osterman, et al.Prognostics of lithium-ion batteries based on Dempster-Shafer theory and the bayesian Monte Carlo method [J].Journal of Power Sources, 2011,196(1):10314-10321.

        [13] Miao Qiang, Xie Lei, Cui Hengjuan, et al.Remaining useful life prediction of lithium-ion battery with unscented particle filter technique [J].Microelectronics Reliability, 2013,53(1):805-810.

        [14] Bhaskar Saha, Kai Goebel, Scott Poll, et al.Prognostics methods for battery health monitoring using a bayesian framework [J].IEEE Transactions on Instrumentation and Measurement, 2009,58(2):291-296.

        [15] Dong Hancheng, Jin Xiaoning, Lou Yangbing, et al.Lithium-ion battery state of health monitoring and remaining useful life prediction based on support vector Regression-Particle filter [J].Journal of Power Sources, 2014,2014(1):114-123.

        [16] 王常虹,董漢成,凌明祥,等.車用鋰離子電池剩余有效工作時間預(yù)測方法[J].汽車工程,2015,37(4):476-479.Wang Changhong, Dong Hancheng, Ling Mingxiang, et al.Remaining useful life prediction of automotive lithium-ion battery [J].Automotive Engineering, 2015,37(4):476-479.

        [17] Priyanka Aggarwal, Deepak Bhatt, Vijay Devabhaktuni.Dempster shafer neural network algorithm for land vehicle navigation application [J].Information Sciences, 2013,253:26-33.

        [18] Kai Goebel Bhaskar Saha, Jon Christophersen.Comparison of prognostic algorithms for estimating remaining useful life of batteries [J].Transactions of the Institute of Measurement and Control, 2009,31(3):293-308.

        2015-04-27;編輯:任志平

        國家自然科學(xué)基金項目(61375046);中央高?;究蒲袠I(yè)務(wù)費專項基金項目(HIT.NSRIF 2014031)

        王常虹(1961-),男,博士,教授,主要從事智能控制與智能系統(tǒng)方面的研究.

        TM912

        A

        2095-4107(2015)04-0109-10

        DOI 10.3969/j.issn.2095-4107.2015.04.014

        猜你喜歡
        充放電鋰離子容量
        V2G模式下電動汽車充放電效率的研究
        高能鋰離子電池的“前世”與“今生”
        科學(xué)(2020年1期)2020-08-24 08:07:56
        基于SG3525的電池充放電管理的雙向DC-DC轉(zhuǎn)換器設(shè)計
        電子制作(2019年23期)2019-02-23 13:21:36
        鋰離子電池充放電保護電路的研究
        SnO2納米片容量異常行為的新解釋
        鋰離子動力電池的不同充電方式
        2015年上半年我國風(fēng)電新增并網(wǎng)容量916萬千瓦
        風(fēng)能(2015年8期)2015-02-27 10:15:12
        2015年一季度我國風(fēng)電新增并網(wǎng)容量470萬千瓦
        風(fēng)能(2015年5期)2015-02-27 10:14:46
        V2G充放電機的設(shè)計及其仿真
        鋰離子電池組不一致性及其彌補措施
        汽車電器(2014年5期)2014-02-28 12:14:15
        老熟妇高潮喷了╳╳╳| 亚洲一区二区三区在线看| 亚洲一区第二区三区四区| 久久99亚洲精品久久久久| 少妇被躁到高潮和人狍大战| 一区二区三区视频免费观看在线 | 国产精品女主播福利在线| 亚洲av一二三四区四色婷婷| 99精品热这里只有精品| 国产精品自产拍在线观看中文| 中文在线最新版天堂av| 日本高清视频在线观看一区二区| 成人无码av免费网站| 九一九色国产| 亚洲最新版无码AV| 国产成人福利av一区二区三区| 日韩av在线亚洲女同| 亚洲av中文无码乱人伦在线视色| 国产国拍亚洲精品mv在线观看| 无码免费人妻超级碰碰碰碰| 亚洲av影片一区二区三区| 日韩av一区二区三区精品久久| 体验区试看120秒啪啪免费| 亚洲免费网站观看视频| 国产精品丝袜在线不卡| 国产青青草自拍视频在线播放| 国产精品成人自拍在线观看| 亚洲av无码乱码国产一区二区| 亚洲综合色自拍一区| www.亚洲天堂.com| 亚洲一区二区日韩精品| 狠狠躁夜夜躁av网站中文字幕| 亚洲精品天堂成人片av在线播放| 无码国产精品一区二区免费16| 青青草一级视频在线观看| 在线国人免费视频播放| 无码人妻精品一区二区三| 男女肉粗暴进来120秒动态图| 宅宅午夜无码一区二区三区| 丝袜美腿一区在线观看| 日韩人妻精品视频一区二区三区|