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

    鋰電池異常狀態(tài)監(jiān)測的FEKF方法研究

    2020-12-30 08:45:50朱浩張瀟
    關(guān)鍵詞:狀態(tài)參數(shù)端電壓充放電

    朱浩,張瀟

    (湖南大學(xué)機(jī)械與運載工程學(xué)院,湖南長沙 410082)

    得益于高能量密度、較高的電壓平臺以及較長的循環(huán)壽命[1],鋰離子電池已被廣泛應(yīng)用于各領(lǐng)域的能源系統(tǒng)中,但鋰離子電池在使用過程中有可能出現(xiàn)過充、過放、溫度過高、溫度過低、電池老化等異常狀態(tài)[2],并由此引起電池?fù)p壞,甚至可能導(dǎo)致燃燒爆炸等安全性事故.因此,對鋰離子電池工作時的異常狀態(tài)進(jìn)行監(jiān)測是十分必要的.

    Singh 等提出了使用卡爾曼濾波(Kalman Filter,KF)生成電池端電壓殘差,然后通過端電壓殘差的零點漂移的大小判斷電池是否有故障的方法[3-4],其中文獻(xiàn)[3]使用一組不同的擴(kuò)展卡爾曼濾波(Extended Kalman Filter,EKF)作為殘差生成器,但EKF 遇到緩變型狀態(tài)參數(shù)異常時可能在故障發(fā)生早期無法辨識甚至發(fā)散;文獻(xiàn)[4]使用了改進(jìn)的自適應(yīng)無跡卡爾曼濾波(Adaptive Unscented Kalman Filter,AUKF)生成端電壓殘差,AUKF 能夠更有效地辨識電池內(nèi)阻R0與時間常數(shù)τ 的緩變故障.Chen 等通過遞推最小二乘法(Recursive Least Square,RLS)與EKF 計算電池內(nèi)阻R0、極化電阻Rp、極化電容Cp,然后得到這些參數(shù)的局部離群因子與故障診斷閾值,判斷電池系統(tǒng)某單體電池電壓的瞬時變化與持續(xù)異常變化,但如果遇到電池狀態(tài)參數(shù)緩慢變化至超過設(shè)定閾值的情況,辨識參數(shù)可能在閾值范圍上下波動,導(dǎo)致出現(xiàn)錯誤判斷[5].Zhao 等引入漸消因子運用自適應(yīng)擴(kuò)展卡爾曼濾波(Adaptive Fading Extended Kalman Filter,AFEKF)對鋰離子電池SOC 進(jìn)行估算,實驗證明漸消因子能夠有效提高卡爾曼濾波的估算精度[6].Gao等使用基于遺傳算法的神經(jīng)網(wǎng)絡(luò)對電池電壓傳感器、溫度傳感器、單體電池故障和電池內(nèi)部短路故障進(jìn)行了辨識,但該方法需要有大量優(yōu)質(zhì)的數(shù)據(jù)作為訓(xùn)練對象[7];Ablay 等分別使用EKF 與龍伯格觀測器(Luenberger Observers)生成電池內(nèi)阻、極化電容電壓與端電壓的殘差,并針對電壓、溫度傳感器的故障做出測試,但是故障狀態(tài)為軟件仿真,沒有實驗驗證[8-9].

    顯然,針對電池的基于端電壓殘差的故障診斷方法能夠很直觀地反映電池的運行狀態(tài).在系統(tǒng)處于正常狀態(tài)時,端電壓殘差為0;當(dāng)端電壓殘差出現(xiàn)較大的零點漂移時,表示電池狀態(tài)參數(shù)出現(xiàn)異常.得益于其較小的計算量與較高的估算精度,EKF 成為使用最多的殘差生成器,但是EKF 自適應(yīng)性和跟蹤性在面對各種不同的工況時無法滿足需求,因此,本文使用遞推最小二乘法對電池內(nèi)部參數(shù)進(jìn)行辨識,以獲得狀態(tài)空間中的未知參數(shù)R0、Rp、Cp;分別使用FEKF 和EKF 生成電池端路電壓的殘差,利用χ2檢驗獲得異常狀態(tài)判斷閾值,并以低溫異常狀態(tài)為例,對方法進(jìn)行實驗驗證.

    1 鋰離子電池建模

    電池模型分為3 類:電化學(xué)模型、經(jīng)驗?zāi)P团c等效電路模型.電化學(xué)模型根據(jù)電池內(nèi)部的電化學(xué)反應(yīng)建立,具有很高的精度,但計算復(fù)雜;經(jīng)驗?zāi)P椭饕揽看罅康臄?shù)據(jù),通過訓(xùn)練來達(dá)到較高的模型精度,對數(shù)據(jù)的數(shù)量和質(zhì)量都有很高要求;等效電路模型使用電容電阻和恒壓源等元器件構(gòu)成電路來模擬電池的運作機(jī)理,結(jié)構(gòu)簡單且靈活,便于在線計算與仿真[10-13].由于鋰離子電池內(nèi)部結(jié)構(gòu)可簡化為歐姆內(nèi)阻和極化內(nèi)阻,且一階模型的精度已經(jīng)能夠滿足需求,因此采用一階RC 等效電路模型來模擬鋰離子電池內(nèi)部結(jié)構(gòu),等效電路如圖1 所示.

    圖1 一階等效電路示意圖Fig.1 First-order equivalent circuit diagram

    根據(jù)基爾霍夫電壓定律,可得到關(guān)系式:

    電池的SOC 可表示為:

    經(jīng)過離散化后,式(3)可變?yōu)椋?/p>

    式中:t0為充放電的開始時刻;Cn為電池容量;T 為采樣周期;Rp為電池極化電阻;Cp為電池極化電容;I為電池工作電流;Vp為極化電容上的電壓;Vocv為電池開路電壓;R0為電池內(nèi)阻;Vt為電池端電壓.開路電壓Vocv與電池SOC 之間的關(guān)系可以通過對充放電實驗的數(shù)據(jù)進(jìn)行曲線擬合得到.

    2 FEKF 濾波

    EKF 在初始條件選取得當(dāng)?shù)那闆r下,能夠得到比較準(zhǔn)確的狀態(tài)估計值,但是任何一個實際系統(tǒng)都具有一定程度的不確定性,此時的EKF 濾波精度可能會下降甚至發(fā)散.為了克服這一缺點,引入漸消因子.

    2.1 狀態(tài)參數(shù)辨識

    使用RLS 對電池的歐姆內(nèi)阻R0、極化電阻Rp、極化電容Cp進(jìn)行在線辨識[14].該系統(tǒng)的離散數(shù)學(xué)模型為:

    式中:zk為系統(tǒng)輸出;μ 為前一時刻的輸入輸出向量,令μ=,Ik,Ik-1];θ0為系統(tǒng)真實狀態(tài)參數(shù)的代換向量,用以表示待辨識的狀態(tài)參數(shù):電池內(nèi)阻R0、極化電阻Rp和極化電容Cp.RLS 的遞推公式表達(dá)為:

    式中:K 為增益因子;P 為協(xié)方差矩陣;u 為遺忘因子;θ 為系統(tǒng)參數(shù)估計值向量,維度與μ 相同.u 越小跟蹤能力越強(qiáng),但波動也越大,一般取0.95 <u <1,為保證參數(shù)辨識的準(zhǔn)確性與穩(wěn)定性,取中間值,令u=0.985.

    首先進(jìn)行參數(shù)初始化,令μ、θ、K 的初始值均為0,P 的初始值則通過多次調(diào)試得到,為I1×3×106,I 為單位矩陣;然后將充放電過程中的V 與I 分別代入式(6),即可得到每一時刻的θ 值,完成參數(shù)辨識.

    2.2 狀態(tài)方程確立

    在早期,電池狀態(tài)參數(shù)的偏差通常較小,由于EKF 與UKF 的魯棒性較差,不易被系統(tǒng)模型識別,造成監(jiān)測診斷的延遲.因此,引入時變漸消因子λk.當(dāng)系統(tǒng)模型存在不確定性時,會首先在端電壓殘差序列rk中表現(xiàn)出來,導(dǎo)致≠0,即輸出端電壓殘差序列不正交.此時,通過漸消因子λk在線實時調(diào)整增益矩陣Kk,可以保證EKF 的輸出端電壓殘差序列相互正交,從而使EKF 具有強(qiáng)跟蹤濾波器的性質(zhì).由于狀態(tài)方程變量為一維,計算量很小,所以采用最優(yōu)漸消因子.

    設(shè)狀態(tài)方程變量xk=SOCk,系統(tǒng)狀態(tài)和觀測方程可表示為:

    式中:w 和v 是均值為0、方差分別為q 和r 的不相關(guān)高斯白噪聲.通過式(4),可以得到g 的表達(dá)式:

    通過式(1)(2)可以得到h 的表達(dá)式:

    引入漸消因子λk,令λ1=1,其它的漸消因子可由式(11)得到:

    式中:φ >0 為迭代步長,取值要求λk能夠快速衰減,而的計算方法為:

    至此可得到卡爾曼濾波的預(yù)測方程:

    迭代方程由式(14)表達(dá):

    式中:G 和H 分別為非線性函數(shù)g 與h 的一階Taylor展開,即得到端電壓殘差值為rk=-zk.

    3 狀態(tài)參數(shù)異常閾值計算

    當(dāng)電池狀態(tài)出現(xiàn)異常時,其狀態(tài)參數(shù)會出現(xiàn)變化,由此引起放電過程中端電壓殘差出現(xiàn)零點漂移,因此可用端電壓殘差數(shù)據(jù)的變化作為狀態(tài)參數(shù)異常判斷的依據(jù).

    影響狀態(tài)參數(shù)異常閾值η 的決定性因素為誤檢率PFA和漏檢率PM.若設(shè)定閾值較高,則漏檢率隨之升高,而誤檢率隨之降低;若設(shè)定閾值較低,則情況相反[15].構(gòu)造檢測函數(shù)為:

    式中:Tk為端電壓殘差序列的方差Rk.由于殘差序列為符合高斯分布的白噪聲,所以Dk符合χ2分布.取滑動窗口寬度為A,將窗口寬度內(nèi)采樣點Dk的平均值作為判斷該窗口內(nèi)數(shù)據(jù)是否出現(xiàn)狀態(tài)參數(shù)異常的依據(jù),即,依此得到誤檢率PFA的計算式[16-18]:

    式中:η 為狀態(tài)參數(shù)異常閾值.可通過初步設(shè)定PFA的值得到閾值大小.本系統(tǒng)中χ2分布的自由度為1,則對應(yīng)不同誤檢率的狀態(tài)參數(shù)異常閾值等同于χ2分布臨界值,如表1 所示.

    表1 誤檢率分布表Tab.1 False detection rate distribution

    當(dāng)滑動窗口內(nèi)檢測函數(shù)ηk大于η 時,判斷為系統(tǒng)狀態(tài)參數(shù)異常,否則系統(tǒng)狀態(tài)參數(shù)正常.判斷邏輯為:

    4 實驗設(shè)計與數(shù)據(jù)分析

    為了模擬緩變型狀態(tài)參數(shù)異常,將電池置入恒溫恒濕箱中緩慢降溫,在動態(tài)應(yīng)力測試工況(Dynamic Stress Test,DST)下進(jìn)行放電循環(huán).電池為億緯鋰能公司容量為2.55 Ah 的18650 型號三元鋰電池,充放電截止電壓分別為4.2 V 與2.75 V;恒溫箱為泓進(jìn)儀器HY-1010 型;充放電機(jī)為德普BTS20-5V/4×200A.實驗平臺原理圖如圖2 所示.

    圖2 充放電實驗平臺Fig.2 Charge and discharge experiment platform

    首先需要得到電池的OCV-SOC 曲線.在25 ℃時,電池充滿至4.2 V 截止電壓的情況下,對電池進(jìn)行恒流放電實驗,三元鋰電池的最大持續(xù)放電電流倍率為3 C,為了保證實驗安全取較小的放電電流倍率0.5 C,放電至電壓低于2.75 V,采樣周期為1 s,實驗測得OCV-SOC 曲線如圖3 所示.

    圖3 OCV-SOC 曲線Fig.3 OCV-SOC curve

    Matlab 的函數(shù)擬合方法主要包括插值擬合、多項式擬合與局部加權(quán)回歸擬合,其中多項式擬合可以通過調(diào)整階數(shù)防止出現(xiàn)過擬合與欠擬合的情況,并且多項式結(jié)構(gòu)簡潔、易于表達(dá).經(jīng)過調(diào)試,最終以七階多項式擬合得到OCV-SOC 對應(yīng)關(guān)系的多項式:

    式中:a7=-20.49,a6=74.31,a5=-106,a4=73.23,a3=-23.36,a2=2.453,a1=0.585 9,a0=0.344 5.

    然后將單體電池置入恒溫恒濕箱,將恒溫恒濕箱內(nèi)溫度穩(wěn)定至室溫25 ℃,開始進(jìn)行放電;1 h 后開始降溫,約在25 min 后降至5 ℃,5 min 后緩慢回溫至25 ℃,直至放電結(jié)束.期間溫度變化曲線如圖4所示,DST 放電循環(huán)的電流曲線如圖5 所示,充放電機(jī)的采樣周期為1 s,圖5 所示的1 000 個采樣點中包括兩個循環(huán)周期.

    實驗結(jié)束后,使用Matlab 做數(shù)據(jù)處理:

    1)將上位機(jī)采集到的實時電流電壓數(shù)據(jù)導(dǎo)入腳本文件,先通過遞推最小二乘法計算得到正常的電池狀態(tài)參數(shù)R0、Rp、Cp,時間常數(shù)τ=RpCp.

    2)將辨識出來的狀態(tài)參數(shù)與上位機(jī)采集到的電池端電壓與電流代入FEKF 計算,得到25 ℃降溫到5 ℃過程中的SOC 估算值以及端電壓殘差值r,并使用電池容量Cn減去每一時刻的電流與時間乘積,得到安時積分的SOC 真實值;所用到的Cn由實驗平臺的充放電機(jī)直接測算得到,為2.469 Ah,溫度變化的影響已被計算在內(nèi).

    圖4 溫度變化曲線Fig.4 Temperature change curve

    圖5 DST 放電電流曲線Fig.5 DST discharge current curve

    3)將端電壓殘差r 變換為χ2分布的Dk,取溫度降低到5 ℃后的1 500 s,即5 000~6 500 s 采樣點的Dk進(jìn)行對比.

    不同算法下,Dk的變化分別如圖6、圖7 所示;不同工況下不同算法的SOC 估算值與實際值分別如圖7~10 所示.

    從圖6、圖7 的對比可知,溫度降低至5 ℃后,端電壓殘差參數(shù)Dk在逐漸增大,在6 200 s 后達(dá)到峰值,但是通過EKF 生成的數(shù)據(jù)在6 200 s 前波動幅值不大,無法在較早的時間檢測出系統(tǒng)參數(shù)異常;而FEKF 生成的端電壓殘差數(shù)據(jù)Dk從5 700 s 開始便出現(xiàn)較大幅值的波動,能夠更及時地檢測出參數(shù)的異常.顯然,漸消因子的引入有效增大了異常狀態(tài)下的端電壓殘差,且未對正常狀態(tài)下的端電壓殘差造成影響.

    圖6 降溫流程結(jié)束時EKF 生成的不同時刻殘差參數(shù)DkFig.6 Residual parameter Dk at different times generated by EKF at the end of the cooling process

    圖7 降溫流程結(jié)束時FEKF 生成的不同時刻殘差參數(shù)DkFig.7 Residual parameter Dk at different times generated by FEKF at the end of the cooling process

    在溫度降至最低5 ℃的時刻約5 000 s 時,由于鋰離子電池在不同溫度下的擴(kuò)散作用的差異,電池的時間常數(shù)τ 與內(nèi)阻R0緩慢增大[19],根據(jù)式(1)(2)(14),在電池容量與電流大小不變的情況下,端電壓殘差會隨之增大,SOC 估算值隨之減小.因此,在溫度由25 ℃降低到5 ℃后,EKF 對SOC 的估算精度開始下降,逐漸偏離由安時積分計算得到的真實SOC曲線,如圖8 所示;而得益于漸消因子的引入,使用FEKF 進(jìn)行估算得到的SOC 與真實值之間的偏差很小,如圖9 所示,證明FEKF 在放大端電壓殘差的同時精度幾乎未受影響.

    分別取滑動窗口寬度A1=100,A2=50.取狀態(tài)參數(shù)異常閾值η=3.84,根據(jù)表1,此時誤檢率PFA=5%.每Ak個點為一組,計算每組Dk的平均值若大于η,則認(rèn)為Ak窗口寬度內(nèi),電池出現(xiàn)狀態(tài)參數(shù)異常,否則正常.取采樣時間為5 000~6 500 s,實時監(jiān)測結(jié)果如圖10~13 所示,概率為1 表示檢測結(jié)果為狀態(tài)參數(shù)異常,概率為0 則表示檢測結(jié)果為狀態(tài)參數(shù)正常.

    圖8 25 ℃降溫至5 ℃過程中使用EKF 的SOC 估算Fig.8 SOC estimation using EKF during cooling from 25 ℃to 5 ℃

    圖9 25 ℃降溫至5 ℃過程中使用FEKF 的SOC 估算Fig.9 SOC estimation using FEKF during cooling from 25 ℃to 5 ℃

    圖10 滑動窗口寬度為100 時EKF 生成殘差的狀態(tài)監(jiān)測結(jié)果Fig.10 Condition monitoring results of residuals generated by EKF when the data window width is 100

    圖11 滑動窗口寬度為100 時FEKF 生成殘差的狀態(tài)監(jiān)測結(jié)果Fig.11 Condition monitoring results of residuals generated by FEKF when the data window width is 100

    圖12 滑動窗口寬度為50 時EKF生成殘差的狀態(tài)監(jiān)測結(jié)果Fig.12 Condition monitoring results of residuals generated by EKF when the data window width is 50

    圖13 滑動窗口寬度為50 時FEKF 生成殘差的狀態(tài)監(jiān)測結(jié)果Fig.13 Condition monitoring results of residuals generated by FEKF when the data window width is 50

    從圖10、圖11 的對比可知,在A=100 時,使用EKF 生成的殘差只在6 200~6 300 s 內(nèi)出現(xiàn)大于閾值的偏移,而使用FEKF 生成的殘差在5 700~5 900 s、6 100~6 400 s 內(nèi)出現(xiàn)大于閾值的偏移;根據(jù)圖12、圖13,在A=50 時,使用EKF 生成的殘差在6 200~6 300 s 內(nèi)出現(xiàn)大于閾值的偏移,而FEKF生成的殘差分別在5 700~5 800 s、5 950~6 000 s、6 100~6 300 s、6 350~6 400 s 內(nèi)出現(xiàn)大于閾值的偏移.很顯然,使用EKF 得到的結(jié)果在6 200 s 才診斷出電池狀態(tài)參數(shù)異常,告警只維持了100 s;而使用FEKF 得到的結(jié)果,最早在5 700 s 便診斷出電池狀態(tài)參數(shù)異常,且在之后的600 s 內(nèi)持續(xù)告警.因此,使用FEKF 生成的殘差診斷出電池異常狀態(tài)的時間更早,告警持續(xù)時間更長,較之EKF 有更好的監(jiān)測性能.

    5 結(jié)論

    1)針對傳統(tǒng)EKF 無法及時識別鋰電池工作過程中出現(xiàn)的緩變型狀態(tài)參數(shù)異常的問題,提出了加入最優(yōu)漸消因子的FEKF,通過FEKF 生成殘差,構(gòu)造了χ2分布數(shù)列來計算誤檢率PFA以及對應(yīng)的狀態(tài)參數(shù)異常閾值η.

    2)將單體電池置于恒溫恒濕箱進(jìn)行緩慢降溫,按設(shè)置的工況,通過實驗測得了電池的電流電壓數(shù)據(jù),將數(shù)據(jù)導(dǎo)入Matlab 的算法程序中進(jìn)行計算,得到EKF 與FEKF 兩種不同算法下實時監(jiān)測的情況.證明了在相同的滑動窗口寬度以及閾值下,F(xiàn)EKF 相較傳統(tǒng)EKF 能夠更及時地識別出鋰電池運行時出現(xiàn)的緩變型狀態(tài)異常.

    猜你喜歡
    狀態(tài)參數(shù)端電壓充放電
    基于CKF的大型拖拉機(jī)狀態(tài)參數(shù)估計研究
    基于HHT算法的分布式電力系統(tǒng)狀態(tài)自動監(jiān)測方法
    堆石料的三維應(yīng)力分?jǐn)?shù)階本構(gòu)模型
    V2G模式下電動汽車充放電效率的研究
    基于SG3525的電池充放電管理的雙向DC-DC轉(zhuǎn)換器設(shè)計
    電子制作(2019年23期)2019-02-23 13:21:36
    The validity analysis of ground simulation test for non-ablative thermal protection materials
    勵磁調(diào)節(jié)器PT回路故障后機(jī)端電壓變化分析
    人工心臟無位置傳感器無刷直流電動機(jī)非導(dǎo)通相端電壓分析
    鋰離子電池充放電保護(hù)電路的研究
    鈦管蒸發(fā)器在酸性溶液中的防腐防垢
    免费看不卡的av| 人体艺术视频欧美日本| 精品一区二区免费观看| 久久久久国产精品人妻一区二区| 成人免费观看视频高清| 国产精品欧美亚洲77777| 天天躁夜夜躁狠狠躁躁| 国产成人精品久久二区二区91 | 电影成人av| 亚洲四区av| 国产男女超爽视频在线观看| 一二三四中文在线观看免费高清| 欧美人与善性xxx| 国产成人啪精品午夜网站| 女人高潮潮喷娇喘18禁视频| bbb黄色大片| 国产欧美日韩一区二区三区在线| 香蕉国产在线看| 秋霞在线观看毛片| 母亲3免费完整高清在线观看| 观看美女的网站| 国产成人一区二区在线| av又黄又爽大尺度在线免费看| 国产乱人偷精品视频| 国产麻豆69| 建设人人有责人人尽责人人享有的| 母亲3免费完整高清在线观看| av不卡在线播放| 免费观看人在逋| 欧美在线一区亚洲| 一本大道久久a久久精品| 欧美日韩综合久久久久久| 久久韩国三级中文字幕| 老汉色av国产亚洲站长工具| 一级毛片电影观看| 精品人妻一区二区三区麻豆| 久久久久精品久久久久真实原创| 欧美亚洲日本最大视频资源| 国产精品偷伦视频观看了| 亚洲欧美中文字幕日韩二区| 欧美日韩一区二区视频在线观看视频在线| 综合色丁香网| 看免费成人av毛片| 成人漫画全彩无遮挡| 久久精品aⅴ一区二区三区四区| 久久久国产一区二区| 18禁裸乳无遮挡动漫免费视频| www.熟女人妻精品国产| 久久久久精品久久久久真实原创| 你懂的网址亚洲精品在线观看| avwww免费| 国产精品久久久久久久久免| 午夜久久久在线观看| 国产探花极品一区二区| 午夜福利视频精品| 欧美xxⅹ黑人| 侵犯人妻中文字幕一二三四区| 91国产中文字幕| 中文字幕亚洲精品专区| 十分钟在线观看高清视频www| 丝袜喷水一区| 国产精品 欧美亚洲| 亚洲成国产人片在线观看| 中文字幕人妻熟女乱码| 女性生殖器流出的白浆| 精品亚洲成国产av| 九色亚洲精品在线播放| 国产99久久九九免费精品| 1024视频免费在线观看| 午夜精品国产一区二区电影| 各种免费的搞黄视频| 1024香蕉在线观看| 热re99久久国产66热| 精品免费久久久久久久清纯 | 另类精品久久| 性少妇av在线| 一级爰片在线观看| 男女高潮啪啪啪动态图| 亚洲精品美女久久久久99蜜臀 | 午夜老司机福利片| 久久婷婷青草| 日韩成人av中文字幕在线观看| 亚洲成色77777| 美女国产高潮福利片在线看| 亚洲精品视频女| 国产一级毛片在线| 久久这里只有精品19| 少妇猛男粗大的猛烈进出视频| 日韩一区二区三区影片| 精品人妻在线不人妻| 交换朋友夫妻互换小说| av网站在线播放免费| 欧美精品亚洲一区二区| 久久久精品免费免费高清| 成年美女黄网站色视频大全免费| 观看美女的网站| 啦啦啦啦在线视频资源| 人妻一区二区av| 欧美日韩成人在线一区二区| 十八禁网站网址无遮挡| 亚洲国产av影院在线观看| 欧美黑人精品巨大| 啦啦啦视频在线资源免费观看| 欧美xxⅹ黑人| 亚洲精品中文字幕在线视频| 日韩一卡2卡3卡4卡2021年| 美女午夜性视频免费| 天美传媒精品一区二区| 99国产综合亚洲精品| 国产精品国产三级国产专区5o| 国产精品无大码| 新久久久久国产一级毛片| 男女边吃奶边做爰视频| 在线精品无人区一区二区三| 精品久久久精品久久久| 最新的欧美精品一区二区| 伊人久久国产一区二区| 成人国语在线视频| 国产色婷婷99| av国产精品久久久久影院| 高清黄色对白视频在线免费看| 蜜桃国产av成人99| 大片电影免费在线观看免费| 久久鲁丝午夜福利片| www日本在线高清视频| 日本av免费视频播放| 亚洲成人免费av在线播放| 天美传媒精品一区二区| 中文字幕另类日韩欧美亚洲嫩草| 777久久人妻少妇嫩草av网站| 日韩伦理黄色片| 国产午夜精品一二区理论片| 中文字幕av电影在线播放| av网站免费在线观看视频| 亚洲专区中文字幕在线 | 久久久久国产精品人妻一区二区| 亚洲国产av新网站| 午夜福利,免费看| 老司机在亚洲福利影院| 日本wwww免费看| 91精品三级在线观看| 欧美日韩精品网址| 操美女的视频在线观看| 成年人免费黄色播放视频| 久久ye,这里只有精品| 国产野战对白在线观看| 国产成人欧美| 9191精品国产免费久久| 亚洲激情五月婷婷啪啪| 久久久久久久大尺度免费视频| 日韩人妻精品一区2区三区| 另类精品久久| 久久这里只有精品19| 国产爽快片一区二区三区| 大香蕉久久成人网| a 毛片基地| 亚洲av日韩在线播放| 久久毛片免费看一区二区三区| 国产爽快片一区二区三区| 国产精品嫩草影院av在线观看| 丝袜人妻中文字幕| 女性生殖器流出的白浆| 亚洲美女搞黄在线观看| av在线老鸭窝| av在线播放精品| 亚洲av中文av极速乱| 亚洲国产精品成人久久小说| 久久久久精品性色| 十八禁人妻一区二区| 黑人巨大精品欧美一区二区蜜桃| 曰老女人黄片| 建设人人有责人人尽责人人享有的| 黄色视频在线播放观看不卡| 亚洲五月色婷婷综合| 国产福利在线免费观看视频| 老司机靠b影院| 女人精品久久久久毛片| 精品一品国产午夜福利视频| 国产成人啪精品午夜网站| 欧美人与善性xxx| 99久国产av精品国产电影| 亚洲人成电影观看| 亚洲专区中文字幕在线 | 欧美少妇被猛烈插入视频| 亚洲男人天堂网一区| 国产精品久久久久成人av| 久久久久久久久久久久大奶| 婷婷色av中文字幕| 又粗又硬又长又爽又黄的视频| 亚洲成av片中文字幕在线观看| 好男人视频免费观看在线| 精品福利永久在线观看| 中文字幕制服av| 国产激情久久老熟女| 亚洲成av片中文字幕在线观看| 日本欧美国产在线视频| 黄片无遮挡物在线观看| 亚洲欧美中文字幕日韩二区| 日本av免费视频播放| 午夜福利乱码中文字幕| 丝袜美足系列| 免费人妻精品一区二区三区视频| 日本91视频免费播放| 丰满迷人的少妇在线观看| 欧美日韩亚洲综合一区二区三区_| 狂野欧美激情性bbbbbb| 丝袜人妻中文字幕| 国产xxxxx性猛交| 最近中文字幕高清免费大全6| 精品酒店卫生间| 日韩免费高清中文字幕av| 久久精品亚洲av国产电影网| 日韩 欧美 亚洲 中文字幕| 国产高清国产精品国产三级| av.在线天堂| 亚洲精品国产区一区二| 国产成人91sexporn| 在线天堂中文资源库| 天天躁日日躁夜夜躁夜夜| 欧美精品一区二区免费开放| 日韩欧美一区视频在线观看| 亚洲欧美一区二区三区久久| 亚洲av中文av极速乱| 亚洲精品国产av蜜桃| 久久ye,这里只有精品| 久久鲁丝午夜福利片| 激情视频va一区二区三区| 777米奇影视久久| 99精国产麻豆久久婷婷| 久久精品亚洲熟妇少妇任你| 观看av在线不卡| 亚洲欧美成人精品一区二区| 亚洲欧美中文字幕日韩二区| 亚洲伊人久久精品综合| 一区二区三区四区激情视频| 中文欧美无线码| 少妇人妻精品综合一区二区| 18禁裸乳无遮挡动漫免费视频| 国产一级毛片在线| 如日韩欧美国产精品一区二区三区| 高清在线视频一区二区三区| 在线观看免费高清a一片| 大香蕉久久成人网| 国产精品国产av在线观看| 一边摸一边抽搐一进一出视频| 伦理电影免费视频| 国产免费一区二区三区四区乱码| 老司机影院成人| 日日爽夜夜爽网站| av天堂久久9| 午夜激情久久久久久久| 精品午夜福利在线看| 老鸭窝网址在线观看| 欧美97在线视频| 国产精品嫩草影院av在线观看| 亚洲国产欧美在线一区| 肉色欧美久久久久久久蜜桃| 一级爰片在线观看| 色婷婷久久久亚洲欧美| 中文精品一卡2卡3卡4更新| 波野结衣二区三区在线| 伊人久久国产一区二区| 国语对白做爰xxxⅹ性视频网站| av国产精品久久久久影院| 一级毛片电影观看| 又粗又硬又长又爽又黄的视频| 免费不卡黄色视频| 亚洲色图 男人天堂 中文字幕| 日韩av不卡免费在线播放| 中文欧美无线码| 老鸭窝网址在线观看| 中文字幕另类日韩欧美亚洲嫩草| 精品少妇黑人巨大在线播放| 亚洲国产毛片av蜜桃av| 午夜精品国产一区二区电影| 国产av国产精品国产| 国产精品人妻久久久影院| av.在线天堂| 日韩精品免费视频一区二区三区| 久久久久精品性色| 国产人伦9x9x在线观看| 香蕉国产在线看| 亚洲精品美女久久av网站| 婷婷色av中文字幕| 在线观看免费高清a一片| 亚洲一级一片aⅴ在线观看| 亚洲欧美一区二区三区黑人| 国产成人系列免费观看| 亚洲精品日本国产第一区| 国产成人啪精品午夜网站| 欧美精品一区二区免费开放| 国产在线免费精品| 国产精品一区二区精品视频观看| 中国国产av一级| 免费观看a级毛片全部| 欧美最新免费一区二区三区| 极品少妇高潮喷水抽搐| 亚洲国产精品成人久久小说| 18禁裸乳无遮挡动漫免费视频| 国产成人精品无人区| 如何舔出高潮| 国产一区二区在线观看av| 亚洲国产毛片av蜜桃av| 国产成人啪精品午夜网站| 欧美黄色片欧美黄色片| 国产99久久九九免费精品| 国产野战对白在线观看| 午夜激情久久久久久久| 人人妻,人人澡人人爽秒播 | 国产精品av久久久久免费| 在线观看一区二区三区激情| 成年人午夜在线观看视频| 一区二区三区乱码不卡18| 国产片内射在线| 熟女少妇亚洲综合色aaa.| 最新的欧美精品一区二区| 亚洲av电影在线进入| 九草在线视频观看| 另类精品久久| 我要看黄色一级片免费的| 亚洲av电影在线观看一区二区三区| 观看av在线不卡| 欧美日韩一区二区视频在线观看视频在线| 三上悠亚av全集在线观看| 热99国产精品久久久久久7| 天堂8中文在线网| 日韩中文字幕视频在线看片| 日韩av在线免费看完整版不卡| 电影成人av| 日韩一本色道免费dvd| 国产女主播在线喷水免费视频网站| 亚洲欧美精品自产自拍| 国产免费福利视频在线观看| 九草在线视频观看| 日韩熟女老妇一区二区性免费视频| 国产精品香港三级国产av潘金莲 | 久久精品亚洲av国产电影网| 精品国产乱码久久久久久男人| 久久天堂一区二区三区四区| 亚洲av在线观看美女高潮| 最近的中文字幕免费完整| 亚洲欧美一区二区三区久久| 青春草亚洲视频在线观看| 亚洲一级一片aⅴ在线观看| 香蕉国产在线看| 母亲3免费完整高清在线观看| 老汉色av国产亚洲站长工具| 久久婷婷青草| 精品人妻一区二区三区麻豆| 国产精品99久久99久久久不卡 | 热99国产精品久久久久久7| 久久99热这里只频精品6学生| 熟女少妇亚洲综合色aaa.| 久久久精品国产亚洲av高清涩受| 国产精品国产三级专区第一集| 欧美激情 高清一区二区三区| 少妇精品久久久久久久| 韩国av在线不卡| 女性被躁到高潮视频| 中文精品一卡2卡3卡4更新| 欧美日韩一级在线毛片| a级毛片在线看网站| 欧美激情高清一区二区三区 | svipshipincom国产片| 久久久国产一区二区| 色视频在线一区二区三区| 一级毛片我不卡| 侵犯人妻中文字幕一二三四区| 三上悠亚av全集在线观看| 黄色 视频免费看| 欧美日韩精品网址| a级毛片黄视频| 我要看黄色一级片免费的| netflix在线观看网站| 91老司机精品| 成人18禁高潮啪啪吃奶动态图| 美女扒开内裤让男人捅视频| 少妇的丰满在线观看| 七月丁香在线播放| 亚洲欧美成人综合另类久久久| 欧美精品一区二区免费开放| 母亲3免费完整高清在线观看| 十八禁人妻一区二区| 国产男女内射视频| 婷婷成人精品国产| 美女高潮到喷水免费观看| 午夜久久久在线观看| 国产精品一二三区在线看| 两个人免费观看高清视频| 国产成人精品久久二区二区91 | 久久午夜综合久久蜜桃| 亚洲熟女毛片儿| 九草在线视频观看| 国产麻豆69| 蜜桃国产av成人99| 人人妻人人澡人人爽人人夜夜| 菩萨蛮人人尽说江南好唐韦庄| 最新在线观看一区二区三区 | 国产高清不卡午夜福利| 一级毛片 在线播放| 久久人人爽av亚洲精品天堂| 一边摸一边抽搐一进一出视频| 免费黄网站久久成人精品| 亚洲视频免费观看视频| 十八禁人妻一区二区| 久久精品久久精品一区二区三区| 极品少妇高潮喷水抽搐| 香蕉丝袜av| 成年女人毛片免费观看观看9 | 乱人伦中国视频| 亚洲精品美女久久久久99蜜臀 | 成年美女黄网站色视频大全免费| av一本久久久久| 国产精品国产三级国产专区5o| 欧美精品一区二区大全| 91老司机精品| 亚洲国产精品国产精品| 免费女性裸体啪啪无遮挡网站| 少妇人妻 视频| 美女视频免费永久观看网站| 日韩熟女老妇一区二区性免费视频| 少妇人妻 视频| 国产精品嫩草影院av在线观看| 久久久久久久久久久免费av| 色综合欧美亚洲国产小说| 老司机靠b影院| 90打野战视频偷拍视频| 色综合欧美亚洲国产小说| 丰满饥渴人妻一区二区三| 天天添夜夜摸| 菩萨蛮人人尽说江南好唐韦庄| 超碰97精品在线观看| 国产伦人伦偷精品视频| 男女高潮啪啪啪动态图| 婷婷色av中文字幕| 最近的中文字幕免费完整| 亚洲国产成人一精品久久久| 欧美激情高清一区二区三区 | 精品国产超薄肉色丝袜足j| 99精国产麻豆久久婷婷| 精品午夜福利在线看| 亚洲国产欧美日韩在线播放| 在线观看免费日韩欧美大片| 国产老妇伦熟女老妇高清| 欧美日韩亚洲综合一区二区三区_| 91精品伊人久久大香线蕉| 99re6热这里在线精品视频| 在线观看国产h片| 日韩电影二区| 中文字幕色久视频| 黑人巨大精品欧美一区二区蜜桃| 欧美日韩一级在线毛片| 老汉色∧v一级毛片| av网站免费在线观看视频| 欧美国产精品va在线观看不卡| 天天添夜夜摸| 午夜日韩欧美国产| 国产成人精品福利久久| 国产精品女同一区二区软件| 搡老岳熟女国产| 一级毛片 在线播放| 免费黄网站久久成人精品| 青草久久国产| 欧美日韩av久久| 你懂的网址亚洲精品在线观看| 中文字幕av电影在线播放| 熟女av电影| 欧美另类一区| 男人添女人高潮全过程视频| 菩萨蛮人人尽说江南好唐韦庄| 蜜桃国产av成人99| 日韩大码丰满熟妇| 欧美日韩av久久| 1024香蕉在线观看| av.在线天堂| 丰满乱子伦码专区| 久热爱精品视频在线9| 国产伦人伦偷精品视频| 久久ye,这里只有精品| 悠悠久久av| 国产黄频视频在线观看| 一区二区三区乱码不卡18| 赤兔流量卡办理| 精品午夜福利在线看| 韩国高清视频一区二区三区| 90打野战视频偷拍视频| 在线观看三级黄色| 无遮挡黄片免费观看| 别揉我奶头~嗯~啊~动态视频 | 国产97色在线日韩免费| 国产亚洲av片在线观看秒播厂| 中文精品一卡2卡3卡4更新| tube8黄色片| 成人国产麻豆网| 九九爱精品视频在线观看| 在线观看免费午夜福利视频| 日本欧美视频一区| 国产精品人妻久久久影院| 宅男免费午夜| 欧美精品av麻豆av| 999久久久国产精品视频| 操出白浆在线播放| 国产片特级美女逼逼视频| 免费日韩欧美在线观看| 久久久久网色| 蜜桃在线观看..| 精品国产一区二区久久| 午夜福利视频精品| 欧美人与善性xxx| 可以免费在线观看a视频的电影网站 | 久久免费观看电影| 啦啦啦视频在线资源免费观看| 日韩视频在线欧美| 一区在线观看完整版| 又黄又粗又硬又大视频| 人体艺术视频欧美日本| 丰满乱子伦码专区| 99久久人妻综合| 午夜福利网站1000一区二区三区| 在现免费观看毛片| 亚洲欧美成人精品一区二区| 一级毛片 在线播放| 夜夜骑夜夜射夜夜干| 最近中文字幕2019免费版| 中文字幕色久视频| 免费看av在线观看网站| 国产av精品麻豆| 国语对白做爰xxxⅹ性视频网站| 午夜免费观看性视频| 波多野结衣一区麻豆| 亚洲熟女毛片儿| 午夜福利,免费看| 亚洲成人国产一区在线观看 | 亚洲av福利一区| 欧美成人午夜精品| 国产探花极品一区二区| 亚洲天堂av无毛| 2021少妇久久久久久久久久久| 我的亚洲天堂| 波多野结衣av一区二区av| 国产极品天堂在线| 亚洲精品第二区| 日韩熟女老妇一区二区性免费视频| 国产精品 欧美亚洲| avwww免费| 国产精品秋霞免费鲁丝片| 亚洲国产精品成人久久小说| 亚洲国产日韩一区二区| 黄色怎么调成土黄色| 亚洲第一av免费看| 国产97色在线日韩免费| 久久久久国产精品人妻一区二区| 亚洲人成77777在线视频| 啦啦啦视频在线资源免费观看| 人体艺术视频欧美日本| 丝袜美足系列| 欧美日韩一区二区视频在线观看视频在线| 黄网站色视频无遮挡免费观看| 91aial.com中文字幕在线观看| 丝袜美足系列| 国产极品天堂在线| 高清欧美精品videossex| 日本av免费视频播放| 悠悠久久av| 亚洲精品自拍成人| 涩涩av久久男人的天堂| 久久久久精品久久久久真实原创| 天堂俺去俺来也www色官网| 国产精品无大码| 妹子高潮喷水视频| 国产成人欧美在线观看 | 久久久久人妻精品一区果冻| 精品一品国产午夜福利视频| 国产激情久久老熟女| 99re6热这里在线精品视频| 亚洲国产精品一区三区| 中文字幕另类日韩欧美亚洲嫩草| 欧美 亚洲 国产 日韩一| 波野结衣二区三区在线| 建设人人有责人人尽责人人享有的| 日本色播在线视频| 在线观看免费视频网站a站| 中文精品一卡2卡3卡4更新| 男女无遮挡免费网站观看| 欧美另类一区| 伊人久久国产一区二区| 自拍欧美九色日韩亚洲蝌蚪91| 国产亚洲av片在线观看秒播厂| e午夜精品久久久久久久| 亚洲三区欧美一区| 两性夫妻黄色片| 18禁动态无遮挡网站| 一级毛片电影观看| 人体艺术视频欧美日本| 热99国产精品久久久久久7| 日韩伦理黄色片| 亚洲视频免费观看视频| 久久青草综合色| 久久免费观看电影| 青春草国产在线视频| 18禁裸乳无遮挡动漫免费视频| 不卡视频在线观看欧美| 丰满少妇做爰视频| 美女中出高潮动态图| 久久久精品国产亚洲av高清涩受| 欧美激情 高清一区二区三区| 老司机亚洲免费影院| 男的添女的下面高潮视频| 99久久人妻综合| 又粗又硬又长又爽又黄的视频| 国产免费现黄频在线看| 亚洲自偷自拍图片 自拍| 久久精品久久久久久噜噜老黄| 一个人免费看片子| 国产1区2区3区精品|