張立軍， 程洪正

(1.同濟(jì)大學(xué) 汽車學(xué)院,上海 201804; 2.同濟(jì)大學(xué) 智能型新能源汽車協(xié)同創(chuàng)新中心,上海 201804)

?

基于相似理論的鋰電池三維電化學(xué)跨尺度建模

張立軍1,2， 程洪正1,2

(1.同濟(jì)大學(xué) 汽車學(xué)院,上海 201804; 2.同濟(jì)大學(xué) 智能型新能源汽車協(xié)同創(chuàng)新中心,上海 201804)

摘要：從拋物線型偏微分方程基本形式出發(fā)，利用量綱為一分析法推導(dǎo)了拋物線型偏微分方程的相似準(zhǔn)則，并推廣得到電化學(xué)過程的相似準(zhǔn)則，從而確定原型與模型之間的電化學(xué)過程參數(shù)相似關(guān)系.利用Comsol Multiphysics軟件進(jìn)行了鋰離子電池單元電化學(xué)過程的有限元建模與計(jì)算，驗(yàn)證了所建立相似準(zhǔn)則和相似關(guān)系的正確性.模型驗(yàn)證結(jié)果表明電化學(xué)有限元建模方法能實(shí)現(xiàn)鋰離子電池三維電化學(xué)過程的一致模擬，同時(shí)，高倍率充放電工況下電池內(nèi)部不均勻性加劇，導(dǎo)致傳統(tǒng)的通用準(zhǔn)二維模型預(yù)測結(jié)果具有較大誤差.因此，在電池單體設(shè)計(jì)及系統(tǒng)集成與管理過程中，有必要對基于準(zhǔn)二維模型的狀態(tài)估計(jì)算法進(jìn)行修正，以避免電池局部過充過放而導(dǎo)致安全隱患.

關(guān)鍵詞：鋰離子電池; 電化學(xué); 相似理論

動(dòng)力電池是電動(dòng)汽車的最核心部件之一，其性能影響到電動(dòng)汽車的綜合性能. 近年來，鋰離子電池因其優(yōu)異的綜合性能倍受關(guān)注[1]. 鋰離子電池性能設(shè)計(jì)包括結(jié)構(gòu)、電極及隔膜材料、電流密度、容量、電壓等的合理設(shè)計(jì)，對于電池各項(xiàng)外特性如功率與能量特性、效率、安全性、可靠性和壽命等都具有重要影響[2-4].

通常利用鋰離子電池電化學(xué)模型進(jìn)行電池電極、隔膜材料選擇、結(jié)構(gòu)設(shè)計(jì)，建立電池內(nèi)部電化學(xué)動(dòng)力學(xué)與外特性(工作電壓、容量、功率特性)間的關(guān)系，在電池設(shè)計(jì)與優(yōu)化方面有著廣泛應(yīng)用[5-6].

國內(nèi)外研究中所提出的鋰離子電池電化學(xué)模型主要分為經(jīng)驗(yàn)?zāi)Ｐ?、原理性模? 經(jīng)驗(yàn)?zāi)Ｐ椭饕械刃щ娐?、阻抗譜模型等，其參數(shù)依賴于試驗(yàn)數(shù)據(jù)擬合，參數(shù)沒有直接的電化學(xué)意義，只能反映電池外特性且預(yù)測精度有限[7]. 原理性模型分為單粒子模型(single particle model)[8-9]和通用準(zhǔn)二維模型(Pseudo-2D)[10-12]. 單粒子模型由Haran等人[8]提出，單粒子模型忽略液相鋰離子濃度和電勢分布，僅在小于1C的低放電倍率時(shí)模型預(yù)測精度尚可[13]. 準(zhǔn)二維模型由Newman， Fuller, Doyle等人[10-11]基于多孔電極理論和濃溶液理論建立. 隨后建立起鋰離子電池準(zhǔn)三維[14]電化學(xué)模型，或考慮電化學(xué)—熱耦合模型[15-17]. 國內(nèi)外學(xué)者所建立的模型一般將電池單元簡化為一維結(jié)構(gòu)[10-12]，只考慮液相濃度、SOC、電流密度、生熱速率等沿厚度方向分布；大多未充分考慮液相濃度、局部SOC(state of charge)、電流密度、電勢、生熱速率等沿極片面內(nèi)分布的不均勻性. 鋰離子電池內(nèi)部的電流密度、局部濃度、SOC、電勢等的極片面內(nèi)方向分布并不均勻.在電動(dòng)車用工況條件下，不均勻性加劇，導(dǎo)致采用以上簡化模型與實(shí)際分布將具有較大偏差[15-17]. 因此，有必要建立嚴(yán)格的三維電化學(xué)模型，考慮極片面內(nèi)電化學(xué)狀態(tài)的不均勻性，提高電池內(nèi)部電化學(xué)狀態(tài)的預(yù)測精度.

電動(dòng)車用鋰離子電池單體一般包含上百個(gè)圖1所示電池單元，每個(gè)電池單元又包括正極集流體、正極、隔膜、負(fù)極、負(fù)極集流體五層堆疊結(jié)構(gòu).由圖1可知，單元各層結(jié)構(gòu)的長度和寬度方向尺寸為分米級，而厚度方向尺寸為微米級，具有顯著的跨尺度特征. 若直接按照電池單元的實(shí)際疊層結(jié)構(gòu)尺寸進(jìn)行電池內(nèi)部電化學(xué)過程有限元建模與計(jì)算，將導(dǎo)致模型單元數(shù)量多、計(jì)算量大、難以實(shí)現(xiàn)的嚴(yán)重困難. 因此，如何解決跨尺度建模問題成為利用有限元方法進(jìn)行電池內(nèi)部電化學(xué)過程的預(yù)測的關(guān)鍵. 本文提出了應(yīng)用于動(dòng)力鋰離子電池三維電化學(xué)過程跨尺度模擬方法并進(jìn)行了驗(yàn)證分析.

圖1鋰離子單電池結(jié)構(gòu)示意圖與三維尺寸(單位:mm)

Fig.1Schematic of lithium-ion battery cell(Unit: mm)

相似理論方法是人們在探索自然規(guī)律的過程中形成的一種?；椒?，對流體力學(xué)、傳熱學(xué)等學(xué)科的發(fā)展曾發(fā)揮了不可估量的作用[18]. 在前期研究中，曾將相似理論應(yīng)用于鋰離子電池的溫度場[19]、電場[20]及熱應(yīng)力場[21]的建模分析，有效地解決鋰離子電池跨尺度有限元建模問題.

本文從拋物線型偏微分方程基本形式出發(fā)，利用量綱為一分析法推導(dǎo)了拋物線型偏微分方程的相似準(zhǔn)則，并推廣得到電化學(xué)過程控制方程的相似準(zhǔn)則. 利用Comsol Multiphysics有限元軟件進(jìn)行了相似準(zhǔn)則和相似關(guān)系正確性的驗(yàn)證，模型驗(yàn)證結(jié)果表明基于相似原理的電化學(xué)有限元建模方法能實(shí)現(xiàn)鋰離子電池三維電化學(xué)過程的一致模擬，從而為鋰離子電池的電性能設(shè)計(jì)、電化學(xué)狀態(tài)估計(jì)與管理奠定良好的基礎(chǔ).

1鋰離子電池電化學(xué)過程基本方程

1.1研究對象

圖2所示為磷酸鐵鋰(LiFePO4)動(dòng)力電池單元內(nèi)部過程示意圖.正極活性材料為磷酸鐵鋰(LiFePO4)，負(fù)極活性材料為石墨(C6)，正、負(fù)極集流體材料分別為鋁和銅，電解液為1 mol·L-1的LiPF6溶液(1∶2 EC/DMC 溶劑).

圖2　鋰離子電池內(nèi)部過程示意圖

電化學(xué)反應(yīng)發(fā)生在多孔電極固、液相界面上，反應(yīng)方程式如式(1)和式(2)所示.

正極反應(yīng)

(1)

負(fù)極反應(yīng)

(2)

式中,x為充放電過程中脫出或嵌入正數(shù)Li的質(zhì)量分?jǐn)?shù).

1.2電化學(xué)過程控制方程

基于多孔電極理論、濃溶液理論、歐姆定律及電極反應(yīng)動(dòng)力學(xué)、物質(zhì)與電荷能量守恒理論建立鋰離子電池電化學(xué)過程模型.建模過程中對鋰離子電池電化學(xué)過程進(jìn)行了合理地簡化假設(shè)：①假設(shè)電極活性材料由均勻球形顆粒構(gòu)成；②電極層中只有固、液兩相進(jìn)行反應(yīng)，無副反應(yīng)發(fā)生; ③電化學(xué)反應(yīng)符合Bulter-Volmer方程；④通過擴(kuò)散和遷移在電解液里傳輸,忽略對流作用；⑤充放電過程體積膨脹忽略不計(jì)，孔隙率為恒定數(shù)值；⑥活性物質(zhì)固體相擴(kuò)散系數(shù)、離子遷移數(shù)與濃度、溫度無關(guān)，為恒定數(shù)值.

1.2.1電荷守恒方程

在多孔活性電極保持電中性條件下，由電流的連續(xù)性可得以下方程

(3)

由歐姆定律，并考慮液相遷移過程有

(4)

(5)

固相電荷守恒方程如下:

(6)

液相電荷守恒方程如下：

(7)

假定鋰離子電池以恒定電流放電，邊界條件定義：

在負(fù)極極耳上

(8-1)

在正極極耳上

(8-2)

在集流體/電極活性材料界面

(8-3)

1.2.2質(zhì)量守恒方程

(1)固相鋰離子守恒方程

假定固相活性粒子為球狀，固相鋰離子嵌入和脫出電極活性粒子的擴(kuò)散過程服從菲克(Fick)第二定律，以球坐標(biāo)形式表達(dá)如下:

(9)

固相擴(kuò)散過程的邊界條件為

(10)

式中:c1為固相鋰離子濃度;D1為固相擴(kuò)散系數(shù);r為活性顆粒徑向坐標(biāo);Rs為顆粒半徑.

(2)液相鋰離子守恒方程

假定鋰離子在電解液中的傳輸過程遵循濃溶液理論.考慮鋰離子的擴(kuò)散、電遷移過程以及界面反應(yīng)產(chǎn)生/消耗，則液相鋰離子守恒方程如下：

(11)

液相傳質(zhì)過程的邊界條件為：正極與正極集流體界面、負(fù)極與負(fù)極集流體界面Ω上

(12)

1.2.3電極反應(yīng)動(dòng)力學(xué)方程

活性電極粒子表面電化學(xué)反應(yīng)過程由Butler-Volmer方程描述如下:

(13)

j0為交換電流密度

(14)

η為超電勢:

(15)

式中:αa,αc為陽極、陰極的傳遞系數(shù);K0為反應(yīng)速率常數(shù)，僅與溫度有關(guān);Ue為平衡電勢;RSEI為SEI膜電阻;c1,max為最大固相濃度;c1,surf為電極活性顆粒表面鋰濃度.

2控制方程相似準(zhǔn)則的推導(dǎo)

2.1拋物線型偏微分方程

三維傳熱、傳質(zhì)、電流傳導(dǎo)過程方程均具有拋物線型偏微分方程的特征，如表1所示. 因此，傳熱、傳質(zhì)、電流傳導(dǎo)過程的相似準(zhǔn)則亦具有相似性，并均可由拋物線型方程的相似準(zhǔn)則推廣得到. 將拋物線型偏微分方程進(jìn)行量綱為一處理，推導(dǎo)其相似準(zhǔn)則.

表1　不同過程的控制方程相似性

三維拋物線型偏微分方程形式如下[22]：

(16)

通常結(jié)合初始條件及邊界條件以確定偏微分方程解的唯一性.初始條件如下:

(17)

三類邊界條件與邊界定義為：

Dirichlet邊界，(第一類)給定邊界上待求變量的分布

(18)

Newman邊界，(第二類)給出了在邊界處待求變量的導(dǎo)數(shù)或偏導(dǎo)數(shù)

(19)

Robin邊界條件，(第三類)定義待求變量與梯度值之間的函數(shù)關(guān)系

(20)

三類邊界條件可統(tǒng)一寫成如下:

(21)

式中，Γ是邊界.

2.2拋物線型控制方程的相似準(zhǔn)則

對式(16)～式(21)進(jìn)行量綱為一處理，得到量綱為一控制方程如下：

(22)

量綱為一初始條件

(23)

量綱為一邊界

(24)

其中，各綱量為一的表達(dá)式及含義如下：

量綱為一位置坐標(biāo)

(25)

量綱為一源項(xiàng)

(26)

量綱為一吸收項(xiàng)系數(shù)

(27)

量綱為一時(shí)間，表征非穩(wěn)態(tài)熱傳導(dǎo)進(jìn)行的程度

(28)

量綱為一邊界約束

(29)

量綱為一邊界換熱系數(shù)，表征非穩(wěn)態(tài)導(dǎo)熱條件下，物體內(nèi)溫度場的分布規(guī)律

(30)

量綱為一待求解量

(31)

由式(22)～式(24)可知量綱為一場量U為量綱為一的量{X,Y,Z,A,F,τ,B,W}的函數(shù)，亦即為拋物線型方程的相似準(zhǔn)則數(shù)，如式(25)～式(31)所示.相似準(zhǔn)則數(shù){X，Y，Z，A，F(xiàn)，τ,B,W}既可作為模型相似的判定準(zhǔn)則，又可基于相似準(zhǔn)則數(shù)構(gòu)造相似模型.

為解決電池物理場有限元建模的跨尺度難題，建立Z方向(厚度方向)尺寸擴(kuò)大N倍的相似模型進(jìn)行等效模擬. 由相似準(zhǔn)則數(shù){X,Y,Z,A,F,τ,B,W}相等的原則，可導(dǎo)出原型和相似模型的參數(shù)相似設(shè)置關(guān)系，如表2所示.

2.3拋物線型控制方程相似準(zhǔn)則的推廣

由于電化學(xué)過程的物理方程具有拋物線型方程的形式，可將拋物線型方程的相似準(zhǔn)則推廣映射到鋰離子電池電化學(xué)過程控制方程中，得到電化學(xué)過程的原型與模型參數(shù)設(shè)置.

由式(9)可知，球形活性粒子域的固相擴(kuò)散過程是一維過程，活性粒子半徑不變的情形下，固相濃度分布及演變規(guī)律僅與粒子表面邊界通量、固相擴(kuò)散系數(shù)相關(guān). 因此，只需原型與模型固相區(qū)域的表面局部電流密度、固相擴(kuò)散系數(shù)滿足以下關(guān)系，其中，以下標(biāo)p表示原型，以下標(biāo)m表示模型.

(32)

(33)

固相濃度即可滿足分布

表2　拋物線型方程參數(shù)相似關(guān)系設(shè)置(原型：相似模型)

(34)

由于固相電荷守恒(6)、液相電荷守恒(7)、液相傳質(zhì)過程(11)的控制方程均為拋物線型偏微分方程，故各過程的原型與相似模型間參數(shù)相似關(guān)系可以參照拋物線型方程的參數(shù)相似關(guān)系設(shè)置表2直接導(dǎo)出，具體如表3所示.

當(dāng)滿足表3條件時(shí)，模型與原型的電化學(xué)過程相似，即

量綱為一固相濃度

(34)

量綱為一液相濃度

(35)

量綱為一固相電勢

(36)

量綱為一液相電勢

(37)

即待求解物理量的實(shí)際分布滿足關(guān)系如下：

(38)

3基于有限元的相似準(zhǔn)則驗(yàn)證

為驗(yàn)證表2、表3以及式(35)～式(38)所示相似準(zhǔn)則和參數(shù)相似關(guān)系的正確性，基于Comsol Multiphysics有限元軟件平臺建立電池電化學(xué)過程的原型和相似模型，并對比了基于相似原理的三維電化學(xué)模型與準(zhǔn)二維通用電化學(xué)模型在預(yù)測結(jié)果上的差異.

3.1有限元模型設(shè)置

鋰離子電池電化學(xué)模型基本簡化假設(shè)如1.2節(jié)所述，模型參數(shù)根據(jù)相似準(zhǔn)則和電池工況條件設(shè)定，具體參數(shù)如表4所示.模型工況根據(jù)電池試驗(yàn)確定，選定為常溫25 ℃下恒流1C放電至截止電壓2.5 V.

表3　電化學(xué)過程原型與模型參數(shù)相似關(guān)系設(shè)置(原型：相似模型)

表4　電化學(xué)特性參數(shù)

在原型和相似模型的基本結(jié)構(gòu)選擇方面，由于疊層結(jié)構(gòu)各層厚度尺寸為微米級別，與長度和寬度方向尺寸相差達(dá)到3個(gè)數(shù)量級.為保證有限元模型單元長寬比在合適的范圍內(nèi)，分別將厚度方向放大1 000倍和放大2 000倍得到的模型作為原型和相似模型，如圖3所示.其中，原型的單元總數(shù)為15 408，相似模型的單元總數(shù)為30 816.利用Comsol Multiphisics軟件分別對原型和相似模型，進(jìn)行相似準(zhǔn)則和相似關(guān)系的驗(yàn)證計(jì)算.

a原型b相似模型

圖3電化學(xué)過程原型和相似模型

Fig.3Prototype and model of electro-chemical process

正負(fù)極平衡電勢Ue表達(dá)式分別如式(39)和式(40)所示[17]

Ue,pos=2.567 5+57.69[1-tan(100θ+2.916 4)]+0.442 953arctan(-65.419 28θ+64.897 41)+0.097 237arctan(-160.905 8θ+154.590)

(39)

Ue,neg=0.637 9+0.541 6exp(-305.530 9θ)+

(40)

3.2相似準(zhǔn)則有限元驗(yàn)證結(jié)果分析

選取了電化學(xué)過程液相濃度場、固相濃度場、液相電勢場的空間分布和對應(yīng)測點(diǎn)的時(shí)間歷程作為原型和相似模型相似關(guān)系的評價(jià)指標(biāo)，對應(yīng)測點(diǎn)的坐標(biāo)如表5所示.

表5　測點(diǎn)坐標(biāo)

圖4～圖7所示為原型和模型計(jì)算結(jié)果.圖4a、圖4b分別為原型放電結(jié)束時(shí)刻液相濃度空間表面和內(nèi)部分布圖，圖4c和圖4d分別為模型放電結(jié)束時(shí)刻液相濃度空間分布圖，圖5a和圖5b分別為原型放電結(jié)束時(shí)刻固相濃度空間表面和內(nèi)部分布圖，圖5c和圖5d為模型放電結(jié)束時(shí)刻固相濃度空間表面和內(nèi)部分布圖.圖6a和圖6b分別為原型和模型放電結(jié)束時(shí)刻液相電勢空間分布圖.圖7a和圖7b分別為原型和模型對應(yīng)測點(diǎn)液相濃度和固相濃度的時(shí)間歷程曲線.

a 放電結(jié)束時(shí)刻原型液相濃度表面分布

b 放電結(jié)束時(shí)刻原型液相濃度內(nèi)部分布

c 放電結(jié)束時(shí)刻模型液相濃度表面分布

d 放電結(jié)束時(shí)刻模型液相濃度內(nèi)部分布

Fig.4Comparison of electrolyte concentration between prototype and scale model in 1C discharge period

由仿真結(jié)果可知，在恒流放電過程中，負(fù)極區(qū)域電解液較正極區(qū)域濃度高，因而鋰離子從負(fù)極區(qū)域向正極區(qū)域擴(kuò)散和遷移，如圖4所示.固相中負(fù)極區(qū)域鋰離子脫出活性粒子，正極區(qū)域鋰鋰離子嵌入電極活性粒子，導(dǎo)致正極固相濃度升高，負(fù)極固相濃度降低，如圖7b所示.

a 放電結(jié)束時(shí)刻原型液相固度表面分布

b 結(jié)束時(shí)刻原型固相濃度內(nèi)部分布

c 結(jié)束時(shí)刻模型固相濃度表面分布

d 放電結(jié)束時(shí)刻模型固相濃度內(nèi)部分布

Fig.5Comparison of solid phase concentration between prototype and scale model in 1C discharge

由圖4～圖7計(jì)算結(jié)果可知，在忽略有限元計(jì)算的誤差情況下，相似模型與原型的電化學(xué)過程完全相似.即電化學(xué)過程液相濃度場、固相濃度場、液相電勢場的空間分布(最大值、最小值和分布規(guī)律)和時(shí)間歷程的完全一致，并與理論推導(dǎo)結(jié)果式(35)～式(38)相符.

a 放電結(jié)束時(shí)刻原型液相電勢表面分布

b 放電結(jié)束時(shí)刻原型液相濃度表面分布

a 原型和模型對應(yīng)測點(diǎn)液相濃度時(shí)間歷程對比

b 原型和模型對應(yīng)測點(diǎn)固相濃度時(shí)間歷程對比

Fig.7Comparison of time history at corresponding points between prototype and scale model

因此，理論推導(dǎo)確定的相似準(zhǔn)則和相似關(guān)系是成立的，進(jìn)一步地，模型驗(yàn)證結(jié)果表明電化學(xué)有限元建模方法能實(shí)現(xiàn)鋰離子電池三維電化學(xué)過程跨尺度建模的一致模擬.

3.3計(jì)算效率的討論

分別采用原始模型、厚度放大100倍模型與200倍的相似模型進(jìn)行仿真，對比其計(jì)算精度與效率，計(jì)算硬件平臺為24核CPU 64G內(nèi)存的工作站，計(jì)算效率如表6所示.

表6　不同模型計(jì)算效率對比

對比結(jié)果顯示，厚度放大100倍模型與200倍的相似模型相對于原始結(jié)構(gòu)在計(jì)算效率上有顯著優(yōu)勢.

3.4通用準(zhǔn)二維(Pseudo-2D)模型與基于相似原理三維模型計(jì)算結(jié)果對比

目前，通用準(zhǔn)二維(Pseudo-2D)模型[10-12]忽略內(nèi)部濃度場、電勢場展向(各層面內(nèi))狀態(tài)的不均勻性，將導(dǎo)致準(zhǔn)二維模型的仿真誤差.為評價(jià)通用準(zhǔn)二維(Pseudo-2D)模型與基于相似原理的三維模型仿真結(jié)果的差異，分別提取兩種模型不同倍率(1/3C、1C、2C、4C)放電過程工作電壓，如圖8所示.

圖8　準(zhǔn)二維模型與三維模型不同倍率放電階段電壓對比

由圖8可知，在不高于1C放電倍率的情況下，準(zhǔn)二維模型與三維模型電池仿真結(jié)果差異很小. 隨著放電倍率的升高，準(zhǔn)二維模型與三維模型電池工作電壓仿真結(jié)果差異加大.4C放電情況放電容量相對三維模型偏高2.6%，且準(zhǔn)二維模型對電池放電容量的預(yù)測產(chǎn)生正偏差.

圖9　三維模型4C放電結(jié)束時(shí)刻液相濃度分布

圖10　4C放電結(jié)束時(shí)刻固相濃度分布

圖11　4C放電結(jié)束時(shí)刻生熱速率分布

準(zhǔn)二維模型與三維模型工作電壓的差異來源于電池內(nèi)部濃度場、電勢場等展向(面內(nèi))分布的不均勻性，如圖9和圖10所示,進(jìn)而導(dǎo)致內(nèi)部生熱速率分布不均勻,如圖11所示.

因此，在電池系統(tǒng)集成和管理過程中，有必要對基于準(zhǔn)二維模型的狀態(tài)估計(jì)算法進(jìn)行修正；在電池單體結(jié)構(gòu)設(shè)計(jì)則應(yīng)合理設(shè)計(jì)極耳結(jié)構(gòu)，減小內(nèi)部狀態(tài)的不均勻性，以避免電池運(yùn)行過程出現(xiàn)局部過充過放等問題導(dǎo)致電池失效及安全問題的發(fā)生.

4結(jié)論

基于相似理論的鋰離子電池跨尺度建模研究主要得到以下結(jié)論：

(1)從拋物線型偏微分方程基本形式出發(fā)，利用量綱為一分析法推導(dǎo)了拋物線型偏微分方程的相似準(zhǔn)則，并推廣得到電化學(xué)過程控制方程的相似準(zhǔn)則，能夠?qū)崿F(xiàn)電化學(xué)過程時(shí)間歷程和空間分布的一致模擬；

(2)在某一方向幾何尺寸縮放N倍，為實(shí)現(xiàn)電化學(xué)過程完全相似，相應(yīng)方向的液相擴(kuò)散系數(shù)、固相電導(dǎo)率、液相電導(dǎo)率相應(yīng)縮小或放大N2倍，而其他參數(shù)保持不變.

(3)利用相似原理有助于解決鋰離子電池疊層結(jié)構(gòu)跨尺度建模問題，并為鋰離子電池多物理場跨尺度建模問題的解決提供思路.

(4)由于高倍率放電工況下電池內(nèi)部不均勻性加劇，導(dǎo)致傳統(tǒng)的通用準(zhǔn)二維模型容量預(yù)測結(jié)果具有正誤差.因此，在電池單體設(shè)計(jì)及系統(tǒng)集成與管理過程中，有必要對基于準(zhǔn)二維模型的狀態(tài)估計(jì)算法進(jìn)行修正，以避免電池局部過充過放而導(dǎo)致安全隱患.

由于電化學(xué)過程與溫度場存在相互耦合，在后續(xù)研究中將考慮這一耦合作用，建立電化學(xué)-熱耦合模型進(jìn)行電池內(nèi)部溫度場、電化學(xué)狀態(tài)的精確預(yù)測.

參考文獻(xiàn):

[1]Salvadori A, Bosco E, Grazioli D. A computational homogenization approach for Li-ion battery cells: part 1—formulation[J]. Journal of the Mechanics and Physics of Solids, 2014, 65: 114.

[2]黃可龍，王兆翔，劉素琴. 鋰離子電池原理與關(guān)鍵技術(shù)[M]. 北京: 化學(xué)工業(yè)出版社, 2011.

HUANG Kelong, WANG Zhaoxiang, LIU Suqin. The principle and key technology of lithium ion battery[M]. Beijing: Chemical Industry Press, 2011.

[3]李相哲，蘇芳，林道勇. 電動(dòng)汽車動(dòng)力電源系統(tǒng)[M]. 北京: 化學(xué)工業(yè)出版社, 2011.

LI Xiangzhe, SU Fang, LIN Daoyong. Power supply system of electric vehicle [M]. Beijing: Chemical Industry Press, 2011.

[4]郭偉春，付艷恕，王丹，等. 車載鋰離子電池放電性能影響因素研究[J]. 電源技術(shù), 2013, 37(3): 372.

[5]Franco A A. Multiscale modelling and numerical simulation of rechargeable lithium ion batteries: concepts, methods and challenges[J]. RSC Advances, 2013, 3(32): 13027.

[6]Wu W, Xiao X, Huang X. The effect of battery design parameters on heat generation and utilization in a Li-ion cell[J].ElectrochimicaActa, 2012, 83: 227.

[7]Hu X, Li S, Peng H. A comparative study of equivalent circuit models for Li-ion batteries[J]. Journal of Power Sources, 2012, 198: 359.

[8]Han X, Ouyang M, Lu L,etal. Simplification of physics-based electrochemical model for lithium ion battery on electric vehicle. Part I: Diffusion simplification and single particle model[J]. Journal of Power Sources, 2015,278:802.

[9]Guo M, Sikha G, White R E. Single-particle model for a lithium-ion cell: thermal behavior[J]. Journal of The Electrochemical Society, 2011, 158(2): 122.

[10]Fuller T F. Simulation and optimization of the dual lithium ion insertion cell[J]. Journal of The Electrochemical Society, 1994, 141(1): 1.

[11]Doyle M, Newman J, Gozdz A S,etal. Comparison of modeling predictions with experimental data from plastic lithium ion cells[J]. Journal of the Electrochemical Society, 1996, 143(6): 1890.

[12]Yan B, Lim C, Yin L,etal. Simulation of heat generation in a reconstructed LiCoO2 cathode during galvanostatic discharge[J].ElectrochimicaActa, 2013, 100: 171.

[13]Kumaresan K, Sikha G, White R E. Thermal model for a Li-ion cell[J]. Journal of the Electrochemical Society, 2008, 155(2): 164.

[14]Gerver R E, Meyers J P. Three-dimension modeling of electrochemical performance and heat generation of lithium-ion batteries in tabbed planar configurations[J]. Journal of The Electrochemical Society, 2011, 158(7): 835.

[15]Guo M, White R E. A distributed thermal model for a Li-ion electrode plate pair[J]. Journal of Power Sources, 2013, 221: 334.

[16]Smith K, Wang C. Power and thermal characterization of a lithium-ion battery pack for hybrid-electric vehicles[J]. Journal of Power Sources, 2006, 160(1): 662.

[17]Saw L H, Ye Y, Tay A A O. Electrochemical-thermal analysis of 18650 lithium ion phosphate cell[J]. Energy Conversion and Management. 2013, 75: 162.

[18]王豐. 相似理論及其在傳熱學(xué)中的應(yīng)用[M]. 北京: 高等教育出版社, 1990.

Wang Feng. Similarity theory and its application to heat transfer theory[M]. Beijing: Higher Education Press, 1990.

[19]程洪正，張立軍，阮丞，等. 鋰離子電池溫度場相似準(zhǔn)則推導(dǎo)與有限元驗(yàn)證[J]. 同濟(jì)大學(xué)學(xué)報(bào):自然科學(xué)版, 2013，41(8): 1249.

Cheng H, Zhang L, Ruan C. Validation of similarity criterions of temperature field inside lithium-ion battery[J]. Journal of Tongji University: Natural Science, 2013, 41(8): 1249.

[20]張立軍，程洪正，阮丞，等. 鋰離子電池電場相似準(zhǔn)則的推導(dǎo)與有限元驗(yàn)證[J]. 汽車工程, 2013，35(8): 691.

ZHANG Lijun, CHENG Hongzheng, RUAN Cheng,etal. Derivation and FEM validation of similarity criterions of electric field inside of lithium-ion battery[J]. Automotive Engineering, 2013, 35(8): 691.

[21]程洪正，張立軍，阮丞，等. 鋰離子電池應(yīng)力應(yīng)變場的相似準(zhǔn)則推導(dǎo)與有限元驗(yàn)證[J]. 中國科學(xué):技術(shù)科學(xué), 2013,43(6): 670.

CHENG Hongzheng, ZHANG Lijun, RUAN Cheng,etal. Derivation and FEM validation of similarity criteria of thermal stress-strain field inside of lithium-ion battery[J].ScientiaSinicaTechnologica, 2013,43(6): 670.

[22]陳恕行. 現(xiàn)代偏微分方程導(dǎo)論[M]. 北京: 科學(xué)出版社,2005.

CHEN Shuxing. Introduction to modern partial differential equation[M]. Beijiing: Science Press, 2005.

Three-dimension Multi-scale Electro-chemical Model for Lithium-ion Battery Based on Similarity Theory

ZHANG Lijun1,2, CHENG Hongzheng1,2

(1. College of Automotive Studies, Tongji University, Shanghai 201804, China; 2. Synergy Innovation Center of Intelligent New-energy Vehicle， Tongji University, Shanghai 201804, China)

Abstract：A dimensional analysis method was performed to derive similarity criteria and the similarity coefficients of parabolic partial differential equation, and then it was further generalized to deduce the similarity coefficients of the electrochemical process. To validate the similarity criteria and the similarity coefficients, three-dimension finite element models of electrochemical process of lithium ion battery under specified conditions were established with the software of Comsol Multiphysics. The simulation results show that the similarity criteria and the similarity coefficients are certified to be correct and the modelling method based on the similarity theory can achieve equivalent simulation. Moreover, the simulation results reveal that the inner state of the battery tends to be uneven under high rate operation, thus corresponding measurements in design and management based on conventional Pseudo-2D model should be modified to avoid safety issues due to a local over-charge or over-discharge.

Key words：lithium ion battery; electro-chemical field; similarity theory

文獻(xiàn)標(biāo)志碼：A

中圖分類號：U463.51

通訊作者：程洪正(1989—)，男，博士生，主要研究方向?yàn)樾履茉雌噭?dòng)力總成與關(guān)鍵部件開發(fā).E-mail：chenghz2007@163.com

基金項(xiàng)目：國家“九七三”重點(diǎn)基礎(chǔ)研究發(fā)展計(jì)劃(2011CB711201)

收稿日期：2015—06—02

第一作者： 張立軍(1972—)，男，教授，工學(xué)博士，博士生導(dǎo)師，主要研究方向?yàn)樾履茉雌噭?dòng)力總成集成與控制、汽車振動(dòng)與噪聲的分析與控制.E-mail：tjedu_zhanglijun@#edu.cn