高 原,黃榮杰,秦東晨,王婷婷

(1.鄭州大學(xué)機(jī)械與動(dòng)力工程學(xué)院,河南 鄭州 450001;2.鄭州輕工業(yè)大學(xué)機(jī)電工程學(xué)院,河南 鄭州 450001)

電池剩余電量、壽命、溫度等狀態(tài)的準(zhǔn)確估算是保證電動(dòng)汽車(chē)安全、高效運(yùn)行的前提[1],并且都依賴(lài)于準(zhǔn)確的電池模型。等效電路模型能很好地反映鋰離子電池的外部特性,被人們廣泛采用。等效電路模型的內(nèi)部參數(shù)通?？梢哉J(rèn)為是關(guān)于溫度[2]和荷電狀態(tài)(SOC)[3]的函數(shù),但是沒(méi)有考慮電流引起的參數(shù)變化。電動(dòng)汽車(chē)的行駛工況復(fù)雜多變,在高速行駛、爬坡等階段,鋰離子電池的放電電流通常會(huì)達(dá)到3C,而在一些特種電動(dòng)車(chē)輛中,峰值放電電流甚至可達(dá)到10C。一方面,電池容量、電阻等參數(shù)在高倍率放電時(shí)會(huì)發(fā)生改變;另一方面,電流的增大會(huì)導(dǎo)致電池溫度明顯上升,影響電池模型的精度。

基于以上分析,本文作者提出一種綜合考慮電流和溫度影響的鋰離子電池電熱耦合模型,以追蹤高倍率電流工況下電池參數(shù)的變化和溫度的改變;通過(guò)混合功率脈沖特性(HPPC)實(shí)驗(yàn),獲取不同倍率電流工況、溫度以及SOC區(qū)間的阻抗參數(shù);最后,對(duì)電熱耦合模型在高倍率電流工況和不同溫度下的動(dòng)態(tài)工況進(jìn)行驗(yàn)證。

1 電池模型

1.1 考慮電流和溫度影響的鋰離子電池電模型

對(duì)傳統(tǒng)等效電路模型進(jìn)行改進(jìn),將電池參數(shù)設(shè)置為電流、溫度和SOC相關(guān)的變量,以模擬鋰離子電池在不同負(fù)載電流、電池溫度和放電區(qū)間下參數(shù)的變化。改進(jìn)后的模型如圖1所示。

圖1 考慮電流和溫度影響的等效電路模型Fig.1 Equivalent circuitmodel considering the effect of current and temperature

圖1中:R0為歐姆內(nèi)阻;R1、C1分別為電化學(xué)極化電阻和電容;R2、C2分別為濃度差極化電阻和電容;I為負(fù)載電流;θ為溫度;Uocv為開(kāi)路電壓;Soc為電池的SOC;U0為歐姆內(nèi)阻的端電壓;U1、U2分別為兩個(gè)RC網(wǎng)絡(luò)的端電壓;UL為電池的端電壓。R0、R1、C1、R2和C2均為關(guān)于電流、溫度和SOC的變量。

根據(jù)基爾霍夫定律,改進(jìn)后等效電路的數(shù)學(xué)模型可以表示為:

1.2 鋰離子電池?zé)崮Ｐ?/h3>

鋰離子電池?zé)崃縼?lái)自于自身內(nèi)部產(chǎn)生的化學(xué)熱和與外界的交換熱。假設(shè)鋰離子電池外表面各處的溫度是一致的,電池與外界環(huán)境的熱交換也是穩(wěn)定的,結(jié)合熱力學(xué)第一定律,鋰離子電池?zé)崮Ｐ涂梢员硎緸?

式(2)中:Qgen為電池產(chǎn)生的化學(xué)熱;Qexc為電池的交換熱;m為電池質(zhì)量;Cs為電池的比熱容;h為電池的熱交換系數(shù);S為電池的表面積;θamb為環(huán)境溫度;θcell為電池的表面溫度;t為時(shí)刻。

1.3 鋰離子電池電熱耦合模型

電模型輸出的電池參數(shù)是熱模型的輸入量,同時(shí),熱模型輸出的電池溫度也是電模型的輸入量,兩者相互耦合,最終輸出電池的端電壓。

根據(jù)式(1)、(2)建立的鋰離子電池電熱耦合框架如圖2所示。

圖2 電熱耦合框架圖Fig.2 Frame diagram of the electro-thermal coupling

該模型需要識(shí)別的參數(shù)有:阻抗參數(shù)R0、R1、C1、R2和C2,以及熱交換系數(shù)h等。

2 高倍率電流工況下電池參數(shù)識(shí)別

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

實(shí)驗(yàn)對(duì)象為電動(dòng)汽車(chē)用20 Ah磷酸鐵鋰正極鋰離子電池(廣東產(chǎn)),主要參數(shù)為:標(biāo)稱(chēng)電壓3.2 V,工作電壓2.50～3.65 V,比熱容1.06 J/(g·K)。

電池測(cè)試平臺(tái)所用儀器包括:控制充放電和監(jiān)控電池狀態(tài)的CT-6001A電池測(cè)試系統(tǒng)(武漢產(chǎn)),控制電池環(huán)境溫度的BLC-300溫控箱(上海產(chǎn)),監(jiān)控電池表面溫度的UTi-380H熱成像儀(東莞產(chǎn)),以及用于數(shù)據(jù)處理的主機(jī)。

2.2 阻抗參數(shù)識(shí)別和分析

分別在5℃、15℃、25℃、35℃和45℃下,在100%～0 SOC范圍內(nèi),每隔10% SOC,分別以0.5C、1.0C、2.0C和3.0C的電流在2.50～3.65 V進(jìn)行HPPC測(cè)試[4],識(shí)別電池阻抗參數(shù)。

鋰離子電池的歐姆內(nèi)阻R0可以根據(jù)放電脈沖開(kāi)始和結(jié)束時(shí)的電壓下降量來(lái)計(jì)算,極化參數(shù)R1、C1、R2和C2可以通過(guò)擬合HPPC實(shí)驗(yàn)中的電壓曲線(xiàn)獲得,具體的計(jì)算步驟參見(jiàn)文獻(xiàn)[5]。

最終識(shí)別出的阻抗參數(shù)與電流倍率、溫度和SOC的關(guān)系如圖3和圖4所示。

圖3 不同電流倍率下鋰離子電池阻抗參數(shù)與溫度的關(guān)系Fig.3 Relation between the impedance parameters and temperatures of battery at different current rates

圖4 不同溫度下鋰離子電池阻抗參數(shù)與SOC的關(guān)系Fig.4 Relation between the impedance parameters and state of charge(SOC)of battery at different temperatures

從圖3、4可知,R0、R1和R2受電流影響的趨勢(shì)相同,均隨著放電電流的增大而減小。巴特勒-福爾默方程表明,電池的電荷轉(zhuǎn)移內(nèi)阻隨著放電電流的增大而減小,同時(shí),電流較小時(shí),Li+的能量較低,躍過(guò)界面壁壘能力較差,因此,極化電阻更大。歐姆內(nèi)阻R0及極化內(nèi)阻R1、R2都隨著溫度的下降而升高,原因是低溫環(huán)境下,活性物質(zhì)的擴(kuò)散能力、電解液的導(dǎo)電率都會(huì)急劇下降,造成電阻下降。

2.3 熱特性參數(shù)識(shí)別

熱模型中的熱交換系數(shù)h可通過(guò)電池高倍率電流放電,溫度急劇升高后,靜置階段的冷卻溫度曲線(xiàn)獲得[6]。實(shí)驗(yàn)記錄的電池中心溫度為38.5℃,最高溫度點(diǎn)為39.1℃,兩者僅相差0.6℃。為方便計(jì)算,假設(shè)電池表面各處溫度一致,在無(wú)負(fù)載電流狀態(tài)下,電池的密度、比熱容和熱交換系數(shù)h及所處的環(huán)境溫度是恒定的。

靜置階段電池的發(fā)熱量Qgen(t)=0,求解式(2),可得到表面溫度隨時(shí)間變化的方程。

式(3)中:a為擬合系數(shù);exp是以自然常數(shù)e為底的指數(shù)函數(shù)。

根據(jù)在25℃環(huán)境下進(jìn)行的電池冷卻實(shí)驗(yàn)及式(3),擬合得到的熱交換系數(shù)h=16.04W/(m2·℃ )。

3 模型驗(yàn)證與分析

根據(jù)電熱耦合模型和識(shí)別的參數(shù),在MATLAB/Simulink環(huán)境下搭建鋰離子電池仿真模型。

為驗(yàn)證電熱耦合模型的準(zhǔn)確性,分別對(duì)模型在不同倍率恒流放電工況和不同溫度下進(jìn)行4個(gè)連續(xù)的新歐洲行駛循環(huán)(NEDC)工況(電流曲線(xiàn)見(jiàn)圖5)仿真,并將得到的端電壓、電池溫度與實(shí)驗(yàn)數(shù)據(jù)進(jìn)行對(duì)比。

圖5 NEDC工況的電流曲線(xiàn)Fig.5 Current curve of new European driving cycle(NEDC)condition

搭建不考慮電流影響的電池模型作為對(duì)比。該電池模型的參數(shù)是在與電熱耦合模型相同的條件下,以0.5C倍率實(shí)驗(yàn)得到的。

采用如式(4)所示的均方根誤差(RMSE,ERMSE)和平均絕對(duì)誤差(MAE,EMAE)來(lái)評(píng)價(jià)模型的準(zhǔn)確性。

式(4)中:y t為仿真值;u t為實(shí)驗(yàn)數(shù)據(jù);k為樣點(diǎn)個(gè)數(shù)。

3.1 高倍率恒流放電工況驗(yàn)證

鋰離子電池在高倍率恒流放電工況下仿真和實(shí)驗(yàn)所得溫度和端電壓數(shù)據(jù)如圖6所示。

圖6 15℃下高倍率恒流放電工況驗(yàn)證結(jié)果Fig.6 Verification results under high rate constant current discharge conditions at 15℃

從圖6可知,電熱耦合模型在1.0～4.0C電流倍率下的仿真結(jié)果均與實(shí)驗(yàn)數(shù)據(jù)吻合良好。在1.0C電流倍率下,是否考慮電流影響的端電壓相差不多,但當(dāng)電流加大時(shí),兩者發(fā)生明顯偏離。電池內(nèi)阻隨著放電倍率的增大而減小,因此在2.0～4.0C電流倍率下,不考慮電流影響的對(duì)比模型(參數(shù)識(shí)別電流為0.5C)所得的仿真電壓,低于實(shí)驗(yàn)數(shù)據(jù)和電熱耦合模型仿真電壓。以3.0C、4.0C倍率放電時(shí),實(shí)驗(yàn)和仿真的電池端電壓都出現(xiàn)了反彈,而不是隨著放電時(shí)間延長(zhǎng)而持續(xù)減小。這是由于電池在高倍率放電工況下,溫度明顯升高,如圖6(a)所示,在3.0C和4.0C倍率下,電池的溫升分別為17.8℃和21.2℃。電池溫度升高,內(nèi)阻顯著降低,因此端電壓會(huì)升高反彈。

3.2 NEDC工況驗(yàn)證

在15℃下,NEDC工況仿真得到的電池端電壓、實(shí)驗(yàn)數(shù)據(jù)及誤差如圖7所示。

圖7 15℃下NEDC工況驗(yàn)證結(jié)果Fig.7 Verification results under NEDC condition at 15℃

從圖7可知,電熱耦合模型仿真結(jié)果與實(shí)驗(yàn)數(shù)據(jù)高度吻合,模型表現(xiàn)整體穩(wěn)定,沒(méi)有隨溫度和電流的改變而出現(xiàn)較大的誤差。仿真結(jié)果的最大誤差為17.1 mV,相比不考慮電流影響的電池模型減少了82.2%。兩種模型與實(shí)驗(yàn)數(shù)據(jù)的走勢(shì)基本相同,所得端電壓都隨著時(shí)間和放電電流的增大而減小,但不同區(qū)間兩種模型的誤差略有不同。在NEDC工況的低倍率電流區(qū)間,電池負(fù)載電流較小,電流的小范圍波動(dòng)不會(huì)引起電池參數(shù)的明顯變化,因此兩種模型的仿真結(jié)果均與實(shí)驗(yàn)數(shù)據(jù)接近。在NEDC工況的高倍率電流區(qū)間,電流的急劇增大引起電池容量的減小和電阻的降低,此時(shí),電流倍率成為影響電池參數(shù)的重要因素,若忽略電流的影響,會(huì)增大仿真誤差。

兩種模型在5～45℃內(nèi)NEDC工況下的RMSE和MAE見(jiàn)表1。

表1 不同模型的端電壓誤差Table 1 Terminal voltage errors of differentmodels

從表1可知,在所有的測(cè)試溫度范圍之內(nèi),電熱耦合模型的精度均有明顯改善。相比不考慮電流影響的電池模型,在5℃、15℃、25℃、35℃和45℃下,RMSE分別減少了45.5%、46.9%、41.6%、50.6%和 32.1%;MAE分別減少了35.4%、41.5%、35.1%、55.0%和30.2%。

綜上所述,考慮電流影響的電熱耦合模型,可以很好地模擬出寬溫度和電流倍率范圍內(nèi)鋰離子電池的電壓響應(yīng)。改進(jìn)后的鋰離子電池模型適用于電動(dòng)汽車(chē),尤其是高倍率電流運(yùn)行場(chǎng)景。

4 結(jié)論

為提升高倍率工況下鋰離子電池模型的精度,本文作者構(gòu)建了一種考慮電流和溫度影響的鋰離子電池電熱耦合模型,主要結(jié)論如下:

高倍率電流會(huì)引起鋰離子電池內(nèi)部參數(shù)的改變,造成電池溫度的急劇上升;同時(shí),溫度也會(huì)影響鋰離子電池的內(nèi)部參數(shù)。在搭建電池模型時(shí),應(yīng)考慮兩者的電熱耦合關(guān)系。

驗(yàn)證結(jié)果表明,與現(xiàn)有模型相比,考慮電流影響的電熱耦合模型可以大幅減小高倍率放電工況下的誤差,在NEDC工況5℃、15℃、25℃、35℃和45℃下的MAE分別減小了35.4%、41.5%、35.1%、55.0%和30.2%,具有更好的模型精度和穩(wěn)定性。