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

    用新陳代謝極限學(xué)習(xí)機(jī)實(shí)現(xiàn)電池健康狀態(tài)估算*

    2021-02-02 08:13:20王惠民李熠婧潘海鴻
    汽車工程 2021年1期
    關(guān)鍵詞:模型

    陳 琳,王惠民,李熠婧,張 沫,黃 江,潘海鴻

    (1.廣西大學(xué)機(jī)械工程學(xué)院,南寧 530004;2.廣西電化學(xué)能源材料重點(diǎn)實(shí)驗(yàn)室(廣西大學(xué)可再生能源材料協(xié)同創(chuàng)新中心),南寧 530004)

    前言

    鋰離子電池因具有能量密度大和循環(huán)壽命長(zhǎng)等優(yōu)勢(shì)作為主要儲(chǔ)能元件被廣泛應(yīng)用于電動(dòng)汽車和可再生能源系統(tǒng)等各個(gè)工業(yè)領(lǐng)域[1-3]。由于鋰離子電池的性能狀態(tài)直接影響到整個(gè)系統(tǒng)的穩(wěn)定性和可靠性,對(duì)其健康狀態(tài)(state of health,SOH)進(jìn)行估算十分必要。SOH主要反映電池當(dāng)前性能相對(duì)于其初始狀態(tài)的變化,它決定電池的維護(hù)和更換時(shí)間,可為電池提供決策性保護(hù)[3-4]。因此,準(zhǔn)確估算SOH 對(duì)于保障儲(chǔ)能系統(tǒng)性能、降低成本和提高安全性等方面具有重要意義[5]。

    近年來(lái)國(guó)內(nèi)外學(xué)者已在SOH 估算領(lǐng)域開(kāi)展了大量研究,主要方法可大致分為3 類:機(jī)理分析法[6-7]、基于濾波法[1,8]和數(shù)據(jù)驅(qū)動(dòng)法[9-14]。其中,數(shù)據(jù)驅(qū)動(dòng)方法因其無(wú)需進(jìn)行深入的機(jī)理研究,僅靠挖掘歷史數(shù)據(jù)的潛在規(guī)律就可以描述復(fù)雜的電池退化過(guò)程而受到廣泛關(guān)注[1]。這類方法主要包括支持向量 機(jī)(support vector machine,SVM)[9-10]、神 經(jīng) 網(wǎng)絡(luò)[11-12]和高斯過(guò)程回歸(gauss progress regression,GPR)[13-14]等。文獻(xiàn)[9]中在考慮多負(fù)載集合和容量測(cè)試耦合性的情況下,將SVM 用于電池容量變化的回歸預(yù)測(cè)來(lái)實(shí)現(xiàn)鋰離子電池的SOH 估算。文獻(xiàn)[10]中引入粒子群優(yōu)化算法對(duì)SVM 的超參數(shù)進(jìn)行優(yōu)化,減少使用SVM 估算SOH 時(shí)帶來(lái)的泛化誤差。文獻(xiàn)[11]中使用電池的動(dòng)態(tài)歷史數(shù)據(jù)分布訓(xùn)練其所構(gòu)建的反向傳播神經(jīng)網(wǎng)絡(luò)框架,結(jié)果表明該框架可用于多工況下的電池SOH 估算。文獻(xiàn)[13]中提出了一種利用多電池?cái)?shù)據(jù)的多輸出GPR 模型,該模型充分利用GPR 的超參數(shù)自適應(yīng)調(diào)節(jié)的優(yōu)點(diǎn),顯著提高SOH 預(yù)測(cè)性能。這些數(shù)據(jù)驅(qū)動(dòng)方法在非線性非高斯的SOH 估計(jì)中具有良好的適應(yīng)性,并且在輸入充足歷史數(shù)據(jù)的情況下可以達(dá)到很高的精度[1,5,12]。但在輸入數(shù)據(jù)量較少的情況下,數(shù)據(jù)驅(qū)動(dòng)算法的泛化性能會(huì)下降[5],影響估算效果。因此,如何在小數(shù)據(jù)樣本情況下對(duì)SOH 進(jìn)行準(zhǔn)確估算成為數(shù)據(jù)驅(qū)動(dòng)方法的研究難點(diǎn)。

    極限學(xué)習(xí)機(jī)(extreme learning machine,ELM)具有結(jié)構(gòu)簡(jiǎn)單和易于調(diào)整參數(shù)等優(yōu)點(diǎn)[15],已被應(yīng)用在其他領(lǐng)域的回歸預(yù)測(cè)案例中。與SVM 和BPNN 等相比,其訓(xùn)練速度更快且不易陷入局部最優(yōu),泛化能力十分優(yōu)秀。因此,基于在線提取的退化特征和容量退化量,引入ELM 來(lái)利用其泛化性優(yōu)勢(shì)構(gòu)建一種描述電池復(fù)雜退化過(guò)程的退化狀態(tài)模型,并結(jié)合新陳代謝機(jī)制對(duì)小數(shù)據(jù)樣本情況下的電池SOH 進(jìn)行估算。此外,采用正極材料分別為鈷酸鋰和鎳鈷鋁三元材料的電池在動(dòng)態(tài)工況下的測(cè)試數(shù)據(jù)對(duì)所提方法的SOH估算效果進(jìn)行驗(yàn)證。

    1 基本方法

    1.1 基于退化特征提取的退化狀態(tài)模型構(gòu)建

    電池的退化是典型的非線性時(shí)變過(guò)程,每個(gè)循環(huán)的退化量與之前循環(huán)的參數(shù)變化都密切相關(guān),而這種退化關(guān)系的詳細(xì)電化學(xué)機(jī)理十分復(fù)雜,很難通過(guò)具體的數(shù)學(xué)公式進(jìn)行量化。為此提出構(gòu)建一種以容量退化為狀態(tài)變量的電池退化狀態(tài)模型,將相鄰循環(huán)之間容量退化量的關(guān)聯(lián)性抽象地表示為

    式中:ΔCk為第k 個(gè)循環(huán)的容量退化量;f(·)為根據(jù)第k-1 個(gè)循環(huán)的容量退化量ΔCk-1建立的描述電池退化機(jī)理的函數(shù)。

    當(dāng)電池容量退化發(fā)生較大變化時(shí),式(1)中的模型可能無(wú)法有效地跟蹤電池容量的退化,為此提出引入其他退化特征對(duì)該公式進(jìn)行修正。因此,根據(jù)電池的放電電流和電壓數(shù)據(jù)建立電池的1 階RC 模型,使用遞推最小二乘算法實(shí)時(shí)辨識(shí)出電池的歐姆內(nèi)阻R0和極化內(nèi)阻R1,并將這兩種內(nèi)阻用于表征電池退化過(guò)程[16],結(jié)合兩種內(nèi)阻和容量退化量構(gòu)建老化模型。內(nèi)阻辨識(shí)的主要步驟如下。

    建立電池的1 階RC 模型(圖1)。圖中:UOC表示電池的開(kāi)路電壓;R1和C1表示電池內(nèi)部的極化內(nèi)阻和極化電容,分別代表電池內(nèi)部電化學(xué)反應(yīng)的極化特性和擴(kuò)散效應(yīng);R0是電池的歐姆內(nèi)阻,主要包括電池隔膜電阻和其他接觸電阻。

    圖1 電池的1階RC模型

    由Kirchhoff 定律可以得到1 階RC 模型的電路表達(dá)式:

    式中:Ut為電池端電壓;U1為RC 網(wǎng)絡(luò)兩端的電壓;IL為負(fù)載電流。

    將式(3)代入式(2)并進(jìn)行z 變換,得到式(2)的離散表達(dá)式:

    式中:k為采樣時(shí)刻;Uoc,k、Ut,k和IL,k分別為電池在k時(shí)刻的開(kāi)路電壓、端電壓和負(fù)載電流;Ts為采樣周期。

    式(4)中的參數(shù)Ut和IL都可以通過(guò)儀器直接測(cè)量獲取,因此矩陣Φk完全已知。之后可通過(guò)式(5)對(duì)矩陣Θk進(jìn)行遞推求解:

    式中:λ 為遺忘因子,本文中取值0.95;Pk和Kk分別為k 時(shí)刻的協(xié)方差矩陣和增益矩陣。根據(jù)式(5)計(jì)算得到Θk之后,即可通過(guò)式(4)求解出R0、R1、C1的值。

    文獻(xiàn)[1]中指出,電動(dòng)汽車動(dòng)力電池的實(shí)際工作荷電狀態(tài)(SOC)區(qū)間主要為[25%,95%];前期研究[17]中發(fā)現(xiàn),當(dāng)電池在其SOC 處于[30%,80%]區(qū)間內(nèi)充放電時(shí),實(shí)時(shí)辨識(shí)出的R0和R1基本保持不變。因此,選擇電池在每個(gè)循環(huán)SOC 的[30%,80%]區(qū)間內(nèi)的平均歐姆內(nèi)阻和平均極化內(nèi)阻來(lái)代表實(shí)際應(yīng)用環(huán)境中的真實(shí)電池內(nèi)阻。平均歐姆內(nèi)阻和平均極化內(nèi)阻的計(jì)算方法見(jiàn)式(6)和式(7)。

    然而,溫度變化對(duì)電池內(nèi)阻存在明顯影響,不同溫度下所辨識(shí)出內(nèi)阻的數(shù)值差異較大,將其作為特征提取會(huì)增加算法復(fù)雜度和計(jì)算處理過(guò)程,影響模型的估算效率和精度。故選擇將歐姆內(nèi)阻和極化內(nèi)阻的增量作為電池的退化特征,計(jì)算方式見(jiàn)式(8)和式(9)。

    式中g(shù)(·)為根據(jù)第k-1個(gè)循環(huán)的容量退化量ΔCk-1和退化特征ΔR0k-1、ΔR1k-1建立的描述電池退化機(jī)理的函數(shù)。

    1.2 ELM基本算法原理

    由于電池內(nèi)部電化學(xué)反應(yīng)過(guò)程十分復(fù)雜,導(dǎo)致式(10)所給出的退化狀態(tài)模型很難得到準(zhǔn)確的解析表達(dá)式,為此提出使用ELM 算法對(duì)退化狀態(tài)模型進(jìn)行擬合。ELM是一種對(duì)輸入權(quán)值和隱含層節(jié)點(diǎn)偏差進(jìn)行隨機(jī)賦值的單隱層前饋神經(jīng)網(wǎng)絡(luò)[15],其學(xué)習(xí)效率高且泛化性能好。算法基本過(guò)程如下。

    式中:j=1,2,…,N;ωi=[ωi1,ωi2,…,ωin]T為連接第i個(gè)隱含層節(jié)點(diǎn)和輸入節(jié)點(diǎn)的輸入權(quán)值;bi為第i 個(gè)隱含層節(jié)點(diǎn)的偏差;βi為連接第i 個(gè)隱含層節(jié)點(diǎn)和輸出節(jié)點(diǎn)的輸出權(quán)值。

    式(11)可以被簡(jiǎn)寫(xiě)成線性方程組:

    其中:

    式中ωi和bi為隨機(jī)賦值。只需再求出βi的值即可完成ELM 訓(xùn)練。βi可通過(guò)求式(12)的最小二乘解來(lái)確定:

    式中H+是H的Moore?Penrose廣義逆矩陣。

    1.3 基于新陳代謝ELM的電池SOH估算架構(gòu)

    采用新陳代謝ELM 進(jìn)行鋰離子電池SOH 估算的主要框架如圖2 所示,主要包括離線構(gòu)建退化狀態(tài)模型和新陳代謝ELM在線估算SOH。

    1.3.1 離線構(gòu)建電池退化狀態(tài)模型

    采用電池全生命周期的測(cè)試數(shù)據(jù)獲取每次老化循環(huán)中的容量退化量序列{ΔC}、歐姆內(nèi)阻增量序列{ΔR0}和極化內(nèi)阻增量序列{ΔR1}。將{ΔC}、{ΔR0}和{ΔR1}作為輸入訓(xùn)練數(shù)據(jù),輸入數(shù)據(jù){ΔC}中每個(gè)循環(huán)對(duì)應(yīng)的下一循環(huán)容量退化量作為目標(biāo)輸出數(shù)據(jù),根據(jù)式(11)~式(14)計(jì)算ELM 的輸出權(quán)值βi,完成基于ELM 的電池退化狀態(tài)模型的訓(xùn)練,即確定式(10)中g(shù)(·)的具體形式,記為gELM(·)。

    圖2 基于新陳代謝ELM的電池SOH估算框架

    1.3.2 新陳代謝ELM在線估算SOH

    實(shí)際工況下的電池狀態(tài)呈現(xiàn)強(qiáng)時(shí)變特性,最新的數(shù)據(jù)才能充分反映電池退化的最新趨勢(shì)和特點(diǎn)。因此,提出利用最新數(shù)據(jù)更新所構(gòu)建的退化狀態(tài)模型對(duì)電池SOH進(jìn)行新陳代謝估算。具體步驟如下。

    (1)假設(shè)已知?dú)v史數(shù)據(jù)的老化循環(huán)為第k-i到第k-1個(gè)循環(huán)(k>i),在線采集這i個(gè)循環(huán)的歷史容量退化數(shù)據(jù)序列{ΔCk-i→k-1},通過(guò)動(dòng)態(tài)工況下實(shí)時(shí)采集的數(shù)據(jù)在線建立電池的1 階RC 模型并辨識(shí)得到這些循環(huán)的歐姆內(nèi)阻增量序列{ΔR0k-i→k-1}和極化內(nèi)阻增量序列{ΔR1k-i→k-1},將其作為退化特征。

    其中:

    (2)將待估算電池的{ΔCk-i→k-1}、{ΔR0k-i→k-1}和{ΔR1k-i→k-1}輸入1.3.1 節(jié)中所構(gòu)建的電池退化狀態(tài)模型gELM(·),模型輸出待估算電池的容量退化量估算值序列該步驟可以簡(jiǎn)單描述為

    (3)根據(jù)式(16)計(jì)算第k個(gè)循環(huán)的SOH估算值:

    式中Cinitial為待估算電池的初始容量。

    將式(18)構(gòu)建的新輸入序列輸入電池退化狀態(tài)模型gELM(·),即返回步驟(2);不斷循環(huán)步驟(2)~步驟(4)來(lái)更新估算,直至容量退化量估算值達(dá)到設(shè)定的截止閾值(電池額定容量的20%左右),停止估算。

    2 實(shí)驗(yàn)設(shè)計(jì)

    2.1 實(shí)驗(yàn)平臺(tái)

    實(shí)驗(yàn)平臺(tái)由高低溫實(shí)驗(yàn)箱GDW-100L、可編程直流電源ITECH IT6523D、電子負(fù)載儀ITECH IT8511、上位計(jì)算機(jī)和鋰離子電池組成(圖3)。其中高低溫實(shí)驗(yàn)箱GDW-100L可實(shí)時(shí)記錄實(shí)驗(yàn)過(guò)程中的環(huán)境溫度;可編程直流電源ITECH IT6523D 和電子負(fù)載儀ITECH IT8511 用于對(duì)鋰離子電池進(jìn)行充放電測(cè)試,以1 Hz 頻率實(shí)時(shí)采集電池在充放電過(guò)程中的電流、電壓和累積充放電容量3 種數(shù)據(jù);上位計(jì)算機(jī)用于充放電工況控制和數(shù)據(jù)儲(chǔ)存。受測(cè)電池為商用動(dòng)力鋰離子電池:正極為鎳鈷鋁三元材料(LiNCA)的松下NCR18650B電池1節(jié)(編號(hào)P01);正極為鈷酸鋰(LCO)的三星ICR18650-26F 電池2 節(jié)(分別編號(hào)S17、S30)。具體型號(hào)和參數(shù)如表1所示。

    圖3 實(shí)驗(yàn)設(shè)備及其連接示意圖

    表1 受測(cè)電池的主要參數(shù)

    2.2 實(shí)驗(yàn)流程

    電池的實(shí)驗(yàn)測(cè)試主要由特性測(cè)試和老化循環(huán)測(cè)試組成(圖4),電池在每次測(cè)試前都以標(biāo)準(zhǔn)恒流-恒壓(constant current and constant voltage,CC-CV)充電方式充滿。

    圖4 老化實(shí)驗(yàn)流程圖

    (1)電池特性測(cè)試:包括一次容量測(cè)試和一次動(dòng)態(tài)工況測(cè)試。①容量測(cè)試:將電池用標(biāo)準(zhǔn)CC-CV 充電法充滿,之后以0.2C 倍率進(jìn)行恒流放電,直至電壓下降至放電截止電壓。②動(dòng)態(tài)工況測(cè)試:使用特定動(dòng)態(tài)工況對(duì)電池進(jìn)行循環(huán)放電直至電池的電壓降至放電截止電壓。其中P01 的動(dòng)態(tài)工況測(cè)試所使用工況為25 ℃下的UDDS(urban dynamometer driving schedule)工況,S17和S30的動(dòng)態(tài)工況測(cè)試分別為25和35 ℃下的NEDC(new european driving cycle)工況。特性測(cè)試不同有利于研究不同種類電池在不同工況激勵(lì)下所反映出的老化特性。

    (2)老化循環(huán)測(cè)試:電池老化測(cè)試則是將電池以標(biāo)準(zhǔn)CC-CV 充電法滿充后,以特定老化循環(huán)工況對(duì)電池進(jìn)行循環(huán)放電直至電壓降至放電截止電壓,該充放電步驟重復(fù)10 次。由于高溫、大倍率、100%放電深度、間隔放電等要素有利于加快電池老化[18],為加速電池老化過(guò)程,早期進(jìn)行的S30 老化實(shí)驗(yàn)在35 ℃下使用較大倍率的脈沖電流[19]作為老化循環(huán)測(cè)試工況;后續(xù)考慮到實(shí)際環(huán)境中電池放電工況并非脈沖放電能夠代替,故在S17和P01的老化循環(huán)測(cè)試中加入了UDDS工況,并提高實(shí)驗(yàn)溫度至45 ℃,在進(jìn)一步加速電池老化的同時(shí)充分體現(xiàn)動(dòng)態(tài)工況對(duì)電池老化的影響。不同的老化循環(huán)測(cè)試有利于探究?jī)煞N電池在不同動(dòng)態(tài)工況下進(jìn)行老化循環(huán)產(chǎn)生的老化差異。3節(jié)電池的實(shí)驗(yàn)流程對(duì)比如表2所示。

    表2 受測(cè)電池的實(shí)驗(yàn)流程對(duì)比

    重復(fù)進(jìn)行電池特性測(cè)試和老化循環(huán)測(cè)試直至電池的容量退化量達(dá)到初始容量的20%。所有實(shí)驗(yàn)步驟均在高低溫實(shí)驗(yàn)箱內(nèi)進(jìn)行,并用上位計(jì)算機(jī)實(shí)時(shí)記錄。

    3 實(shí)驗(yàn)結(jié)果與分析

    3.1 退化特征與電池性能退化

    根據(jù)老化實(shí)驗(yàn)所得容量數(shù)據(jù),不同電池的容量隨著充放電循環(huán)周期的增加都呈現(xiàn)出非線性下降趨勢(shì)(圖5)。對(duì)于正極材料為L(zhǎng)iNCA 的電池P01,其初始容量較高,容量退化速率在整個(gè)老化過(guò)程中一直保持相對(duì)穩(wěn)定,退化軌跡的斜率基本保持不變;而對(duì)于正極材料為L(zhǎng)CO 的兩節(jié)電池S17和S30,它們的退化過(guò)程則存在更多容量波動(dòng),其中S30 在退化后期表現(xiàn)出明顯的加速退化趨勢(shì)。不同電池在退化過(guò)程中表現(xiàn)出的不同特性為SOH 估算帶來(lái)一定困難,因此尋找有效的退化特征來(lái)反映不同電池的退化特性十分必要。

    按照1.1 節(jié)中遞推最小二乘算法辨識(shí)電池的歐姆內(nèi)阻和極化內(nèi)阻。以S30 為例,其辨識(shí)出的歐姆內(nèi)阻在不同充放電循環(huán)中的變化曲線如圖6所示。

    圖5 P01、S17和S30的容量退化圖

    圖6 S30辨識(shí)得到的歐姆內(nèi)阻隨充放電循環(huán)變化曲線

    隨著電池充放電循環(huán)周數(shù)的增加,電池老化程度加深,辨識(shí)得到的歐姆內(nèi)阻逐漸變大,表明歐姆內(nèi)阻可反映電池的老化狀態(tài)。此外,在電池的單次放電過(guò)程中,當(dāng)SOC 處于30%~80%區(qū)間內(nèi)時(shí),其歐姆內(nèi)阻整體變化趨勢(shì)相對(duì)平穩(wěn),將該區(qū)間內(nèi)歐姆內(nèi)阻平均值作為當(dāng)前循環(huán)的標(biāo)準(zhǔn)歐姆內(nèi)阻可有效減少隨機(jī)誤差和異常值的影響。極化內(nèi)阻的變化趨勢(shì)與歐姆內(nèi)阻類似,此處不再贅述。根據(jù)所得內(nèi)阻求其增量構(gòu)建退化特征,分析其與容量退化量之間的相關(guān)程度。圖7 給出P01、S17 和S30 電池的內(nèi)阻增量與容量退化量的關(guān)系。

    圖7 中3 節(jié)電池的兩種內(nèi)阻增量(ΔR0和ΔR1)都和其容量退化量ΔC 存在明顯正相關(guān)。但松下和三星兩種電池的ΔC隨兩種內(nèi)阻增量變化的趨勢(shì)不同:正極為L(zhǎng)iNCA 三元材料的P01 的ΔC 隨ΔR0和ΔR1變化的趨勢(shì)較為相似;而正極材料為L(zhǎng)CO 的S17/S30,其ΔC 與兩種內(nèi)阻增量的線性關(guān)系更明顯且變化趨勢(shì)不同(圖7(b)和圖7(c))。其內(nèi)阻增量變化由于測(cè)試溫度不同而有所差異,在容量退化量達(dá)到0.4 A·h 時(shí),25 ℃下測(cè)試的S43 內(nèi)阻增量略大于35 ℃下測(cè)試的S30 內(nèi)阻增量,這是因?yàn)闇囟葧?huì)影響鋰離子遷移速度,且鋰離子脫出速度和嵌入速度的變化也不同[20]。雖然不同溫度對(duì)內(nèi)阻增量有一定影響,但電池容量退化量和內(nèi)阻增量仍然存在近似線性關(guān)系。

    使用皮爾遜(Pearson)積矩相關(guān)系數(shù)rp和斯皮爾曼(Spearman)等級(jí)相關(guān)系數(shù)rs來(lái)定量地評(píng)判退化特征對(duì)容量退化量的影響程度。rp是量化兩個(gè)變量之間線性相關(guān)程度的參數(shù)標(biāo)準(zhǔn),rs則是用來(lái)度量?jī)蓚€(gè)變量之間單調(diào)相關(guān)性的指標(biāo),按式(19)和式(20)計(jì)算。兩者的取值范圍都在[-1,1]之間,絕對(duì)值越接近1說(shuō)明變量間相關(guān)性越強(qiáng)。

    圖7 不同電池的內(nèi)阻增量與容量退化量的關(guān)系

    式中X和Y為待分析變量序列。

    分別將3 節(jié)電池的內(nèi)阻增量ΔR0和ΔR1代入X,容量退化量代入Y,計(jì)算得到相關(guān)系數(shù)值,見(jiàn)表3。

    表3 退化特征與容量退化量相關(guān)性分析結(jié)果

    3 節(jié)電池的ΔR0和ΔR1與容量退化量的rp和rs都大于0.94,表明以電池內(nèi)阻增加量構(gòu)建的表征電池退化的退化特征和電池容量退化量間具有良好的線性關(guān)系。以在線辨識(shí)得到的鋰離子電池內(nèi)阻增量作為退化特征描述電池性能退化狀態(tài),用兩者之間的近似線性關(guān)系修正容量退化量構(gòu)建的電池退化狀態(tài)模型,最終用于估算電池SOH是合理的。

    3.2 SOH估算結(jié)果分析

    將基于新陳代謝ELM 的電池SOH 估算框架用于不同工況和溫度下的電池SOH 估算,并對(duì)估算精度進(jìn)行評(píng)估。測(cè)試軟件環(huán)境為MATLAB R2017a,ELM 的隱含層節(jié)點(diǎn)數(shù)L=30,激勵(lì)函數(shù)選用Sigmoid函數(shù)。

    在離線構(gòu)建狀態(tài)模型的參考電池和待估算電池型號(hào)及材料不同的情況下對(duì)新陳代謝SOH 估算框架進(jìn)行有效性驗(yàn)證:使用LiNCA 電池P01 離線構(gòu)建電池退化狀態(tài)模型,再用材料和測(cè)試工況都不同的LCO 電池S17、S30 分別進(jìn)行在線SOH 估算驗(yàn)證?;贓LM 的電池退化狀態(tài)模型對(duì)P01全生命周期數(shù)據(jù)中提取退化特征和容量退化數(shù)據(jù)離線訓(xùn)練的結(jié)果如圖8所示。

    圖8 電池退化狀態(tài)模型離線訓(xùn)練結(jié)果

    所構(gòu)建的退化狀態(tài)模型都能夠隨著P01 老化過(guò)程中的容量退化量變化輸出準(zhǔn)確的容量退化量估算值,估算誤差不超過(guò)0.03 A·h,說(shuō)明離線訓(xùn)練得到的退化狀態(tài)模型能夠準(zhǔn)確描述容量退化量與退化特征之間的關(guān)系,估算結(jié)果表現(xiàn)出良好的跟蹤性。

    根據(jù)1.3 節(jié)中新陳代謝ELM 估算SOH 流程,用P01 離線訓(xùn)練的退化狀態(tài)模型對(duì)S17 和S30 在線SOH 估算。設(shè)置S17 和S30 已知的歷史數(shù)據(jù)點(diǎn)個(gè)數(shù)i=4,SOH 估算起始循環(huán)數(shù)k=5,估算結(jié)果如圖9所示。

    圖9 基于新陳代謝ELM的SOH估算結(jié)果

    基于P01數(shù)據(jù)離線構(gòu)建的新陳代謝ELM 估算的兩節(jié)電池(S17 和S30)SOH 與真實(shí)SOH 的退化軌跡都有較高的重合度。其中S17(圖9(a))的真實(shí)SOH退化軌跡除老化初始階段外相對(duì)平緩,因而估算誤差最大值出現(xiàn)在退化初期,后續(xù)SOH 估算誤差大致呈波動(dòng)收斂趨勢(shì);S30(圖9(b))的真實(shí)SOH 在退化中期階段出現(xiàn)一次明顯的波動(dòng),這是由于老化實(shí)驗(yàn)暫停休息期間電池容量出現(xiàn)暫時(shí)回升,導(dǎo)致該波動(dòng)點(diǎn)對(duì)應(yīng)的SOH 估算誤差達(dá)到極大值,之后隨著老化實(shí)驗(yàn)的繼續(xù)進(jìn)行,其真實(shí)SOH 變化軌跡逐漸趨于平緩,SOH 估算值亦能很好地跟蹤真實(shí)軌跡。S30 在退化中后期出現(xiàn)了明顯的加速退化現(xiàn)象,SOH 退化速率明顯加快,SOH估算誤差相較老化前期稍大,但其均未超過(guò)1%。表明不同工況下兩節(jié)電池的真實(shí)SOH 軌跡的退化過(guò)程中出現(xiàn)的任何波動(dòng),新陳代謝ELM都能準(zhǔn)確跟隨其變化趨勢(shì)。

    為進(jìn)一步驗(yàn)證本文所提出的電池SOH 估算框架的有效性,將未引入退化狀態(tài)模型和新陳代謝機(jī)制的傳統(tǒng)ELM 算法[19]用于SOH 對(duì)比估算:將P01 的退化特征序列{ΔR0}和{ΔR1}作為輸入數(shù)據(jù)、容量退化量序列{ΔC}作為目標(biāo)輸出訓(xùn)練傳統(tǒng)ELM,之后分別將S17 和S30 的退化特征序列輸入訓(xùn)練完成的傳統(tǒng)ELM 估算其SOH。傳統(tǒng)ELM 參數(shù)設(shè)置與新陳代謝ELM 完全一致(隱含層節(jié)點(diǎn)數(shù)L=30,激勵(lì)函數(shù)為Sigmoid),估算結(jié)果如圖10所示。

    圖10 基于傳統(tǒng)ELM的SOH估算結(jié)果

    由圖可見(jiàn),傳統(tǒng)ELM 的估算SOH 與真實(shí)SOH 的整體變化趨勢(shì)雖然在一定程度上比較相近,但估算值與真實(shí)值在電池退化中后期逐漸出現(xiàn)較大偏差。用最大絕對(duì)誤差、平均絕對(duì)誤差和均方根誤差來(lái)量化評(píng)估新陳代謝ELM 與傳統(tǒng)ELM 兩種算法對(duì)2 節(jié)電池的SOH估算精度,結(jié)果見(jiàn)表4。

    表4 SOH估算精度評(píng)估

    由表可見(jiàn),新陳代謝ELM 估算2 節(jié)電池SOH 的全局最大絕對(duì)誤差不超過(guò)2.18%,平均絕對(duì)誤差不超過(guò)0.6%且最大均方根誤差僅為0.009 7,說(shuō)明新陳代謝ELM 的SOH 估算結(jié)果不僅很接近真實(shí)值而且其分布的離散程度也很低。新陳代謝ELM 的各項(xiàng)誤差指標(biāo)都遠(yuǎn)小于傳統(tǒng)ELM 算法,表明電池退化狀態(tài)模型的構(gòu)建和新陳代謝機(jī)制的引入對(duì)不同類型電池的SOH 估算精度提升明顯,驗(yàn)證了所提出SOH估算方法相較于傳統(tǒng)ELM方法的優(yōu)越性。

    上述結(jié)果表明,基于所提取退化特征構(gòu)建的ELM電池退化狀態(tài)模型能夠準(zhǔn)確描述電池退化特性和規(guī)律,結(jié)合新陳代謝機(jī)制構(gòu)建的電池SOH 估算框架可在少數(shù)據(jù)樣本的情況下實(shí)現(xiàn)對(duì)不同工況和不同材料的電池SOH 準(zhǔn)確估算,具有良好的精度和泛化性。

    4 結(jié)論

    針對(duì)鋰離子電池健康狀態(tài)估算問(wèn)題,提出基于ELM 的電池退化狀態(tài)模型方法,并引入新陳代謝機(jī)制改進(jìn)估算流程,最終使用動(dòng)態(tài)工況下兩種材料的鋰離子電池實(shí)驗(yàn)數(shù)據(jù)進(jìn)行驗(yàn)證,得到如下結(jié)論。

    (1)在動(dòng)態(tài)工況下提取用于表征鋰離子電池內(nèi)部狀態(tài)退化規(guī)律的退化特征,并與電池容量退化進(jìn)行關(guān)聯(lián)性分析,驗(yàn)證所提取退化特征用于構(gòu)建電池退化狀態(tài)模型的合理性和有效性。

    (2)基于鋰離子電池全生命周期的退化特征和容量退化量,通過(guò)離線訓(xùn)練ELM 來(lái)建立用于描述電池退化規(guī)律的鋰離子電池退化狀態(tài)模型。在此基礎(chǔ)上引入新陳代謝機(jī)制,充分反映電池退化的最新趨勢(shì)和特征,構(gòu)建新陳代謝ELM估算電池SOH框架。

    (3)使用所提出的SOH 估算框架對(duì)不同正極材料、動(dòng)態(tài)放電工況和測(cè)試溫度的電池進(jìn)行測(cè)試,結(jié)果表明新陳代謝ELM 估算電池SOH 具有良好的精度和泛化性,在起始只用4 個(gè)歷史數(shù)據(jù)點(diǎn)的情況下估算得到的SOH 估算最大絕對(duì)誤差不超過(guò)2.18%,精度顯著優(yōu)于基于傳統(tǒng)ELM的SOH估算方法。

    猜你喜歡
    模型
    一半模型
    一種去中心化的域名服務(wù)本地化模型
    適用于BDS-3 PPP的隨機(jī)模型
    提煉模型 突破難點(diǎn)
    函數(shù)模型及應(yīng)用
    p150Glued在帕金森病模型中的表達(dá)及分布
    函數(shù)模型及應(yīng)用
    重要模型『一線三等角』
    重尾非線性自回歸模型自加權(quán)M-估計(jì)的漸近分布
    3D打印中的模型分割與打包
    亚洲一区二区三区不卡视频| 噜噜噜噜噜久久久久久91| 国产欧美日韩精品一区二区| 日韩欧美在线二视频| 久久午夜福利片| 国产精品人妻久久久久久| 色视频www国产| 夜夜夜夜夜久久久久| 麻豆成人午夜福利视频| 夜夜爽天天搞| 18禁黄网站禁片免费观看直播| 欧美绝顶高潮抽搐喷水| 久久久久性生活片| 18禁在线播放成人免费| 男女下面进入的视频免费午夜| 久久精品人妻少妇| 亚洲美女搞黄在线观看 | 女人被狂操c到高潮| 久久久久久大精品| 午夜福利欧美成人| 动漫黄色视频在线观看| 在线a可以看的网站| 91麻豆av在线| 99国产精品一区二区蜜桃av| 在线免费观看的www视频| 国内少妇人妻偷人精品xxx网站| 国产午夜精品论理片| 久久久久精品国产欧美久久久| 国产真实乱freesex| 欧美成人性av电影在线观看| 村上凉子中文字幕在线| 日本成人三级电影网站| 欧美+日韩+精品| 亚洲美女搞黄在线观看 | 精品欧美国产一区二区三| 免费观看人在逋| 五月伊人婷婷丁香| 免费一级毛片在线播放高清视频| 久久久久性生活片| 国产黄色小视频在线观看| 久久欧美精品欧美久久欧美| 亚洲成av人片在线播放无| 真人做人爱边吃奶动态| 最近最新中文字幕大全电影3| 色视频www国产| 免费在线观看成人毛片| 国产精华一区二区三区| 国产一区二区三区在线臀色熟女| 国产av不卡久久| 美女高潮喷水抽搐中文字幕| 欧美在线黄色| 看十八女毛片水多多多| 国产精品一区二区免费欧美| 老女人水多毛片| 真实男女啪啪啪动态图| 成人无遮挡网站| 俄罗斯特黄特色一大片| 最新中文字幕久久久久| 亚洲成人中文字幕在线播放| 变态另类成人亚洲欧美熟女| 亚洲第一电影网av| 精品无人区乱码1区二区| 中亚洲国语对白在线视频| 久久久久久国产a免费观看| 中文字幕人成人乱码亚洲影| 国产精品综合久久久久久久免费| 日韩精品中文字幕看吧| 丰满乱子伦码专区| 国产三级在线视频| 国产精品久久久久久久电影| 久久久久久久久久成人| 亚洲在线自拍视频| 国产精品亚洲av一区麻豆| 久久久国产成人免费| 99久久无色码亚洲精品果冻| a级毛片a级免费在线| 久久久久久久久大av| 日韩成人在线观看一区二区三区| 亚洲中文字幕一区二区三区有码在线看| 在线十欧美十亚洲十日本专区| 国产亚洲精品综合一区在线观看| 九色国产91popny在线| 一本精品99久久精品77| 精品久久久久久久久久久久久| 日韩 亚洲 欧美在线| 一边摸一边抽搐一进一小说| 我要搜黄色片| 国产高清有码在线观看视频| 丰满的人妻完整版| 日本黄大片高清| 欧美不卡视频在线免费观看| av福利片在线观看| 永久网站在线| av在线蜜桃| 亚洲经典国产精华液单 | 精品久久久久久,| 国产精品综合久久久久久久免费| 色哟哟哟哟哟哟| 色在线成人网| 亚洲成人精品中文字幕电影| 男女做爰动态图高潮gif福利片| 精品人妻视频免费看| 真人做人爱边吃奶动态| 内地一区二区视频在线| 精品人妻一区二区三区麻豆 | 中亚洲国语对白在线视频| 午夜亚洲福利在线播放| 欧美xxxx黑人xx丫x性爽| а√天堂www在线а√下载| 黄色女人牲交| 神马国产精品三级电影在线观看| 欧美三级亚洲精品| 黄色日韩在线| 欧美精品啪啪一区二区三区| 亚洲国产日韩欧美精品在线观看| 午夜精品一区二区三区免费看| 丰满人妻熟妇乱又伦精品不卡| 老司机午夜福利在线观看视频| 欧美丝袜亚洲另类 | 国产一区二区在线av高清观看| 有码 亚洲区| 中文字幕熟女人妻在线| av国产免费在线观看| 亚洲国产高清在线一区二区三| 99久久99久久久精品蜜桃| 天堂网av新在线| 国产av一区在线观看免费| 亚洲第一欧美日韩一区二区三区| 久久午夜亚洲精品久久| 亚洲国产精品sss在线观看| 白带黄色成豆腐渣| 内射极品少妇av片p| 在线观看66精品国产| av在线蜜桃| 午夜福利在线在线| 老熟妇仑乱视频hdxx| 日韩人妻高清精品专区| av在线蜜桃| 麻豆久久精品国产亚洲av| 日韩精品中文字幕看吧| 婷婷精品国产亚洲av| av在线蜜桃| 国产黄片美女视频| 午夜免费男女啪啪视频观看 | av国产免费在线观看| 欧美三级亚洲精品| 变态另类成人亚洲欧美熟女| 国产熟女xx| 91九色精品人成在线观看| 搡老熟女国产l中国老女人| 别揉我奶头 嗯啊视频| 美女高潮的动态| 欧美日韩福利视频一区二区| 禁无遮挡网站| 一本久久中文字幕| 国产精品一及| 亚洲乱码一区二区免费版| 在线观看av片永久免费下载| 久久久久久九九精品二区国产| 日韩欧美精品v在线| 欧美午夜高清在线| 97碰自拍视频| 亚洲成人久久性| 九九在线视频观看精品| 99视频精品全部免费 在线| 久久久久亚洲av毛片大全| 1000部很黄的大片| 欧美不卡视频在线免费观看| 亚洲国产精品合色在线| 国产午夜福利久久久久久| av在线观看视频网站免费| 88av欧美| 中文字幕熟女人妻在线| av在线天堂中文字幕| 搡老妇女老女人老熟妇| 性插视频无遮挡在线免费观看| 欧美激情在线99| 一边摸一边抽搐一进一小说| 久久久久性生活片| 黄色丝袜av网址大全| 看黄色毛片网站| 亚洲国产精品合色在线| 精品国产亚洲在线| 此物有八面人人有两片| 久久精品国产自在天天线| 亚洲av电影在线进入| 99国产精品一区二区蜜桃av| 丝袜美腿在线中文| 国产亚洲精品久久久com| 亚洲人成伊人成综合网2020| 淫妇啪啪啪对白视频| 国产欧美日韩精品一区二区| 人妻制服诱惑在线中文字幕| 在线a可以看的网站| 免费看a级黄色片| 一二三四社区在线视频社区8| 亚洲av第一区精品v没综合| 麻豆成人午夜福利视频| 国产综合懂色| 国产乱人伦免费视频| 精品久久久久久,| 我要看日韩黄色一级片| 69av精品久久久久久| 能在线免费观看的黄片| 成年女人毛片免费观看观看9| 亚洲av成人精品一区久久| 成年版毛片免费区| 又粗又爽又猛毛片免费看| 深夜精品福利| 嫁个100分男人电影在线观看| 国产高清三级在线| 一区二区三区激情视频| 他把我摸到了高潮在线观看| 亚洲欧美日韩卡通动漫| 久久草成人影院| 精品久久久久久久人妻蜜臀av| 精品一区二区三区人妻视频| 我的老师免费观看完整版| 日韩人妻高清精品专区| 色吧在线观看| 我要搜黄色片| 国产精品美女特级片免费视频播放器| 欧美国产日韩亚洲一区| 国产乱人伦免费视频| 亚洲综合色惰| 日韩av在线大香蕉| 日本黄色片子视频| 麻豆一二三区av精品| 国产高清视频在线播放一区| 欧美一区二区国产精品久久精品| 国内精品久久久久久久电影| 亚洲中文日韩欧美视频| 亚洲精品亚洲一区二区| 免费看美女性在线毛片视频| 亚洲精华国产精华精| 特大巨黑吊av在线直播| 国内精品久久久久久久电影| 国产久久久一区二区三区| 国产精品99久久久久久久久| 亚洲精品成人久久久久久| 毛片一级片免费看久久久久 | 久久久久国内视频| 老女人水多毛片| 亚洲国产日韩欧美精品在线观看| 黄色日韩在线| 国产精品一区二区性色av| 老熟妇仑乱视频hdxx| 日韩欧美国产一区二区入口| 婷婷六月久久综合丁香| 97超视频在线观看视频| 琪琪午夜伦伦电影理论片6080| 欧美日韩乱码在线| 欧美成人一区二区免费高清观看| 亚洲成人久久爱视频| 女人被狂操c到高潮| 国产aⅴ精品一区二区三区波| 国产在线男女| 国产乱人伦免费视频| 又粗又爽又猛毛片免费看| 国产成年人精品一区二区| 一级毛片久久久久久久久女| 国模一区二区三区四区视频| 亚洲人与动物交配视频| 亚洲国产精品久久男人天堂| 18+在线观看网站| 小说图片视频综合网站| 自拍偷自拍亚洲精品老妇| 色综合亚洲欧美另类图片| a级毛片a级免费在线| 欧美+日韩+精品| 少妇高潮的动态图| 日本黄大片高清| 国产一区二区在线观看日韩| 国产成人福利小说| 成年人黄色毛片网站| 久久精品国产99精品国产亚洲性色| 99热这里只有是精品50| 特级一级黄色大片| 国产高清视频在线播放一区| av欧美777| 丰满乱子伦码专区| 亚洲五月天丁香| 免费观看的影片在线观看| 男人狂女人下面高潮的视频| 人妻丰满熟妇av一区二区三区| 亚洲经典国产精华液单 | 日韩欧美精品免费久久 | 又粗又爽又猛毛片免费看| 久久99热这里只有精品18| 国产精品精品国产色婷婷| 日韩亚洲欧美综合| 别揉我奶头 嗯啊视频| 午夜免费男女啪啪视频观看 | 噜噜噜噜噜久久久久久91| 精品熟女少妇八av免费久了| 国产精品久久电影中文字幕| 亚洲无线观看免费| 中文字幕av成人在线电影| 亚州av有码| 国模一区二区三区四区视频| 亚洲第一欧美日韩一区二区三区| 在线十欧美十亚洲十日本专区| 亚洲真实伦在线观看| 88av欧美| 欧美成人一区二区免费高清观看| 欧美黄色淫秽网站| 99久久精品一区二区三区| 日韩中文字幕欧美一区二区| 在线观看一区二区三区| 极品教师在线免费播放| 亚洲欧美日韩高清专用| 国产久久久一区二区三区| 一本综合久久免费| 好男人电影高清在线观看| 亚洲av.av天堂| 嫩草影院入口| 亚洲avbb在线观看| 熟妇人妻久久中文字幕3abv| 午夜福利高清视频| 五月玫瑰六月丁香| 亚洲五月婷婷丁香| 中文在线观看免费www的网站| 麻豆成人av在线观看| 久久久久久久久大av| 亚洲精华国产精华精| 婷婷丁香在线五月| 老司机深夜福利视频在线观看| 色在线成人网| 少妇裸体淫交视频免费看高清| 最新在线观看一区二区三区| 免费一级毛片在线播放高清视频| 日韩欧美 国产精品| АⅤ资源中文在线天堂| 一边摸一边抽搐一进一小说| 国产精品永久免费网站| 在线观看免费视频日本深夜| 久久婷婷人人爽人人干人人爱| 亚洲七黄色美女视频| 国产亚洲欧美98| 看黄色毛片网站| 亚洲精品亚洲一区二区| 深爱激情五月婷婷| 中亚洲国语对白在线视频| 久久久久久久亚洲中文字幕 | 久久久久久久精品吃奶| 校园春色视频在线观看| 两性午夜刺激爽爽歪歪视频在线观看| 99国产极品粉嫩在线观看| 波多野结衣高清作品| 日韩高清综合在线| 亚洲无线观看免费| 久久精品夜夜夜夜夜久久蜜豆| 久久久久国内视频| 国产av麻豆久久久久久久| 欧美三级亚洲精品| 少妇的逼好多水| 欧美日韩福利视频一区二区| 精品久久久久久久人妻蜜臀av| 好男人电影高清在线观看| 久久久久久国产a免费观看| 午夜精品久久久久久毛片777| 免费在线观看日本一区| 久久精品91蜜桃| 波多野结衣高清作品| 久久久久九九精品影院| 一本精品99久久精品77| 一区二区三区激情视频| 亚洲av免费在线观看| 成年人黄色毛片网站| 免费在线观看日本一区| 一进一出抽搐gif免费好疼| 亚洲国产欧美人成| 中文字幕精品亚洲无线码一区| 欧美日本视频| 午夜福利18| 日日干狠狠操夜夜爽| 一级av片app| 一个人看视频在线观看www免费| 少妇被粗大猛烈的视频| 色精品久久人妻99蜜桃| 天堂影院成人在线观看| 精品不卡国产一区二区三区| 黄片小视频在线播放| a在线观看视频网站| 成人美女网站在线观看视频| 国产乱人视频| 嫩草影院精品99| 欧美黄色淫秽网站| 一个人看视频在线观看www免费| 搡老岳熟女国产| 亚洲av第一区精品v没综合| 欧美另类亚洲清纯唯美| 婷婷亚洲欧美| 久久久久久久久大av| 国产成人影院久久av| 日日夜夜操网爽| 中文字幕av成人在线电影| 国内毛片毛片毛片毛片毛片| 蜜桃亚洲精品一区二区三区| 亚洲欧美日韩东京热| 精品国产三级普通话版| 一区二区三区免费毛片| 国产大屁股一区二区在线视频| 九九在线视频观看精品| 床上黄色一级片| 很黄的视频免费| 午夜免费男女啪啪视频观看 | 最后的刺客免费高清国语| 成人美女网站在线观看视频| 老熟妇仑乱视频hdxx| 亚洲性夜色夜夜综合| 久久久国产成人精品二区| 日韩精品中文字幕看吧| 婷婷色综合大香蕉| 成人无遮挡网站| 免费一级毛片在线播放高清视频| 欧美又色又爽又黄视频| 男女之事视频高清在线观看| 国产又黄又爽又无遮挡在线| 国产av在哪里看| 69人妻影院| av黄色大香蕉| 国产爱豆传媒在线观看| 超碰av人人做人人爽久久| 婷婷丁香在线五月| 亚洲一区二区三区色噜噜| 国产一区二区三区在线臀色熟女| 婷婷丁香在线五月| 欧美一区二区精品小视频在线| 嫩草影院精品99| 丰满人妻熟妇乱又伦精品不卡| 老熟妇仑乱视频hdxx| 两个人视频免费观看高清| 国产精品电影一区二区三区| 一进一出抽搐gif免费好疼| 偷拍熟女少妇极品色| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 亚洲精品亚洲一区二区| 亚洲精华国产精华精| 日日夜夜操网爽| 欧美在线黄色| 欧美乱妇无乱码| 国产亚洲精品综合一区在线观看| 少妇丰满av| 午夜福利视频1000在线观看| 欧美激情国产日韩精品一区| 国产精品一及| 精品无人区乱码1区二区| 白带黄色成豆腐渣| 亚洲熟妇熟女久久| 欧美3d第一页| 在线天堂最新版资源| 国产一区二区激情短视频| 亚洲第一欧美日韩一区二区三区| 国产老妇女一区| 51午夜福利影视在线观看| 亚洲在线自拍视频| 97人妻精品一区二区三区麻豆| 久久久久国内视频| 91九色精品人成在线观看| 久久婷婷人人爽人人干人人爱| 别揉我奶头~嗯~啊~动态视频| 国产成+人综合+亚洲专区| 久久久成人免费电影| 国产麻豆成人av免费视频| 波野结衣二区三区在线| 欧美性猛交黑人性爽| 色av中文字幕| 人妻夜夜爽99麻豆av| 乱人视频在线观看| 国产国拍精品亚洲av在线观看| x7x7x7水蜜桃| 国产精品不卡视频一区二区 | 最近最新免费中文字幕在线| 日韩国内少妇激情av| 亚洲熟妇熟女久久| 老司机午夜福利在线观看视频| 中文字幕av在线有码专区| 国产午夜精品久久久久久一区二区三区 | 中文字幕久久专区| 久久精品久久久久久噜噜老黄 | 12—13女人毛片做爰片一| 久久6这里有精品| 亚洲中文日韩欧美视频| 香蕉av资源在线| 欧美中文日本在线观看视频| 成人特级黄色片久久久久久久| 91在线精品国自产拍蜜月| 啦啦啦韩国在线观看视频| 精品久久久久久久久亚洲 | 亚洲五月天丁香| 国产亚洲精品综合一区在线观看| 亚洲美女视频黄频| 色视频www国产| 国内精品美女久久久久久| а√天堂www在线а√下载| 看片在线看免费视频| 搡女人真爽免费视频火全软件 | 又爽又黄无遮挡网站| 欧美日韩黄片免| 啪啪无遮挡十八禁网站| 99热精品在线国产| 久久99热这里只有精品18| 亚洲一区高清亚洲精品| 亚洲精品在线美女| 69av精品久久久久久| 看黄色毛片网站| 国产精品,欧美在线| 婷婷丁香在线五月| 久久精品人妻少妇| 欧美xxxx黑人xx丫x性爽| 内射极品少妇av片p| 国产亚洲欧美在线一区二区| 国产高清有码在线观看视频| 看黄色毛片网站| 十八禁网站免费在线| 欧美高清性xxxxhd video| 无遮挡黄片免费观看| 精品无人区乱码1区二区| 麻豆久久精品国产亚洲av| 丝袜美腿在线中文| 国产精品99久久久久久久久| 国产一区二区激情短视频| 天堂影院成人在线观看| 丁香欧美五月| 欧美3d第一页| 国产精品久久久久久亚洲av鲁大| 给我免费播放毛片高清在线观看| 亚洲欧美日韩东京热| 舔av片在线| 日日夜夜操网爽| 国产91精品成人一区二区三区| 亚洲成a人片在线一区二区| 国产高潮美女av| 十八禁人妻一区二区| 国产精品国产高清国产av| 90打野战视频偷拍视频| 日日夜夜操网爽| 亚洲精品色激情综合| 午夜久久久久精精品| 午夜a级毛片| 中出人妻视频一区二区| 亚洲精品亚洲一区二区| 一本一本综合久久| 亚洲中文字幕日韩| 欧美xxxx黑人xx丫x性爽| 日本熟妇午夜| 欧美一区二区国产精品久久精品| 夜夜躁狠狠躁天天躁| 久久亚洲精品不卡| 日本黄色片子视频| 91午夜精品亚洲一区二区三区 | 高清毛片免费观看视频网站| 黄色配什么色好看| 少妇裸体淫交视频免费看高清| 亚洲av中文字字幕乱码综合| 国产成年人精品一区二区| 日韩欧美国产在线观看| 亚洲成人中文字幕在线播放| 精品人妻熟女av久视频| 色哟哟哟哟哟哟| 俺也久久电影网| 国产老妇女一区| 欧美成人一区二区免费高清观看| 人妻制服诱惑在线中文字幕| 搡老妇女老女人老熟妇| 搞女人的毛片| 超碰av人人做人人爽久久| 中文资源天堂在线| 18禁在线播放成人免费| 国产精品人妻久久久久久| 国产三级黄色录像| 日日夜夜操网爽| 变态另类丝袜制服| 国内精品久久久久精免费| 亚洲国产欧洲综合997久久,| 中文字幕久久专区| 波多野结衣高清作品| 欧美潮喷喷水| 国产亚洲欧美在线一区二区| 人人妻人人看人人澡| 欧美一区二区国产精品久久精品| 免费一级毛片在线播放高清视频| 99在线人妻在线中文字幕| 精品久久久久久久人妻蜜臀av| 日本精品一区二区三区蜜桃| 99热这里只有是精品50| 香蕉av资源在线| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 亚洲自偷自拍三级| 久久久精品大字幕| 欧美精品国产亚洲| av专区在线播放| 色综合站精品国产| 麻豆一二三区av精品| 亚洲狠狠婷婷综合久久图片| 小蜜桃在线观看免费完整版高清| 色综合站精品国产| 9191精品国产免费久久| 天美传媒精品一区二区| 老熟妇乱子伦视频在线观看| 国产成人影院久久av| 日本黄大片高清| 简卡轻食公司| 午夜免费成人在线视频| 国产精品伦人一区二区| 国产一区二区在线观看日韩| www.熟女人妻精品国产| 欧美黑人欧美精品刺激| 美女被艹到高潮喷水动态| 欧美日韩中文字幕国产精品一区二区三区| netflix在线观看网站| 偷拍熟女少妇极品色| 伦理电影大哥的女人| 国产精华一区二区三区|