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

    自動引導(dǎo)車電池SOC估算方法研究

    2017-11-01 07:18:11吳鐵洲曾藝師熊金龍
    計算機(jī)測量與控制 2017年8期
    關(guān)鍵詞:積分法初值卡爾曼濾波

    吳鐵洲,張 敏,曾藝師,熊金龍

    (1.湖北工業(yè)大學(xué) 太陽能高效利用湖北省協(xié)同創(chuàng)新中心, 武漢 430068;2.國網(wǎng)湖北省電力公司隨州供電公司,湖北 隨州 441300;3.國網(wǎng)湖北省電力公司檢修公司,武漢 430051)

    自動引導(dǎo)車電池SOC估算方法研究

    吳鐵洲1,張 敏1,曾藝師2,熊金龍3

    (1.湖北工業(yè)大學(xué) 太陽能高效利用湖北省協(xié)同創(chuàng)新中心, 武漢 430068;2.國網(wǎng)湖北省電力公司隨州供電公司,湖北 隨州 441300;3.國網(wǎng)湖北省電力公司檢修公司,武漢 430051)

    自動引導(dǎo)車(AGV車)工況特殊,電流積分法估算電池剩余容量(SOC)誤差較大,而且存在累積誤差;為了提高AGV車電池剩余容量估算的準(zhǔn)確度,對擴(kuò)展卡爾曼濾波法估算AGV車電池剩余容量進(jìn)行了研究,分析了AGV車特殊工況,提出將擴(kuò)展卡爾曼濾波法的濾波增益改進(jìn)為動態(tài)調(diào)整濾波增益,有效提高擴(kuò)展卡爾曼濾波法的跟蹤效果;實驗表明使用擴(kuò)展卡爾曼濾波法估算AGV車電池剩余容量精度較高,采用動態(tài)校正的濾波增益提高了估算過程的跟蹤效果,解決了AGV車電池剩余容量估算不準(zhǔn)確的問題。

    電池剩余容量;自動引導(dǎo)車;卡爾曼濾波

    0 引言

    自動引導(dǎo)車(AGV)是一種自動化的搬運設(shè)備,主要應(yīng)用在工業(yè)生產(chǎn)中[1]。AGV車一般使用蓄電池作為電源,在AGV車能量管理系統(tǒng)中準(zhǔn)確估算SOC是至關(guān)重要的,可以保證AGV車高效安全地運行,同時也能夠避免由于過充或過放影響AGV車電池的使用壽命。但是蓄電池的剩余容量和很多因素有關(guān),而且AGV車工況特殊,這就使得蓄電池SOC值估算變得復(fù)雜,準(zhǔn)確度難以得到保證[2]。

    SOC估算方法主要有開路電壓法,電流積分法,卡爾曼濾波法,神經(jīng)網(wǎng)絡(luò)法。AGV車中常用的是電流積分法,但是AGV車工況特殊,電流積分法估算SOC值的準(zhǔn)確度并不理想,而且會產(chǎn)生累積誤差[3]。雖然在實際工業(yè)應(yīng)用中對電流積分法估算的SOC值會進(jìn)行校對,但是并沒有從根本上解決電流積分法依賴初值,累積誤差較大的問題。開路電壓法無法在線實時估算SOC值,神經(jīng)網(wǎng)絡(luò)法對訓(xùn)練方法和訓(xùn)練數(shù)據(jù)依賴性很大,目前還沒有得到很好的使用??柭鼮V波法估算電池SOC不依賴于初值,不會產(chǎn)生誤差[4],適合應(yīng)用在電流變化幅度較大的工況中,綜合分析卡爾曼濾波法更加適合應(yīng)用在AGV車SOC估算中。

    本文首先分析了AGV車在工業(yè)生產(chǎn)中特殊的工況,在此基礎(chǔ)上分析了電流積分法估算SOC值不準(zhǔn)的根本原因,提出使用擴(kuò)展卡爾曼濾波法估算AGV車SOC值,針對擴(kuò)展卡爾曼濾波法在工業(yè)生產(chǎn)中跟蹤效果差的問題,本文對其進(jìn)行優(yōu)化,改進(jìn)其濾波增益為動態(tài)變化增益,提高擴(kuò)展卡爾曼濾波法應(yīng)用在AGV車SOC估算上的跟蹤效果。最后通過一系列實驗驗證了擴(kuò)展卡爾曼濾波算法的準(zhǔn)確性和有效性,驗證了改進(jìn)擴(kuò)展卡爾曼濾波算法良好的跟蹤效果。

    1 AGV車工況分析

    AGV車運行在工業(yè)生產(chǎn)線中,其工作節(jié)奏較快,一般不會預(yù)留專門的充電時間,只有在等待工序完成的時間里才進(jìn)行充電,而且充電電流大,充電倍率可以達(dá)到1~2 C,充電時間短,一般是幾秒或者是幾十秒。而AGV車在正常工作時電流較小,一般是在0.5C之內(nèi)。所以AGV車的工況就可以概括為充電電流大,放電電流小,充電時間短,充電頻率高。圖1描述了AGV車在12分鐘內(nèi)充放電變化情況,放電電流為正,充電電流為負(fù)。從圖中可以看出AGV車在10 min內(nèi)就充電兩次,而且每次充電電流達(dá)到了100 A,而放電電流僅有5 A左右。在這種電流波動幅度大,充電頻率高的情況下如何準(zhǔn)確估算SOC值是一個值得研究的問題。

    圖1 AGV車實際工況下的電流變化圖

    電流積分法是估算SOC值的常用方法,但是電流積分法依賴初值,對電流的測量精度要求較高,AGV車在特殊工況下SOC初始值和電流測量精度都難以得到保證,所以電流積分法估算AGV車SOC值的準(zhǔn)確度不是很高??柭鼮V波法的基本思想是一個不斷加權(quán)迭代的過程,在整個濾波的過程中有兩個值,分別是模型的預(yù)測值和儀器的測量值,然后預(yù)測值和測量值進(jìn)行綜合加權(quán)得到了最優(yōu)估計值。并且卡爾曼濾波法不依賴初始值,不會產(chǎn)生累積誤差,能夠保證特殊工況下SOC值的準(zhǔn)確性。但是卡爾曼濾波法依賴電池模型,電池模型的精度直接影響SOC值估算的精度,所以準(zhǔn)確建立電池模型是卡爾曼濾波法估算SOC值的基礎(chǔ)。

    2 蓄電池模型選取

    本文采用Thevenin模型對AGV車電池進(jìn)行等效。該模型可以將鋰離子電池充放電過程中的突變特性和漸變特性都描述出來,其等效電路原理圖如圖2所示。其中R0描述了鋰電池充放電過程中突變的電阻特性,Rp和Cp模擬了電池漸變的電容特性。電池等效電路模型中的參數(shù)可以通過對電池進(jìn)行間歇脈沖充放電實驗進(jìn)行辨識。通過實驗中電池在充放電過程中的電壓升降及充放電時間計算電路模型中的內(nèi)阻和電容等參數(shù)[5]。

    圖2 電池等效電路原理圖

    圖2中,E是電池電動勢,R0是電池內(nèi)阻,Rp和Cp是極化電阻和極化電容,Uoc是電池開路電壓,I是電池總電流,Ip是流過極化電阻的電流,Up是極化電阻和極化電容兩端的電壓。由卡爾曼濾波法的原理可知在濾波的過程中需要一個狀態(tài)量和一個觀測量,而SOC值并不能直接測量得到,所以選取Up為狀態(tài)變量,經(jīng)過線性化和離散化處理狀態(tài)方程如式(1)所示。選取檢測到的電池端電壓Uoc為觀測量,建立的觀測方程如式(2)所示:

    (1)

    (2)

    式(1)中,Up,k表示的是k時刻Up的值,Up,k-1表示的是k-1時刻Up的值,同樣SOCk表示的是k時刻電池SOC值,SOCk-1表示的是k-1時刻SOC值,T是系統(tǒng)的采樣時間,w1,k-1和w2,k-1為系統(tǒng)噪聲。由式(1)和(2)可得卡爾曼濾波所需的矩陣如式(3):

    (3)

    3 改進(jìn)擴(kuò)展卡爾曼濾波算法估算AGV車SOC值

    卡爾曼濾波法是針對線性系統(tǒng)進(jìn)行最優(yōu)估計,擴(kuò)展卡爾曼濾波算法是對傳統(tǒng)卡爾曼濾波算法的推廣,適用于非線性系統(tǒng)的最優(yōu)估計。蓄電池是典型的非線性系統(tǒng),所以需要使用擴(kuò)展卡爾曼濾波法來估算其剩余容量[6]。通過蓄電池的建模得到了系統(tǒng)的狀態(tài)方程和觀測方程,結(jié)合卡爾曼濾波法的特征方程就可以得到如式(4)所示的卡爾曼濾波過程。

    Pk/k=(I-KkCk)Pk/k-1

    (4)

    濾波方程的初始條件有:SOC的初始值SOC0,Cp兩端的初始電壓Up(0)以及預(yù)估誤差的協(xié)方差P0。SOC0通過讀取上一次的SOC值獲得,Up(0)在初始化時一般為0,預(yù)估誤差的協(xié)方差初始值P0根據(jù)工程實踐經(jīng)驗選取合適的值。通常情況下,這些初值的準(zhǔn)確度對估算的SOC值影響并不是很大,因為卡爾曼濾波器對精確度較差的初值有很好的魯棒性,能夠很快收斂到真實值,如果初值取值合理那么卡爾曼濾波收斂到真實值的速度會比較快,如果取值不合理收斂速度會相對比較慢,需要多次迭代才能收斂到真實值,但是最終都會消除初始值不準(zhǔn)確帶來的誤差[7],這就是擴(kuò)展卡爾曼濾波法估算AGV車SOC值不依賴初值的根本原因。擴(kuò)展卡爾曼濾波法估算電池SOC就是當(dāng)k=1,2,3……對式(4)進(jìn)行循環(huán)迭代,使?fàn)顟B(tài)量收斂到真實值。

    具體估算過程如下:

    1) 根據(jù)式(4)中的第一個公式,通過電池等效模型進(jìn)行預(yù)測,通過k-1時刻的狀態(tài)量預(yù)測出k時刻的狀態(tài)量;

    2) 根據(jù)式(4)中的第二個公式,計算預(yù)測值的誤差協(xié)方差矩陣,確定預(yù)測的準(zhǔn)確度,為后邊進(jìn)行權(quán)重分配做準(zhǔn)備;

    3) 根據(jù)式(4)中的第三個公式計算卡爾曼濾波增益,也就是給實測值和預(yù)測值分配權(quán)重,進(jìn)而通過卡爾曼濾波增益修正狀態(tài)量;

    4) 根據(jù)式(4)中的第四個公式,通過該式進(jìn)行卡爾曼的最優(yōu)估計,利用預(yù)測狀態(tài)信息反算觀測量,然后計算實測值Zk與反算值之差,再進(jìn)行綜合加權(quán);

    5) 根據(jù)式(4)中的第五個公式,該式就是對計算出來的最優(yōu)結(jié)果進(jìn)行評價,計算一下最優(yōu)結(jié)果的準(zhǔn)確度如何,也是為了下一時刻做準(zhǔn)備,保證后續(xù)可以一直循環(huán)下去。

    卡爾曼濾波法估算電池SOC值就是通過上邊5個步驟不停地進(jìn)行迭代,從而使SOC值收斂于真實值,而且估算的過程中對初值沒有依賴性,也不存在累積誤差,估算的SOC值準(zhǔn)確度和穩(wěn)定性高于電流積分法,但是AGV車工況特殊,電流突變頻率高,在實際工況下擴(kuò)展卡爾曼濾波法估算SOC在電流突變時跟蹤效果較差,影響SOC值估算的準(zhǔn)確性和實時性。所以需要對擴(kuò)展卡爾曼濾波法進(jìn)行改進(jìn),加快算法在電流突變時的預(yù)測速度,提高其在特殊工況下估算SOC值的準(zhǔn)確性和實時性。

    加快算法在電流突變時的預(yù)測速度就要求算法在電流突變剛開始的時刻加大SOC值的修正幅度。由式(4)得出增大濾波增益Kk的值可以提高狀態(tài)量的更新幅度,加快算法的跟蹤速度。所以要改善特殊工況下擴(kuò)展卡爾曼濾波算法的預(yù)測速度就需要在突變開始時增大濾波增益,突變過后減小濾波增益,從而達(dá)到動態(tài)修正SOC的作用。

    令改進(jìn)后的濾波增益為Kk’則擴(kuò)展卡爾曼濾波算法的狀態(tài)更新如式(5)所示:

    (5)

    由式(5)可得:

    (6)

    由式(4)可得:

    (7)

    令濾波增益修正函數(shù)為W(t),則由式(6)和式(7)可得:

    (8)

    濾波增益在突變開始時增大,提高修正幅度,隨著突變的進(jìn)行濾波增益逐漸減小,在突變結(jié)束后恢復(fù)正常值,為了找出特殊工況下濾波增益修正函數(shù)的這種變化規(guī)律,本文選取實驗數(shù)據(jù)然后通過數(shù)據(jù)擬合找出修正函數(shù)。實驗數(shù)據(jù)如表1所示。

    使用Matlab最小二乘法擬合實驗數(shù)據(jù)可得如式(9)所示的函數(shù)關(guān)系。

    W(t)=1.011e-0.4909t+0.9998

    (9)

    將式(9)化簡可得:

    W(t)=0.6121t+0.9998 (10)

    在Matlab中繪制實驗數(shù)據(jù)曲線圖,如圖3所示。在圖3中前6秒是電流的突變時間,在突變剛開始時修正函數(shù)使得濾波增益值較大,對SOC值的修正幅度也比較大,突變過后濾波增益值恢復(fù)正常值。這樣濾波增益就可以根據(jù)突變的情況實現(xiàn)動態(tài)調(diào)整,加快算法的跟蹤速度。

    圖3 濾波增益修正函數(shù)變化曲線

    通過實驗數(shù)據(jù)的擬合可以得出修正函數(shù)符合指數(shù)函數(shù)的變化規(guī)律,令其動態(tài)修正函數(shù)如式(11)所示:

    δ=1+φβt

    (11)

    則改進(jìn)后擴(kuò)展卡爾曼濾波算法中的增益Kk’如式(12)所示

    Kk′ = (1 + ?βt)Pk/k-1CkT(CkPk/k-1CkT+Rk)-1

    (12)

    式(12)中β是一個調(diào)整因子,其值是0-1的數(shù),該參數(shù)可以調(diào)整電流突變持續(xù)的時間,如果突變時間持續(xù)較短則該參數(shù)取值較小,如果突變時間持續(xù)時間較長該參數(shù)取值較大,具體取值可根據(jù)實際工程確定。φ是調(diào)節(jié)幅度因子,該參數(shù)可以跟蹤系統(tǒng)的突變強度,如果突變強度較大該參數(shù)取值較大,突變強度較小取值較小。式中t是系統(tǒng)突變的持續(xù)時間,設(shè)t0是突變開始的時間,而t’是突變結(jié)束的時間,則t=t0-t’。由此可計算出改進(jìn)擴(kuò)展卡爾曼濾波算法中增益變化如式(13)所示。

    δ=1+φ(t→0)

    δ=1(t→∞)

    Kk'=(1+φ)Kk(t→0)

    Kk'=Kk(t→∞)

    (13)

    改進(jìn)擴(kuò)展卡爾曼濾波增益后,系統(tǒng)在突變剛開始的時候?qū)V波增益擴(kuò)大到?jīng)]有改進(jìn)前的1+φ倍,系統(tǒng)SOC值的修正幅度變大,迅速收斂到真實值;隨著突變的結(jié)束,濾波增益逐漸變小,最終恢復(fù)為Kk。這樣就實現(xiàn)了擴(kuò)展卡爾曼濾波算法中濾波增益動態(tài)調(diào)整,修正幅度動態(tài)變化,有利于系統(tǒng)在突變情況下更快地收斂到真是值,減小擴(kuò)展卡爾曼濾波算法中觀測值滯后影響SOC估算精度的影響,提高了系統(tǒng)的跟蹤能力。

    4 實驗與分析

    為了驗證擴(kuò)展卡爾曼濾波法估算AGV車SOC值的效果,同時還要驗證改進(jìn)后的擴(kuò)展卡爾曼濾波算法的跟蹤效果,本研究進(jìn)行了一系列的實驗,實驗對比分析了擴(kuò)展卡爾曼濾波法與安時法估算AGV車SOC值的效果差異,對比分析了擴(kuò)展卡爾曼濾波算法改進(jìn)前后的跟蹤效果。

    AGV車在正常運行的情況下,電池的各種數(shù)據(jù)是實時記錄的,所以本實驗使用的數(shù)據(jù)都是AGV車實時運行的現(xiàn)場數(shù)據(jù)。本文在Visual Studio 2012平臺上編寫了算法驗證程序,通過讀取AGV車實時運行的現(xiàn)場數(shù)據(jù)來模擬出AGV車的實際工況,然后分別使用電流積分法和擴(kuò)展卡爾曼濾波法估算AGV車的SOC值,其對比分析曲線如圖4所示。

    圖4 SOC值估算對比曲線

    在實驗中為了保證安時法的準(zhǔn)確性,給定其初值是準(zhǔn)確的,而卡爾曼濾波法的初值有較大誤差。通過實驗結(jié)果可知安時法誤差較大,最大誤差超過7%,而且這種誤差范圍是在初值準(zhǔn)確的情況下產(chǎn)生的,如果初值不準(zhǔn)確其誤差范圍還會進(jìn)一步擴(kuò)大。而卡爾曼濾波法在初值不準(zhǔn)確的情況下,隨著迭代次數(shù)的增加其估算結(jié)果逐漸收斂到真實值,誤差范圍基本可以控制在4%以內(nèi),準(zhǔn)確度明顯高于安時積分法。但是從結(jié)果對比曲線上也可以看出擴(kuò)展卡爾曼濾波法也存在一定的問題。在圖5a中標(biāo)出了兩個個圈,在每一個標(biāo)識圈內(nèi)可以看出卡爾曼濾波法估算的SOC值都要滯后于真實值,導(dǎo)致誤差比較大,這是因為當(dāng)AGV車電流發(fā)生突變時,擴(kuò)展卡爾曼濾波法的觀測電壓是滯后的,從而導(dǎo)致SOC的估算也出現(xiàn)滯后。為了解決這個問題,對擴(kuò)展卡爾曼濾波法進(jìn)行了改進(jìn),即將擴(kuò)展卡爾曼濾波算法中的調(diào)整增益設(shè)置成動態(tài)增益,在該實驗中β=0.5,φ=1,對改進(jìn)后的算法進(jìn)行實驗,結(jié)果如圖5(b)所示。

    圖5 擴(kuò)展卡爾曼濾波算法跟蹤效果圖

    從改進(jìn)后的實驗結(jié)果可知SOC估算的精度進(jìn)一步提高,相比改進(jìn)前的卡爾曼濾波算法其跟蹤效果更好,在電流發(fā)生突變時也能較快地收斂到真實值。為進(jìn)一步說明優(yōu)化后的卡爾曼濾波算法其跟蹤效果更好,將SOC真實值,改進(jìn)前卡爾曼濾波法估算的SOC值,改進(jìn)后卡爾曼濾波法估算的SOC值進(jìn)行放大比較,結(jié)果如圖6所示。從圖6中可以很清楚地看到改進(jìn)后的卡爾曼濾波法其跟蹤效果相比沒有改進(jìn)的卡爾曼濾波法提升了30%左右,所以在AGV車這種電流變化幅度大的工況下,改進(jìn)后的卡爾曼濾波法估算的SOC值更加準(zhǔn)確。

    圖6 跟蹤效果對比圖

    5 結(jié)論

    本文立足解決工業(yè)生產(chǎn)實際問題,充分分析了AGV車特殊的工況,指出安時積分法在估算其SOC值不準(zhǔn)確的原因所在,提出使用擴(kuò)展卡爾曼濾波法估算其電池SOC值,針對傳統(tǒng)擴(kuò)展卡爾曼濾波法跟蹤效果差的問題,改進(jìn)其濾波增益為動態(tài)變化增益,提高估算SOC值的實時性,保證AGV車SOC值的準(zhǔn)確性在5%以內(nèi),有效解決了工程中AGV車SOC值估算不準(zhǔn)確的問題。

    [1] 朱從民, 黃玉美, 上官望義,等. AGV多傳感器導(dǎo)航系統(tǒng)研究[J]. 儀器儀表學(xué)報, 2008, 29(11):2419-2423.

    [2] 朱玉玉, 劉福兵, 李朋飛. AGV車用鋰離子電池管理系統(tǒng)設(shè)計[J]. 電源技術(shù), 2015, 39(5):917-920.

    [3] 劉勝永, 李 源, 趙振森,等. 基于DSP的電動汽車鋰電池荷電狀態(tài)估算的研究與實現(xiàn)[J]. 計算機(jī)測量與控制, 2015, 23(10):3533-3535.

    [4] 吳鐵洲. HEV鋰離子電池組管理關(guān)鍵技術(shù)研究[D].武漢:華中科技大學(xué), 2010.

    [5] He H, Xiong R, Zhang X, et al. State-of-charge estimation of the lithium-ion battery using an adaptive extended Kalman filter based on an improved the venin model[J]. IEEE Transactions on Vehicular Technology, 2011, 60(4):1461-1469.

    [6] Aung H, Soon Low K, Ting Goh S. State-of-Charge Estimation of Lithium-Ion Battery Using Square Root Spherical Unscented Kalman Filter (Sqrt-UKFST) in Nanosatellite[J]. IEEE Transactions on Power Electronics, 2014, 30(9):4774-4783.

    [7] Chen Z, Fu Y, Mi C C. State of charge estimation of lithium-ion batteries in electric drive vehicles using extended kalman filtering[J]. IEEE Transactions on Vehicular Technology, 2013, 62(3):1020-1030.

    Research on SOC Estimation Algorithm Applied to AGV vehicle

    Wu Tiezhou2, Zhang Min2, Zeng Yishi2, Xiong Jinlong3

    (1.Hubei University of Technology, Wuhan 430068,China;2.State Grid Hubei Electric Power Company Suizhou Power Supply Company, Suizhou 441300,China;3.State Grid Hubei Electric Power Company Maintenance Company, Wuhan 430051,China)

    The working condition of autonomous guided vehicle (AGV) is special, The error of battery SOC is estimated by the current integration method, and also there is a cumulative error. The SOC estimation accuracy can be improved by using the EKF method. Aiming at the special working conditions of the AGV vehicle, the filtering gain of the EKF method is improved for the dynamic adjustment of the filter gain. Effectively improve the tracking performance of the EKF method. The experimental results show that using the EKF method to estimate the SOC of the AGV car battery is higher, The tracking effect of the estimation process is improved by using the filter gain of the dynamic correction, which solves the problem of inaccurate estimation of the SOC of the battery in the special condition of the AGV vehicle.

    battery remaining capacity;autonomous guided vehicle;kalman filter

    2017-02-14;

    2017-03-05。

    國家自然科學(xué)基金(51677058)。

    吳鐵洲(1966-),男,湖北天門人,博士,教授,主要從事儲能控制及能量管理方向的研究。

    1671-4598(2017)08-0166-04

    10.16526/j.cnki.11-4762/tp.2017.08.043

    TM912

    A

    猜你喜歡
    積分法初值卡爾曼濾波
    具非定常數(shù)初值的全變差方程解的漸近性
    一種適用于平動點周期軌道初值計算的簡化路徑搜索修正法
    三維擬線性波方程的小初值光滑解
    基于遞推更新卡爾曼濾波的磁偶極子目標(biāo)跟蹤
    巧用第一類換元法求解不定積分
    基于模糊卡爾曼濾波算法的動力電池SOC估計
    基于擴(kuò)展卡爾曼濾波的PMSM無位置傳感器控制
    隨機(jī)結(jié)構(gòu)地震激勵下的可靠度Gauss-legendre積分法
    基于積分法的軸對稱拉深成形凸緣區(qū)應(yīng)力、應(yīng)變數(shù)值解
    探討不定積分分部積分法
    河南科技(2014年15期)2014-02-27 14:12:50
    a级毛片在线看网站| 又大又爽又粗| 狠狠婷婷综合久久久久久88av| 色网站视频免费| 久久精品人人爽人人爽视色| 亚洲av日韩精品久久久久久密 | 亚洲图色成人| 丁香六月天网| 免费少妇av软件| 国产成人午夜福利电影在线观看| 久久精品国产亚洲av涩爱| 波野结衣二区三区在线| 欧美日韩一区二区视频在线观看视频在线| 精品福利永久在线观看| 亚洲精品自拍成人| 久久99精品国语久久久| 亚洲激情五月婷婷啪啪| 一区二区三区精品91| 男人操女人黄网站| 极品人妻少妇av视频| 欧美日韩视频高清一区二区三区二| av国产久精品久网站免费入址| 在现免费观看毛片| 久久久久精品性色| 久久久久久久精品精品| 精品久久久精品久久久| 侵犯人妻中文字幕一二三四区| 伊人久久大香线蕉亚洲五| 久久久久久免费高清国产稀缺| av在线老鸭窝| 51午夜福利影视在线观看| 99久久精品国产亚洲精品| 美女福利国产在线| 在线观看人妻少妇| 亚洲精品aⅴ在线观看| 久久久久视频综合| 丁香六月欧美| 精品人妻一区二区三区麻豆| 在线看a的网站| 成年女人毛片免费观看观看9 | 国产一区二区三区av在线| av福利片在线| 18禁裸乳无遮挡动漫免费视频| 夫妻午夜视频| 我的亚洲天堂| 国产成人免费无遮挡视频| 美国免费a级毛片| 亚洲精品第二区| 999久久久国产精品视频| 国产精品蜜桃在线观看| 亚洲成人免费av在线播放| 考比视频在线观看| videosex国产| 国产精品久久久久久人妻精品电影 | 日韩电影二区| 18禁国产床啪视频网站| 亚洲av日韩精品久久久久久密 | 大话2 男鬼变身卡| av一本久久久久| 最近的中文字幕免费完整| 日韩欧美一区视频在线观看| 在线观看免费日韩欧美大片| 在线观看国产h片| 99热全是精品| 免费看不卡的av| 亚洲成人av在线免费| 男女下面插进去视频免费观看| 久久性视频一级片| av片东京热男人的天堂| 国产成人欧美在线观看 | 亚洲精品自拍成人| 国产熟女午夜一区二区三区| 一本—道久久a久久精品蜜桃钙片| 激情五月婷婷亚洲| 亚洲欧美清纯卡通| 亚洲四区av| 久久人人爽人人片av| www.自偷自拍.com| 搡老岳熟女国产| 国产午夜精品一二区理论片| 久久女婷五月综合色啪小说| 国产一区亚洲一区在线观看| 亚洲欧美清纯卡通| 日韩电影二区| 亚洲美女黄色视频免费看| 国产成人av激情在线播放| 国产无遮挡羞羞视频在线观看| 国产在线视频一区二区| 国产熟女午夜一区二区三区| 亚洲综合精品二区| 另类亚洲欧美激情| 最黄视频免费看| 少妇猛男粗大的猛烈进出视频| 亚洲图色成人| 伊人亚洲综合成人网| 男女国产视频网站| 亚洲熟女毛片儿| 欧美黑人欧美精品刺激| 最新在线观看一区二区三区 | 人人澡人人妻人| 大片免费播放器 马上看| 一边摸一边抽搐一进一出视频| 欧美日韩亚洲高清精品| 欧美精品高潮呻吟av久久| 男女床上黄色一级片免费看| 久久国产精品大桥未久av| 国产麻豆69| 一级毛片我不卡| 国产精品久久久av美女十八| 亚洲精品久久久久久婷婷小说| 人人澡人人妻人| 久久久久久久国产电影| 久久精品国产亚洲av高清一级| 2018国产大陆天天弄谢| 久久久精品区二区三区| 成人亚洲欧美一区二区av| 美女主播在线视频| 国产精品亚洲av一区麻豆 | 久久久精品区二区三区| h视频一区二区三区| 叶爱在线成人免费视频播放| 欧美日韩福利视频一区二区| 人人妻人人添人人爽欧美一区卜| 下体分泌物呈黄色| 亚洲成人免费av在线播放| 国产在线免费精品| 一级片免费观看大全| 日韩一卡2卡3卡4卡2021年| 久久99热这里只频精品6学生| a 毛片基地| 乱人伦中国视频| 女性生殖器流出的白浆| 99精品久久久久人妻精品| av国产精品久久久久影院| 999精品在线视频| a级片在线免费高清观看视频| 欧美在线一区亚洲| 国产一区二区在线观看av| 精品一区二区三区av网在线观看 | 日韩欧美精品免费久久| 国产精品免费视频内射| 韩国精品一区二区三区| 另类亚洲欧美激情| 天天躁夜夜躁狠狠久久av| 一级毛片 在线播放| 久久综合国产亚洲精品| 免费日韩欧美在线观看| 国产精品二区激情视频| 女人精品久久久久毛片| 天堂俺去俺来也www色官网| 国产在线免费精品| 男女国产视频网站| 黑人巨大精品欧美一区二区蜜桃| 国产伦人伦偷精品视频| 国产av国产精品国产| videos熟女内射| 中文天堂在线官网| 国产免费视频播放在线视频| 日韩 欧美 亚洲 中文字幕| 大香蕉久久成人网| 婷婷色综合www| 另类亚洲欧美激情| 最近最新中文字幕免费大全7| 最近的中文字幕免费完整| 成年人免费黄色播放视频| 女人精品久久久久毛片| 青春草国产在线视频| 国产精品.久久久| 精品久久蜜臀av无| 黄色视频在线播放观看不卡| 欧美在线黄色| 老司机影院毛片| 91aial.com中文字幕在线观看| 久久久精品94久久精品| 一边摸一边抽搐一进一出视频| 波多野结衣av一区二区av| 波野结衣二区三区在线| av福利片在线| av在线播放精品| 婷婷色av中文字幕| 国产精品久久久av美女十八| 人妻 亚洲 视频| 欧美日韩综合久久久久久| 91国产中文字幕| 色吧在线观看| 亚洲国产精品一区三区| 久热爱精品视频在线9| 嫩草影视91久久| 宅男免费午夜| 午夜福利网站1000一区二区三区| 99精品久久久久人妻精品| 久久ye,这里只有精品| 一本大道久久a久久精品| 亚洲国产欧美日韩在线播放| 两个人免费观看高清视频| 亚洲成色77777| 亚洲一卡2卡3卡4卡5卡精品中文| 制服人妻中文乱码| 老司机深夜福利视频在线观看 | 久久毛片免费看一区二区三区| 老鸭窝网址在线观看| 中文字幕最新亚洲高清| 国产黄色视频一区二区在线观看| 91成人精品电影| 国产免费现黄频在线看| 哪个播放器可以免费观看大片| 搡老乐熟女国产| 少妇 在线观看| 午夜免费观看性视频| 大话2 男鬼变身卡| 亚洲综合色网址| 国产精品嫩草影院av在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 国产有黄有色有爽视频| 精品久久蜜臀av无| 色婷婷久久久亚洲欧美| 九九爱精品视频在线观看| 成人免费观看视频高清| 免费在线观看黄色视频的| 亚洲,一卡二卡三卡| 国产精品久久久久久久久免| 亚洲,欧美,日韩| 国产 一区精品| 校园人妻丝袜中文字幕| 国产成人欧美| 国产一区二区三区综合在线观看| 一区在线观看完整版| 亚洲成人手机| 涩涩av久久男人的天堂| 最近中文字幕2019免费版| 99久久精品国产亚洲精品| 中文字幕精品免费在线观看视频| 日本91视频免费播放| 亚洲国产精品一区二区三区在线| 精品午夜福利在线看| 久久精品亚洲熟妇少妇任你| 国产xxxxx性猛交| 久久精品国产亚洲av高清一级| 亚洲av福利一区| a级毛片黄视频| 精品一品国产午夜福利视频| 伊人亚洲综合成人网| 成人三级做爰电影| 五月天丁香电影| 午夜免费观看性视频| 汤姆久久久久久久影院中文字幕| 精品一区二区三区av网在线观看 | 午夜精品国产一区二区电影| 19禁男女啪啪无遮挡网站| 最黄视频免费看| 精品卡一卡二卡四卡免费| 乱人伦中国视频| 国产有黄有色有爽视频| 777米奇影视久久| 老司机在亚洲福利影院| 精品一品国产午夜福利视频| 麻豆av在线久日| 国产精品一区二区在线不卡| 亚洲成人一二三区av| 在线天堂中文资源库| 久久久久精品久久久久真实原创| 午夜日韩欧美国产| videos熟女内射| 国产成人精品无人区| √禁漫天堂资源中文www| 欧美成人午夜精品| 欧美日本中文国产一区发布| 国产极品粉嫩免费观看在线| 国产免费一区二区三区四区乱码| 老熟女久久久| 国产熟女欧美一区二区| 啦啦啦中文免费视频观看日本| 一区二区日韩欧美中文字幕| 国产精品人妻久久久影院| 女的被弄到高潮叫床怎么办| 久久精品亚洲av国产电影网| 黄色视频在线播放观看不卡| 欧美精品一区二区大全| 女性被躁到高潮视频| 久久久国产精品麻豆| 国产一区二区 视频在线| 美女中出高潮动态图| 日韩一卡2卡3卡4卡2021年| 悠悠久久av| 18在线观看网站| 亚洲一区中文字幕在线| 久久99一区二区三区| 亚洲欧美一区二区三区久久| 国产精品久久久人人做人人爽| 如日韩欧美国产精品一区二区三区| 亚洲av电影在线观看一区二区三区| 国产精品免费视频内射| 欧美xxⅹ黑人| 一本久久精品| 国产av精品麻豆| 日日爽夜夜爽网站| 岛国毛片在线播放| 亚洲精品久久午夜乱码| 欧美黑人欧美精品刺激| 国产精品女同一区二区软件| 亚洲精华国产精华液的使用体验| 777米奇影视久久| 亚洲少妇的诱惑av| 国产麻豆69| 国产一级毛片在线| 免费人妻精品一区二区三区视频| 免费黄色在线免费观看| 一级片免费观看大全| 国产精品熟女久久久久浪| 欧美日韩综合久久久久久| 久久 成人 亚洲| 亚洲精品美女久久久久99蜜臀 | 国产精品一二三区在线看| 老鸭窝网址在线观看| 免费在线观看视频国产中文字幕亚洲 | 天堂俺去俺来也www色官网| 亚洲av欧美aⅴ国产| 99九九在线精品视频| 啦啦啦 在线观看视频| 制服人妻中文乱码| 免费在线观看完整版高清| 一级毛片电影观看| 国产精品久久久人人做人人爽| 亚洲综合色网址| 一个人免费看片子| 十分钟在线观看高清视频www| 日韩制服丝袜自拍偷拍| 精品一区二区免费观看| videos熟女内射| 精品一区二区三区四区五区乱码 | 国产亚洲av片在线观看秒播厂| 国产黄色免费在线视频| 我的亚洲天堂| 亚洲av在线观看美女高潮| 男的添女的下面高潮视频| 观看美女的网站| 在线观看人妻少妇| 国产精品 国内视频| 啦啦啦在线免费观看视频4| 无遮挡黄片免费观看| 国产精品国产三级专区第一集| 丝袜美足系列| 国产淫语在线视频| xxx大片免费视频| 最近的中文字幕免费完整| 日韩精品免费视频一区二区三区| 18在线观看网站| 热re99久久国产66热| 人妻 亚洲 视频| bbb黄色大片| 欧美日韩av久久| 黄色一级大片看看| 精品久久蜜臀av无| 热re99久久精品国产66热6| 国产熟女欧美一区二区| 一级毛片电影观看| 99久久人妻综合| 久久久亚洲精品成人影院| 五月开心婷婷网| 色婷婷久久久亚洲欧美| 久久天堂一区二区三区四区| 操出白浆在线播放| 欧美日韩精品网址| 久久国产精品大桥未久av| 亚洲国产中文字幕在线视频| 欧美日韩亚洲高清精品| 亚洲欧美精品自产自拍| 中文字幕色久视频| 日本欧美国产在线视频| 高清视频免费观看一区二区| 亚洲精品日本国产第一区| 丝袜在线中文字幕| 日本黄色日本黄色录像| 亚洲色图 男人天堂 中文字幕| 日韩av不卡免费在线播放| 国产精品欧美亚洲77777| 我要看黄色一级片免费的| 亚洲成人一二三区av| 不卡视频在线观看欧美| 日韩一卡2卡3卡4卡2021年| 亚洲国产最新在线播放| 亚洲国产欧美日韩在线播放| 亚洲精品美女久久av网站| 国产精品无大码| 卡戴珊不雅视频在线播放| 自线自在国产av| 日本黄色日本黄色录像| 午夜老司机福利片| 午夜日韩欧美国产| 国产在线一区二区三区精| 色网站视频免费| 日日撸夜夜添| av不卡在线播放| 日韩电影二区| 亚洲av国产av综合av卡| 91老司机精品| 久久久久久人人人人人| 美女高潮到喷水免费观看| 精品国产一区二区久久| 日日摸夜夜添夜夜爱| 国产欧美日韩一区二区三区在线| 成人免费观看视频高清| 99久久精品国产亚洲精品| 午夜福利影视在线免费观看| 亚洲av在线观看美女高潮| 色网站视频免费| 在线观看三级黄色| 亚洲欧美中文字幕日韩二区| 亚洲精品久久成人aⅴ小说| 18在线观看网站| 在线观看免费高清a一片| 亚洲av成人不卡在线观看播放网 | 国产精品电影一区二区三区| 久久久久久免费高清国产稀缺| 亚洲av日韩精品久久久久久密| 久久久久久亚洲精品国产蜜桃av| 亚洲专区字幕在线| 99re在线观看精品视频| 亚洲人成电影观看| 真人做人爱边吃奶动态| 午夜福利欧美成人| 久99久视频精品免费| 纯流量卡能插随身wifi吗| 国产精品爽爽va在线观看网站 | 女生性感内裤真人,穿戴方法视频| 精品一区二区三区四区五区乱码| 日韩精品中文字幕看吧| 99热只有精品国产| 午夜老司机福利片| 757午夜福利合集在线观看| 久久精品国产清高在天天线| 亚洲一卡2卡3卡4卡5卡精品中文| 国产亚洲精品一区二区www| 婷婷六月久久综合丁香| 在线观看日韩欧美| 欧美性长视频在线观看| 色综合亚洲欧美另类图片| 欧美日韩福利视频一区二区| 久久精品亚洲精品国产色婷小说| 国产99久久九九免费精品| 自线自在国产av| 两个人视频免费观看高清| 久久午夜综合久久蜜桃| 日韩免费av在线播放| 亚洲av熟女| 久久久久久久久免费视频了| av天堂在线播放| 欧美成人一区二区免费高清观看 | 日韩欧美在线二视频| 国产一级毛片七仙女欲春2 | 亚洲成av片中文字幕在线观看| 一进一出抽搐动态| 不卡av一区二区三区| 91麻豆精品激情在线观看国产| 久久香蕉激情| 两性午夜刺激爽爽歪歪视频在线观看 | 午夜视频精品福利| 99国产综合亚洲精品| 日韩欧美国产在线观看| 一进一出抽搐gif免费好疼| 亚洲精品中文字幕一二三四区| 久久热在线av| 亚洲美女黄片视频| 18美女黄网站色大片免费观看| 国产成人免费无遮挡视频| 美女大奶头视频| 国产伦人伦偷精品视频| 男人操女人黄网站| 女人爽到高潮嗷嗷叫在线视频| 精品一品国产午夜福利视频| 91国产中文字幕| 丝袜美足系列| 男女做爰动态图高潮gif福利片 | 久久久久久亚洲精品国产蜜桃av| 九色亚洲精品在线播放| 国产一区二区三区综合在线观看| 最新在线观看一区二区三区| 一级黄色大片毛片| 女生性感内裤真人,穿戴方法视频| 成年女人毛片免费观看观看9| 欧美日韩黄片免| 一夜夜www| 看免费av毛片| 婷婷六月久久综合丁香| 日韩一卡2卡3卡4卡2021年| 亚洲最大成人中文| 99在线人妻在线中文字幕| 午夜福利成人在线免费观看| 精品久久久久久,| 午夜精品久久久久久毛片777| 欧美 亚洲 国产 日韩一| 岛国视频午夜一区免费看| 免费在线观看亚洲国产| 女性被躁到高潮视频| 国产精品日韩av在线免费观看 | 亚洲三区欧美一区| 天堂动漫精品| 亚洲精品国产一区二区精华液| 亚洲国产欧美日韩在线播放| 1024视频免费在线观看| 亚洲色图综合在线观看| 国产一区在线观看成人免费| 国产欧美日韩一区二区精品| 国产乱人伦免费视频| av在线天堂中文字幕| 人妻丰满熟妇av一区二区三区| 男人舔女人下体高潮全视频| 精品一区二区三区四区五区乱码| 国产亚洲av高清不卡| 欧美成人性av电影在线观看| 日日摸夜夜添夜夜添小说| 18禁国产床啪视频网站| 狂野欧美激情性xxxx| 亚洲av成人不卡在线观看播放网| 久久性视频一级片| 亚洲国产精品合色在线| 啦啦啦 在线观看视频| 9色porny在线观看| avwww免费| 国产亚洲欧美98| 午夜精品国产一区二区电影| 一边摸一边抽搐一进一小说| 欧美 亚洲 国产 日韩一| 国产成人精品在线电影| 欧美久久黑人一区二区| 日本欧美视频一区| 无人区码免费观看不卡| 亚洲国产欧美日韩在线播放| 日韩高清综合在线| 777久久人妻少妇嫩草av网站| 大陆偷拍与自拍| 窝窝影院91人妻| 亚洲久久久国产精品| 国产高清videossex| 国产精品久久久av美女十八| 女警被强在线播放| 999久久久精品免费观看国产| av超薄肉色丝袜交足视频| 欧美绝顶高潮抽搐喷水| 久久精品影院6| 两个人看的免费小视频| 国产高清激情床上av| 99久久综合精品五月天人人| 亚洲av成人不卡在线观看播放网| 一区在线观看完整版| 亚洲欧美日韩另类电影网站| 夜夜爽天天搞| www日本在线高清视频| 日韩欧美免费精品| 男男h啪啪无遮挡| 久久香蕉国产精品| 人人澡人人妻人| 免费av毛片视频| 久久久久精品国产欧美久久久| 国产高清激情床上av| 亚洲第一欧美日韩一区二区三区| 成人国产一区最新在线观看| 国产精华一区二区三区| 国产精品99久久99久久久不卡| 欧美成人性av电影在线观看| 中文字幕精品免费在线观看视频| 成在线人永久免费视频| 香蕉久久夜色| 久久中文字幕人妻熟女| 国产亚洲精品综合一区在线观看 | 欧洲精品卡2卡3卡4卡5卡区| 国产又色又爽无遮挡免费看| 亚洲 欧美一区二区三区| 国产午夜精品久久久久久| 麻豆av在线久日| 成人永久免费在线观看视频| 日韩欧美免费精品| 在线十欧美十亚洲十日本专区| 亚洲人成77777在线视频| 麻豆一二三区av精品| 国产三级黄色录像| 欧美最黄视频在线播放免费| 久久精品影院6| 真人做人爱边吃奶动态| 露出奶头的视频| 99久久99久久久精品蜜桃| 丝袜美腿诱惑在线| 国产伦一二天堂av在线观看| 欧美日韩乱码在线| 国产伦一二天堂av在线观看| av视频在线观看入口| 成年版毛片免费区| 免费人成视频x8x8入口观看| 少妇熟女aⅴ在线视频| 免费观看人在逋| 两性午夜刺激爽爽歪歪视频在线观看 | cao死你这个sao货| 精品福利观看| 母亲3免费完整高清在线观看| 岛国在线观看网站| 久久久久久久久久久久大奶| 97人妻天天添夜夜摸| 啪啪无遮挡十八禁网站| 精品福利观看| 视频在线观看一区二区三区| 亚洲熟女毛片儿| 制服诱惑二区| 午夜精品国产一区二区电影| 三级毛片av免费| 韩国精品一区二区三区| aaaaa片日本免费| 亚洲精品一区av在线观看| 成人国语在线视频| 欧美av亚洲av综合av国产av| 亚洲欧美精品综合久久99| 日本 av在线| 高清在线国产一区| 欧美老熟妇乱子伦牲交| 夜夜躁狠狠躁天天躁| 女生性感内裤真人,穿戴方法视频|