中國儲能網(wǎng)訊：
摘 要 鋰離子電池內(nèi)短路誘因復(fù)雜，為深入研究內(nèi)短路引起的電池失效問題須構(gòu)建合適的精細化仿真模型。本工作以NCM/石墨電池為研究對象，圍繞電池內(nèi)短路失效機理，基于電化學(xué)-熱耦合物理場，建立了考慮熱失控放熱副反應(yīng)的三維單層電芯內(nèi)短路模型，探究了熱失控觸發(fā)邊界，并從內(nèi)外部特征討論了單層電芯內(nèi)短路-熱失控的演變過程。首先利用Arrhenius公式得到內(nèi)短路觸發(fā)的四種放熱副反應(yīng)產(chǎn)熱量與反應(yīng)速率，探究對電池溫升影響最大的副反應(yīng)類別，結(jié)果表明內(nèi)短路過程放熱副反應(yīng)中負極與電解液反應(yīng)總熱量最大。進一步分析單層電芯內(nèi)四種典型內(nèi)短路形式的熱失控觸發(fā)特性，綜合考慮組分材料導(dǎo)電性和導(dǎo)熱性，得到鋁-陽極內(nèi)短路危險程度最高，其短路電阻值與熱失控觸發(fā)時間呈現(xiàn)正相關(guān)趨勢，且臨界短路電阻的高溫?zé)狳c區(qū)域面積值約為30 mm2。模擬結(jié)果獲得了四種形式內(nèi)短路臨界短路電阻值，并揭示了單層電芯內(nèi)短路-熱失控觸發(fā)時內(nèi)部鋰離子濃度和溫度分布的空間演變規(guī)律，相關(guān)結(jié)果可為研究內(nèi)短路失效機制和設(shè)計安全鋰離子電池提供理論指導(dǎo)。
關(guān)鍵詞 鋰離子電池；單層電芯；內(nèi)短路模型；熱失控觸發(fā)
在環(huán)境污染和化石能源危機的雙重壓力下，節(jié)能減排已成為世界各國的共同責(zé)任。電動汽車的發(fā)展和大規(guī)模應(yīng)用成為未來趨勢。在此背景下，汽車行業(yè)不斷尋求安全、高效、低碳的新能源儲能器件。鋰離子電池具有高功率、高能量密度、低自放電和長循環(huán)壽命等優(yōu)點，已成為動力電池的首選。但是，鋰離子電池主要由易燃電解質(zhì)和活性電極材料組成，在極端條件下存在潛在的安全問題。機械濫用、電濫用、熱濫用是電池?zé)崾Э氐闹饕T因，而內(nèi)短路一般被認(rèn)為是熱失控誘因的共性環(huán)節(jié)。因此，解決內(nèi)短路引起的電池失效問題，保證鋰離子電池的安全使用，受到國內(nèi)外學(xué)者廣泛關(guān)注。
鋰離子電池內(nèi)短路是典型的熱電耦合過程，涉及內(nèi)部電化學(xué)反應(yīng)等多尺度多物理場問題。采用計算機斷層掃描成像(CT)等先進測試技術(shù)可以獲取電池內(nèi)部狀態(tài)，但呈現(xiàn)的現(xiàn)象基本上是靜態(tài)和局部的，不能完整、動態(tài)地反映電池的內(nèi)短路演變過程。目前，研究內(nèi)短路機制主要通過電特性和熱特性推斷，但實際內(nèi)短路的熱電響應(yīng)結(jié)果并不能與誘因一一對應(yīng)，而是存在“一對多”“多對一”“多對多”的關(guān)系。所以多因素觸發(fā)的內(nèi)短路機制尚不準(zhǔn)確，有待進一步研究。構(gòu)建考慮電和熱力學(xué)特性的仿真模型是深入研究內(nèi)短路失效機制的有效手段。常見的電池模型例如等效電路模型，用電阻等電子元件來描述電池的電氣特性；而熱力學(xué)模型描述了內(nèi)短路的熱特性，由產(chǎn)熱模型和傳熱耗散模型組成。一些研究工作使用耦合模型來模擬內(nèi)短路的觸發(fā)，等效電路-內(nèi)短路耦合模型耦合簡單，通過連接等效內(nèi)短路電阻來模擬內(nèi)短路過程中的產(chǎn)熱與熱耗散，Ouyang等將此模型應(yīng)用于實際電池管理系統(tǒng)(BMS)中的內(nèi)短路檢測，然而由于該模型沒有考慮內(nèi)短路的熱力學(xué)特性，導(dǎo)致內(nèi)短路演變過程的表征不夠準(zhǔn)確。為了描述內(nèi)短路模型的熱力學(xué)特性，Hu等構(gòu)建了等效電路-熱-內(nèi)短路耦合模型，將一個簡單的一維熱力學(xué)模型添加到等效電路-內(nèi)短路耦合模型中，由于此模型考慮了低維熱力學(xué)模型，對內(nèi)短路電學(xué)和熱物性評估更全面。但上述模型仍然未涉及電池內(nèi)部復(fù)雜的化學(xué)反應(yīng)，以及無法描述電池?zé)崃糠植嫉牟痪鶆蛐浴oyle等構(gòu)建的偽二維(P2D)模型是一種常用的簡化電化學(xué)模型，很好地描述了鋰離子在電池中的擴散和遷移以及活性顆粒表面的電化學(xué)反應(yīng)，因此，電化學(xué)-熱-內(nèi)短路耦合模型被廣泛用于研究內(nèi)短路過程。Fang等利用此3D模型研究了特定形式內(nèi)短路(鋁-陽極、陰-陽極)的溫度響應(yīng)結(jié)果，然而對于其他形式內(nèi)短路沒有進行綜合分析。Zhao等采用多物理三維模型研究了針刺內(nèi)短路的過程，揭示了電池?zé)犴憫?yīng)與電化學(xué)行為的耦合，但在該研究中未考慮熱失控放熱反應(yīng)的影響。進一步地，Liu等通過將多個Arrhenius公式擬合，得到兩個集總表達式構(gòu)建了新簡化模型，并研究了電池電極尺寸與層數(shù)對內(nèi)短路的影響。前人建模工作為鋰離子電池內(nèi)短路研究奠定了基礎(chǔ)，由于內(nèi)短路觸發(fā)熱失控是從局部熱點演變的，因此熱失控放熱副反應(yīng)對觸發(fā)熱失控至關(guān)重要。此外，不同形式內(nèi)短路導(dǎo)致的熱失控結(jié)果不同，探究內(nèi)短路到熱失控的演變過程需要對特定形式進行研究。
針對以上問題，本文從電池內(nèi)短路微觀電化學(xué)-宏觀產(chǎn)熱規(guī)律出發(fā)，利用熱力學(xué)分析，以P2D為基礎(chǔ)建立適用于不同形式內(nèi)短路的單層電芯電化學(xué)-熱耦合模型。在單層內(nèi)短路模型的基礎(chǔ)上，耦合熱失控放熱副反應(yīng)模型，來探究熱失控的安全邊界。通過模型比較了四種形式內(nèi)短路的危險程度，并以危險程度最高的鋁-陽極形式梳理內(nèi)短路到熱失控的觸發(fā)過程，總結(jié)了單層內(nèi)短路-熱失控觸發(fā)時序特征和空間演變規(guī)律。
1 模型構(gòu)建
單層電芯由正極集流體(20 μm)、正極材料(52 μm)、隔膜(25 μm)、負極材料(59 μm)、負極集流體(9 μm)構(gòu)成。在電芯內(nèi)部設(shè)置內(nèi)短路電極(幾何形狀為圓柱體)，其半徑為R0，橫截面積為A，長度為L。內(nèi)短路形成的電子導(dǎo)通發(fā)生在正極與負極部分，將在電極內(nèi)部形成閉環(huán)電路，由于內(nèi)短路電極處溫度變化最明顯，為確保求解的精確性對其附近區(qū)域網(wǎng)格作細化處理，如圖1所示。影響內(nèi)短路熱電響應(yīng)的主要因素是短路電極電導(dǎo)率σ [式(1)]，本模型通過改變電極長度L實現(xiàn)四種形式內(nèi)短路的模擬，并調(diào)整變量電阻rshort對導(dǎo)電特性進行設(shè)置。本工作選取三元鋰離子電池作為模擬對象，正極材料選取Li(NiCoMn)1/3O2，負極材料為石墨，電解質(zhì)材料為LiPF6 in EC∶EMC(3∶7)。模型所使用參數(shù)均來自COMSOL Multiphysics數(shù)據(jù)庫以及文獻[9, 34-38]，表1是內(nèi)短路電極相關(guān)參數(shù)，表2列出了電芯電化學(xué)、熱力學(xué)參數(shù)。
圖1 幾何模型與內(nèi)短路模型
表1 內(nèi)短路電極參數(shù)
表2 電化學(xué)及熱力學(xué)參數(shù)
1.1 耦合關(guān)系和控制方程
圖2是本文單層電芯耦合模型，耦合模型由4個子模型構(gòu)成。內(nèi)短路模型為電化學(xué)模型提供相應(yīng)的邊界條件，電化學(xué)模型模擬電化學(xué)響應(yīng)，例如電壓、電流以及電化學(xué)反應(yīng)。相關(guān)參數(shù)列于表2，變量cs、ce分別表示固相嵌入鋰濃度、液相電解質(zhì)鹽濃度，j是摩爾通量。利用熱模型來表示電池溫度分布，正常工作條件下溫度的變化來源于內(nèi)部的產(chǎn)熱Qtotal，以及與外界的熱量交換，總產(chǎn)熱Qtotal由內(nèi)短路產(chǎn)熱Qshort、熱失控副反應(yīng)熱Qex、電化學(xué)熱Qchem共同構(gòu)成，而電化學(xué)產(chǎn)熱包含可逆熱Qrev、不可逆熱Qirrv與歐姆熱Qohm。相關(guān)控制方程是式(10)~(12)，其中Tsurf是電芯表面溫度，U是電芯開路電壓。
圖2 單層電芯3D耦合模型
隨著內(nèi)短路的演變，當(dāng)電芯溫度上升到一定程度后，會觸發(fā)副反應(yīng)產(chǎn)生巨大熱量，因此進一步耦合放熱副反應(yīng)模型，副反應(yīng)產(chǎn)熱為Qex，僅考慮SEI膜分解反應(yīng)，負極與電解液反應(yīng)，正極與電解液反應(yīng)，以及電解液分解反應(yīng)。放熱副反應(yīng)表達式以及參數(shù)值(表3)來源于文獻[40]，耦合模型控制方程均列于表4，其中?s、?e分別表示固相電勢、液相電勢。
表3 副反應(yīng)表達式參數(shù)值
表4 模型控制方程
1.2 邊界條件與模型驗證
電化學(xué)模型與內(nèi)短路模型之間通過電池短路電流Ishort耦合，Ishort通過對整個陽極或陰極域上的反應(yīng)電流密度進行積分獲得，如式(19)所示：
式中，a.和c.分別為陽極和陰極域；j Li為電極反應(yīng)電流密度。
內(nèi)短路產(chǎn)生的焦耳熱由電化學(xué)-熱多物理場耦合接口通過局部電流密度和短路電阻計算：
電芯其他區(qū)域作正常放電處理，正常放電過程中，陰、陽兩個電極的平衡電位與電解液中鋰離子濃度相關(guān)，因此設(shè)置電解液鋰濃度初始值如下：
由于僅在隔膜與電極內(nèi)部存在電解液，因此對電極與集流體界面處發(fā)生的式(4)與式(6)設(shè)置零通量邊界條件，如式(22)所示：
同時在其他邊界設(shè)置：
電化學(xué)模型是耦合模型的關(guān)鍵，可以正確反映電池內(nèi)部鋰離子狀態(tài)。因此，為了驗證模型的電化學(xué)性能，將無內(nèi)短路狀態(tài)的電池在不同的充放電倍率(0.5C、0.33C、0.2C)下電化學(xué)模擬結(jié)果(虛線)與Feng等實驗數(shù)據(jù)(實線)對比，如圖3所示。兩者吻合度高，驗證了電化學(xué)模型和參數(shù)的可靠性。
圖3 無內(nèi)短路時模型與實驗對比結(jié)果
1.3 數(shù)值模擬求解
使用多物理場軟件COMSOL Multiphysics (6.0)版本，采用有限元數(shù)值求解法，添加“鋰離子電池”和“固體傳熱”接口，最后耦合多物理場電化學(xué)熱來模擬內(nèi)短路過程。為確保計算精度，進行了網(wǎng)格無關(guān)性驗證。以短路電阻r=5 Ω的鋁-陽極形式內(nèi)短路為例作網(wǎng)格無關(guān)性檢驗。對內(nèi)部中心有短路電極的單層電芯生成四邊形結(jié)構(gòu)網(wǎng)格，并對內(nèi)短路電極附近區(qū)域作網(wǎng)格細分化處理。圖4給出了當(dāng)單層電芯發(fā)生內(nèi)短路后其內(nèi)部最高溫度。當(dāng)相對容差設(shè)置為10-3時，計算結(jié)果隨著網(wǎng)格數(shù)的增加能夠趨于穩(wěn)定，當(dāng)網(wǎng)格數(shù)為26108時，最高溫度為560 K。因此，選擇具有該網(wǎng)格數(shù)和相對容差的模型作為后續(xù)的數(shù)值計算。模擬計算時長設(shè)置為50 s，時間步長設(shè)置為0.1 s，初始時刻內(nèi)短路電極電導(dǎo)率設(shè)置為極低值，通過引入平滑階躍函數(shù)使模型在1 ms計算后將電導(dǎo)率從極小升高到原始值，增加模型可收斂性。
圖4 網(wǎng)格無關(guān)性檢驗
2 典型內(nèi)短路形式
2.1 熱失控觸發(fā)特性
不同形式內(nèi)短路風(fēng)險程度不同，而內(nèi)短路劇烈程度決定后續(xù)電池是否演變?yōu)闊崾Э?。?nèi)部短路會促進區(qū)域高溫和系列放熱副反應(yīng)的發(fā)生，從低溫到高溫排序，將依次經(jīng)歷如SEI膜分解反應(yīng)，負極與電解液反應(yīng)，正極與電解液反應(yīng)，電解液分解反應(yīng)四個過程。以內(nèi)短路電阻值小于1 Ω的內(nèi)短路為例，計算其放熱副反應(yīng)的產(chǎn)熱量與反應(yīng)速率。圖5(a)是平均產(chǎn)熱量的變化，在t=0~4.5 s時間段內(nèi)，負極與電解液反應(yīng)產(chǎn)生的總熱量最大。然而，在t=0.5~1 s的短暫時間內(nèi)，正極與電解液反應(yīng)迅速增長，顯著促進了產(chǎn)熱過程。此外，正極與電解液反應(yīng)與電解液分解反應(yīng)幾乎同時發(fā)生，且反應(yīng)速率迅速達到最大值。對于熱失控的預(yù)防，盡管冷卻措施可以改變熱失控發(fā)生的時間和溫度，但并不能從根本上抑制其發(fā)生。因此，預(yù)防熱失控的關(guān)鍵在于在熱失控安全臨界值之前控制電池的溫度。當(dāng)電池溫度達到正極與電解液反應(yīng)的初始溫度時，由于反應(yīng)釋放的大量熱量和電解液分解反應(yīng)的迅速進行，此時已難以有效遏制熱失控的進一步發(fā)展。因此，本文后續(xù)內(nèi)短路是否發(fā)展到熱失控的臨界條件以正極與電解液的反應(yīng)觸發(fā)溫度[443.15 K，式(17)]來判斷。
圖5 副反應(yīng)平均產(chǎn)熱量與副反應(yīng)進行速率
2.2 四種形式內(nèi)短路特征
根據(jù)短路電極實際接觸方式將內(nèi)短路分為四種形式，即陰極-陽極、陰極-銅、鋁-銅、鋁-陽極。內(nèi)短路引發(fā)的溫度變化由產(chǎn)熱和散熱能力綜合決定，在放熱副反應(yīng)和對流換熱條件不變時，產(chǎn)熱主要為內(nèi)短路阻抗差異帶來的歐姆熱變化，而散熱則由材料固有導(dǎo)熱性能主導(dǎo)。圖6說明了不同內(nèi)短路電阻時電芯最高溫度變化結(jié)果，并介紹了四種形式內(nèi)短路特征。對于陰極-陽極形式內(nèi)短路[圖6(a)]，由于正負極材料的導(dǎo)電能力遠低于鋁銅金屬，其內(nèi)短路電阻在四種形式中最高，產(chǎn)熱速率低導(dǎo)致溫升較慢。陰極-銅短路接觸時[圖6(b)]，正極材料高阻抗同樣具有低產(chǎn)熱特性，而且銅的熱導(dǎo)率高達400 W/(m·K)，導(dǎo)致陰極-銅形式內(nèi)短路散熱能力強于陰極-陽極形式，表現(xiàn)出溫升進一步降低趨勢。因此當(dāng)電池發(fā)生局部陰極-陽極形式內(nèi)短路或陰極-銅形式內(nèi)短路，其溫度緩慢升高模式一般不會觸發(fā)熱失控。
圖6 四種形式內(nèi)短路不同短路阻值下的熱響應(yīng)曲線
當(dāng)鋁與銅或負極短路接觸時發(fā)生鋁-銅[圖6(c)]或鋁-陽極形式內(nèi)短路[圖6(d)]，鋁金屬材料良好的導(dǎo)電性決定了兩者低阻抗的特性，產(chǎn)熱能力較強。然而銅金屬同時具有高導(dǎo)熱能力，導(dǎo)致鋁-銅內(nèi)短路散熱能力在四種形式中最高，可以及時將熱量散去。對于鋁-陽極形式，由于負極材料的低導(dǎo)熱率，其散熱能力遠差于前者，導(dǎo)致小短路電阻值(r<8 Ω)時，溫度呈現(xiàn)更高上升趨勢，且可以明顯觀察到r=1 Ω時，溫升速率更大。綜上所述，當(dāng)電池發(fā)生鋁-陽極形式內(nèi)短路時，其溫度驟升模式擁有最大概率觸發(fā)熱失控。因此本工作針對熱失控危險程度最高的鋁-陽極內(nèi)短路，探究內(nèi)短路到熱失控的觸發(fā)時序特征，并總結(jié)內(nèi)短路到熱失控的空間演變規(guī)律。
3 單層內(nèi)短路-熱失控觸發(fā)時序特征
電池發(fā)生內(nèi)短路時，溫升由局部內(nèi)短路處向外蔓延，之后內(nèi)短路的演變過程也取決于放熱副反應(yīng)的觸發(fā)與進行，因此探究副反應(yīng)觸發(fā)過程對研究內(nèi)短路-熱失控觸發(fā)有重要作用。本節(jié)通過監(jiān)測放熱副反應(yīng)徑向歸一化濃度變化來探究副反應(yīng)的觸發(fā)時間，如圖7所示，點A至點B徑向長度上每隔0.2 mm設(shè)置一個歸一化濃度監(jiān)測點，當(dāng)點A的歸一化濃度初始值開始變化時，則證明電池開始放熱反應(yīng)，處于內(nèi)短路到熱失控的演變過程。結(jié)合前文分析以正極歸一化濃度ca>0.04作為熱失控的觸發(fā)條件，將電池分為正常區(qū)域和熱失控觸發(fā)區(qū)域。當(dāng)任一點ca值大于0.04判定此處已達到熱失控觸發(fā)條件，將從正常區(qū)域演變到熱失控觸發(fā)區(qū)域，且熱失控會傳播到下一個正常區(qū)域。
圖7 內(nèi)短路放熱副反應(yīng)徑向歸一化濃度監(jiān)測
3.1 放熱副反應(yīng)的觸發(fā)過程
鋁-陽極是危險程度最高的內(nèi)短路形式，圖8是其短路電阻為5 Ω時的放熱副反應(yīng)觸發(fā)過程。SEI膜歸一化濃度監(jiān)測結(jié)果如圖8(a)所示，可以看到在0.6 s后內(nèi)短路電極處負極區(qū)域點A達到SEI膜分解反應(yīng)的起始溫度，即開始了放熱反應(yīng)的過程，1 s時點B處濃度也發(fā)生變化，證明1 s以后由點A到點B的區(qū)域都開始了放熱副反應(yīng)，且將擴散到下一個區(qū)域。圖8(b)是負極歸一化濃度監(jiān)測結(jié)果，可以觀察到即使放熱副反應(yīng)已進行到B區(qū)域，A區(qū)域負極與電解液反應(yīng)也不會達到反應(yīng)極限(負極歸一化濃度消耗到0)。圖8(c)是正極歸一化濃度曲線圖，在2 s左右，內(nèi)短路電極處正極區(qū)域點A達到正極與電解液分解反應(yīng)的起始溫度(熱失控觸發(fā)溫度)，可判定此時局部短路處進入熱失控觸發(fā)階段。放大圖可知在5~10 s內(nèi)，點B所在區(qū)域也達到熱失控觸發(fā)溫度，證明此時整個A到B區(qū)域都演變到熱失控觸發(fā)區(qū)域。而電解液的分解反應(yīng)如圖8(d)所示，在約2~5 s內(nèi)點A處才達到反應(yīng)觸發(fā)溫度。這說明電解液分解反應(yīng)進一步產(chǎn)生的熱量由內(nèi)短路區(qū)域向電池正常區(qū)域傳遞，促進了正極與電解液的反應(yīng)，即促進了熱失控的傳播。
圖8 鋁-陽極內(nèi)短路r=5 Ω放熱副反應(yīng)觸發(fā)過程
3.2 放熱副反應(yīng)的觸發(fā)時序特征
將活性物質(zhì)濃度變化0.0001判定為開始進行放熱副反應(yīng)，對上述4個副反應(yīng)的觸發(fā)時間進行總結(jié)，如圖9(a)所示。對于SEI膜分解反應(yīng)和負極與電解液反應(yīng)，相鄰監(jiān)測點觸發(fā)時間差異不會太大。比較正極與電解液反應(yīng)和電解液分解反應(yīng)的觸發(fā)時間，在相鄰監(jiān)測點上正極與電解液反應(yīng)的觸發(fā)間隔明顯更短，即電池相鄰區(qū)域間熱失控的觸發(fā)間隔較短，而電池整體呈現(xiàn)越靠近內(nèi)短路區(qū)域處熱失控觸發(fā)時間越短，遠離內(nèi)短路電極處的觸發(fā)時間更長的趨勢。這也說明若電池內(nèi)短路附近區(qū)域達到了熱失控觸發(fā)溫度，熱失控演變趨勢向正常區(qū)域擴散概率極大。對比圖9(a)和(b)不同內(nèi)短路電阻值影響下的副反應(yīng)觸發(fā)時間，當(dāng)內(nèi)短路阻值r=5 Ω時，從內(nèi)短路處至距離1 mm位置，觸發(fā)時間由2 s增加至5 s，而當(dāng)r=1 Ω時，觸發(fā)時間由1 s增至2 s。因此，隨著短路電阻值的縮小，觸發(fā)時間縮短，但總體觸發(fā)趨勢仍符合上述分析。
圖9 鋁-陽極內(nèi)短路放熱副反應(yīng)觸發(fā)時間
4 單層內(nèi)短路-熱失控觸發(fā)空間演變
4.1 Li+輸運與電流分布
單層內(nèi)短路-熱失控演變內(nèi)部特征主要體現(xiàn)在Li+的輸運分布，如圖10所示。在t=0 s時，電池處于平衡狀態(tài)，假設(shè)Li+濃度分布均勻，當(dāng)發(fā)生內(nèi)短路的瞬間電池放電，陰極電解液中Li+嵌入正極，負極Li+脫入陽極電解液中。由于鋁-陽極內(nèi)短路電極與鋁集流體接觸，使電極與鋁直接導(dǎo)通，高導(dǎo)電率的鋁導(dǎo)致此區(qū)域流通更多的電子，與陰極電解液中的Li+快速結(jié)合，導(dǎo)致在短時間(t=0.5 s)內(nèi)，隔膜/陰極電解液中Li+濃度幾乎消耗至0，Li+濃度下降趨勢呈現(xiàn)從陰極/隔膜界面處擴散至整個陰極的現(xiàn)象。
圖10 鋁-陽極內(nèi)短路r=5 Ω徑向電解液Li+濃度變化
一方面離子電導(dǎo)率隨Li+濃度的降低而降低，另一方面交換電流密度隨著電解液中Li+濃度的降低而減小[式(7)]，圖11是內(nèi)短路區(qū)域電流密度分布，受多孔電極材料阻值Li+輸運限制，內(nèi)短路電流主要集中在短路位置周圍的較小區(qū)域。由于內(nèi)短路電極與鋁直接導(dǎo)通，因此鋁集流體靠近內(nèi)短路電極邊緣處電流密度更大。電流密度越大，在未熱失控的情況下電池產(chǎn)生的歐姆熱就越大，這也解釋了局部短路處的顯著加熱。
圖11 鋁-陽極內(nèi)短路r=5 Ω，t=0.5 s內(nèi)短路區(qū)域電流密度分布
4.2 溫度邊界演化過程
單層內(nèi)短路-熱失控演變外部特征主要體現(xiàn)在溫度變化，其三維溫度分布如圖12所示。局部加熱導(dǎo)致電芯內(nèi)短路附近溫度最高，由于電池三個方向熱導(dǎo)率不同，電池表面溫度存在圓環(huán)等溫線。圖13是溫度變化趨勢，當(dāng)發(fā)生內(nèi)短路后，最高溫度急劇上升，在達到最大溫度峰值后又平穩(wěn)降低。溫度由中心短路處向外擴散，其他大部分區(qū)域卻保持溫度較低。總的來說，電芯呈現(xiàn)左右區(qū)域溫度對稱分布但整體溫度分布不均的現(xiàn)象。
圖12 鋁-陽極內(nèi)短路r=5 Ω，t=30 s內(nèi)部等溫分布
圖13 鋁-陽極內(nèi)短路r=5 Ω溫度變化趨勢
前述分析將443.15 K設(shè)置為觸發(fā)熱失控的安全邊界，因此本工作將溫度大于443.15 K的電芯區(qū)域判定為高溫?zé)狳c區(qū)域，當(dāng)足夠大的高溫?zé)狳c區(qū)域達到觸發(fā)溫度時，觸發(fā)熱失控的趨勢越明顯。如圖14所示是鋁-陽極形式內(nèi)短路的熱點區(qū)域變化，觀察短路電阻為5 Ω的面積變化趨勢，0~4 s斜率增長不穩(wěn)定，推測此時發(fā)熱的主導(dǎo)因素是內(nèi)短路產(chǎn)生的焦耳熱，而在4~18 s熱點區(qū)域增長斜率幾乎不變，表明此階段放熱副反應(yīng)占主導(dǎo)作用，使電芯溫度發(fā)展到最高階段。18 s以后由于放熱副反應(yīng)的減弱，熱點區(qū)域增長斜率也明顯減小，電芯處于溫升穩(wěn)定階段。觀察8 Ω電阻值(鋁-陽極內(nèi)短路臨界電阻值)熱點區(qū)域面積變化，其發(fā)展分為三個階段，在早期(0~8 s)內(nèi)短路產(chǎn)生的焦耳熱不足以引起大幅度電池的溫升，其后(8~10 s)隨著焦耳熱增大引起的連鎖副反應(yīng)，熱點區(qū)域迅速增長到最大值，最后被限制在30 mm2左右，該值也可以作為內(nèi)短路觸發(fā)熱失控的安全邊界。
圖14 鋁-陽極內(nèi)短路熱點區(qū)域面積變化
5 結(jié) 論
本工作構(gòu)建了3D單層電芯內(nèi)短路模型，系統(tǒng)分析了單層電芯內(nèi)短路到熱失控的觸發(fā)和演變過程，探究了熱失控觸發(fā)邊界值，通過模型分析發(fā)現(xiàn)：
（1）內(nèi)短路會促進區(qū)域高溫和系列放熱副反應(yīng)的發(fā)生，正極與電解液的分解反應(yīng)對電池溫升影響最大，本工作觸發(fā)熱失控的臨界條件以此反應(yīng)觸發(fā)溫度(443.15 K)判定。
（2）以鋁-陽極內(nèi)短路梳理放熱副反應(yīng)的觸發(fā)過程，結(jié)果表明靠近內(nèi)短路區(qū)域的觸發(fā)時間明顯小于電芯其他區(qū)域，同時短路電阻值越小觸發(fā)越快。
（3）鋁-陽極內(nèi)短路由于電極與高導(dǎo)電率鋁直接導(dǎo)通，導(dǎo)致接觸區(qū)域流通更多的電子，影響陰極內(nèi)部Li+輸運演變，最終導(dǎo)致電芯溫度分布不均勻性，其觸發(fā)熱失控的臨界溫度區(qū)域為30 mm2左右。
本工作建模分析彌補了單層內(nèi)短路到熱失控的觸發(fā)和演變過程的不足，提供了一種熱失控觸發(fā)安全邊界的判定方式。后續(xù)可將模型從單層電芯擴展到多層電芯，并泛化以適用于不同內(nèi)短路模式，例如改變電池SOC狀態(tài)，內(nèi)短路電極尺寸或位置等。該熱失控觸發(fā)研究也可針對不同種類的電池，以建立大規(guī)模熱失控觸發(fā)安全邊界數(shù)據(jù)庫，為未來電池?zé)崾Э匕踩A(yù)警提供指導(dǎo)。
