把P2D電化學(xué)模型和熱模型放在同一個一維框架里再用實驗數(shù)據(jù)去做電壓與溫度曲線的聯(lián)合擬合校準這件事我前后做過好幾個電芯項目踩過不少坑但也是我認為在COMSOL里做電池仿真最值得投資的一條技術(shù)路線。這篇文章我把從模型搭建、熱源拆解、分步擬合到求解器配置的完整流程寫出來參數(shù)表和排查方法都是可以直接拿去用的適合正在做鋰電池仿真、BMS策略驗證、熱管理設(shè)計的朋友參考。1. P2D模型一維化到底是怎么回事為什么要這樣搭1.1 P2D模型的物理畫面與一維坐標體系很多人第一次聽到一維P2D會覺得矛盾P2D不是偽二維嗎怎么又一維了這里需要先厘清模型的空間維度劃分。P2D模型的完整畫面包含兩個尺度一個是宏觀的電極厚度方向也就是從負極集流體到正極集流體的這個方向通常記為x軸另一個是微觀的活性材料顆粒內(nèi)部鋰離子從顆粒表面向中心擴散的徑向方向通常記為r方向。這兩個方向疊加在一起才是P2D里2D的來源也就是偽二維——它并不是真的在空間上畫了一個二維幾何而是宏觀的x方向加上微觀的顆粒徑向r方向。在COMSOL里實現(xiàn)P2D時幾何建模其實只需要一條從負極集流體到正極集流體的線段把這條線段分成六個區(qū)域負集流體、負極多孔電極、隔膜、正極多孔電極、正集流體再加上最外側(cè)的邊界層。顆粒內(nèi)部的徑向擴散則通過多孔電極節(jié)點的內(nèi)置功能來描述不需要在幾何里把每個顆粒畫出來。這就是一維電化學(xué)模型的準確含義宏觀坐標只有一個維度微觀尺度以內(nèi)部離散方式嵌入。這個架構(gòu)的妙處在于你不需要真的把顆粒球體建成三維CAD也不需要在宏觀方向上劃分密密麻麻的網(wǎng)格就能得到足夠可信的固相濃度分布、液相濃度分布、固相電位、液相電位和Butler-Volmer電化學(xué)反應(yīng)速率。對于電壓仿真和熱源計算來說這些量已經(jīng)是核心輸入了。1.2 一維模型在處理熱源與電壓擬合時的天然優(yōu)勢實際做仿真的人都知道三維電化學(xué)-熱耦合模型雖然看起來很完整、渲染圖很漂亮但真正拿來擬合實驗數(shù)據(jù)時往往非常痛苦。三維模型動輒幾十萬網(wǎng)格一次倍率放電仿真要跑幾小時甚至幾天優(yōu)化算法每迭代一步都要重跑一次根本做不起參數(shù)校準。一維P2D模型則完全不同。它的自由度通常只有幾千到幾萬一次完整的1C恒流放電仿真在普通工作站上幾分鐘就能跑完這給參數(shù)擬合留出了巨大的空間。我常用的做法是先用參數(shù)化掃描快速摸清每個參數(shù)對電壓曲線和溫度曲線的影響趨勢再用優(yōu)化模塊做幾十步迭代擬合這個工作量在三維模型里是不可想象的。另一個容易被忽略的優(yōu)勢是熱源計算的一致性。一維P2D模型輸出的產(chǎn)熱項包括極化熱、可逆熱和歐姆熱這些熱源項在空間上都是沿x方向分布的。如果你后續(xù)想升級成三維熱模型完全可以把一維計算得到的體積產(chǎn)熱率映射到三維幾何的對應(yīng)區(qū)域。但在做擬合校準的階段一維模型加上集總熱容近似已經(jīng)能夠把電壓曲線和表面溫度曲線擬合到相當(dāng)好的精度。1.3 COMSOL中搭建一維P2D模型的基礎(chǔ)配置在COMSOL 6.4里搭建一維P2D模型最直接的方式是使用鋰離子電池模塊的一維電池設(shè)計接口。幾何建立一個長度為幾百微米的一維線段然后按照負極集流體、負極、隔膜、正極、正集流體的順序設(shè)置各個區(qū)域的材料參數(shù)和電極參數(shù)。關(guān)鍵設(shè)置項包括多孔電極的固相體積分數(shù)、顆粒半徑、最大固相鋰濃度、初始荷電態(tài)、交換電流密度、固相擴散系數(shù)、液相擴散系數(shù)、電導(dǎo)率、孔隙率、彎曲因子等。這些參數(shù)并不是隨便填的它們共同決定了模擬電壓和實驗電壓之間的偏差形態(tài)。比如顆粒半徑增大等效擴散路徑變長放電末期電壓跌落會更早出現(xiàn)交換電流密度減小初始極化過電位增大電壓曲線整體會下移。搭建完成后先用一組簡單的參數(shù)跑通模型輸出放電電壓曲線確認腳本流程沒有問題再進入?yún)?shù)校準階段。我第一次做的時候沒有先跑通基準模型就直接懟實驗數(shù)據(jù)結(jié)果擬合失敗半天才發(fā)現(xiàn)是幾何區(qū)域設(shè)置的單位錯了一位浪費了不少時間。先跑通再擬合這是必須養(yǎng)成的習(xí)慣。2. 熱源計算與熱模型耦合的關(guān)鍵拆解2.1 電化學(xué)熱源的四項物理拆解鋰離子電池在充放電過程中的產(chǎn)熱來源可以分為四類極化熱、可逆熱、歐姆熱和混合熱。對于一維P2D模型來說混合熱通常數(shù)值較小工程上更多關(guān)注前三項。極化熱來源于電化學(xué)反應(yīng)過電位偏離平衡態(tài)也就是Butler-Volmer動力學(xué)中過電位與法拉第電流密度的乘積。這部分熱是不可逆的方向永遠是發(fā)熱。過電位越大極化熱越大所以高倍率放電時極化熱占總產(chǎn)熱的比例會顯著上升。可逆熱來源于鋰離子的嵌入和脫出引起的熵變化數(shù)學(xué)表達式是電流密度乘以溫度再乘以熵?zé)嵯禂?shù)dU/dT。這部分熱有一個特點放電時可能是吸熱也可能是放熱取決于活性材料的熵變符號。對磷酸鐵鋰來說放電過程的熵?zé)嵩谀硞€SOC區(qū)間甚至是負值會在局部出現(xiàn)微弱的吸熱效應(yīng)實際擬合溫度曲線時如果忽略這一項數(shù)值上會明顯偏大。歐姆熱包括固相歐姆熱和液相歐姆熱分別來自電子在固相導(dǎo)電網(wǎng)絡(luò)中的輸運損耗和鋰離子在電解液中的輸運損耗。在COMSOL的多孔電極接口中這幾項產(chǎn)熱源是自動計算的但你需要確認自己在熱模型中收集了哪些項避免重復(fù)相加或者漏項。2.2 熱模型的幾何構(gòu)建與邊界條件選擇做一維電化學(xué)模型與熱模型耦合時熱模型的幾何本身仍然是一維的但物理意義是電芯內(nèi)部的集總溫度或者沿厚度方向的一維溫度分布。最簡做法是把整個電芯簡化為一個集總熱容系統(tǒng)只用一個溫度變量描述平均溫度這個溫度的動力學(xué)方程寫作體積熱容乘以溫度變化率等于總產(chǎn)熱減去表面對流散熱。更精細一點的做法是沿厚度方向解一維熱傳導(dǎo)方程考慮極片、隔膜、集流體不同層之間的熱導(dǎo)率和熱容差異。對于常規(guī)的軟包電芯和方形電芯這個做法已經(jīng)能給出不錯的表面溫度預(yù)測而且計算量幾乎可以忽略不計。邊界條件中對流換熱系數(shù)h是最難確定也最容易成為擬合自由參量的量。自然對流情況下h大約在5到15瓦每平方米開爾文強制風(fēng)冷時可以達到20到50。我做擬合時通常先把h作為一個未知參數(shù)留到最后一輪再擬合因為表面溫度對h非常敏感過早擬合會掩蓋電化學(xué)參數(shù)帶來的誤差。2.3 電化學(xué)參數(shù)的溫度依賴耦合回路與收斂問題真正的電化學(xué)-熱耦合并不是簡單地把電化學(xué)產(chǎn)熱當(dāng)作熱模型的源項就夠了因為電解液電導(dǎo)率、固相擴散系數(shù)、反應(yīng)速率常數(shù)都隨溫度呈阿倫尼烏斯關(guān)系變化。也就是說電化學(xué)計算需要溫度熱計算需要電化學(xué)產(chǎn)熱兩者構(gòu)成了閉環(huán)耦合。在實際操作中溫度依賴的引入方式是在材料節(jié)點的參數(shù)定義里寫入Arrhenius表達式。比如把固相擴散系數(shù)寫成D_s乘以exp(-Ea_s除以R乘以(1/T減去1/T_ref))的形式其中Ea_s是固相擴散的活化能T_ref是參考溫度。這樣做之后仿真結(jié)果會更貼近低溫或高溫工況但代價是方程組的非線性增強瞬態(tài)求解時時間步長可能被迫縮小。我個人的做法是在參數(shù)擬合階段先不考慮溫度依賴把等溫P2D模型的產(chǎn)熱輸出到熱模型用實驗數(shù)據(jù)擬合出產(chǎn)熱總量和熱容等基礎(chǔ)參數(shù)穩(wěn)定之后再逐步引入溫度依賴并利用不同環(huán)境溫度下的充放電數(shù)據(jù)去標定活化能。一步到位反而容易出現(xiàn)多個參數(shù)互相補償、誰都擬合不準的情況。3. 電壓與熱源數(shù)據(jù)的擬合校準全流程3.1 擬合前必須準備好的實驗數(shù)據(jù)與預(yù)處理工作想要把P2D模型和實驗數(shù)據(jù)擬合得靠譜光靠一組1C恒流放電的電壓曲線是不夠的。合理的實驗數(shù)據(jù)集至少應(yīng)該包含三部分。第一部分是極低倍率的充放電曲線倍率通常在0.02C到0.05C之間用來近似開路電壓OCV。為什么需要它因為P2D模型中正負極的平衡電位曲線是輸入函數(shù)而你手里通常只有電芯端電壓沒有正負極各自的平衡電位。通過低倍率充放電數(shù)據(jù)和已知正負極材料體系可以反推出正負極平衡電位的大致曲線以及正負極容量配比。這個過程天然就是一個擬合過程但它的優(yōu)先級排在其它參數(shù)之前。第二部分是至少兩個不同倍率下的恒流放電電壓-容量曲線推薦0.5C和1C的組合如果有條件再加一個2C。不同倍率的數(shù)據(jù)能提供不同過電位水平的約束交換電流密度和擴散系數(shù)才能被區(qū)分出來。如果只擬合一個倍率很容易出現(xiàn)參數(shù)多個解都能對上同一條曲線的情況換一個倍率就露餡。第三部分是溫度-時間曲線最好在電池表面布置熱電偶同時記錄環(huán)境溫度。溫度數(shù)據(jù)對熱容和對流系數(shù)的標定至關(guān)重要。嚴格來說要解耦熱容和產(chǎn)熱率最理想的情況是有絕熱熱失控儀或恒溫箱強制對流條件下的兩組數(shù)據(jù)分別用來標定熱容和散熱系數(shù)。數(shù)據(jù)的預(yù)處理同樣不能馬虎。我遇到過實驗數(shù)據(jù)的SOC定義和模型的SOC定義不一致導(dǎo)致的整體偏移也遇到過放電截止后電壓回升段的數(shù)據(jù)被錯誤納入擬合區(qū)間。一般建議把恒流放電段單獨截取電壓低于截止電壓之后的數(shù)據(jù)直接丟棄容量單位統(tǒng)一換算成Ah或者換算成歸一化容量后再送入擬合。3.2 分步擬合從OCV到動力學(xué)參數(shù)再到溫度參數(shù)參數(shù)擬合最忌諱的就是所有參數(shù)一把梭一次性交給優(yōu)化器同時擬合。因為P2D模型的參數(shù)之間存在嚴重的相關(guān)性比如交換電流密度和固相擴散系數(shù)都能影響過電位但它們的特征時間尺度不同完全可以在分步擬合中逐一確定。我的擬合順序是這樣的。第一步先擬合正負極平衡電位曲線和容量配比。這需要用極低倍率充放電數(shù)據(jù)把模型設(shè)置為幾乎無極化狀態(tài)也就是把交換電流密度調(diào)到很大、擴散系數(shù)調(diào)到很大讓端電壓直接反映平衡電位之差。通過微調(diào)正負極初始荷電態(tài)和最大固相濃度把低倍率曲線對齊。第二步固定平衡電位用半電池或?qū)ΨQ電池數(shù)據(jù)擬合正負極的交換電流密度初始值和固相擴散系數(shù)。如果沒有半電池數(shù)據(jù)就退而求其次用全電池低倍率放電的電壓降來約束總極化再通過不同倍率的電壓差來解耦正負極動力學(xué)。第三步是引入熱模型先固定電化學(xué)參數(shù)和產(chǎn)熱計算用溫度曲線的上升斜率標定體積熱容再用降溫段或者穩(wěn)態(tài)段標定對流換熱系數(shù)。這一步的好處是溫度數(shù)據(jù)對熱容和散熱系數(shù)的敏感度極高而對電化學(xué)參數(shù)不那么敏感可以解耦處理。最后一步才是多參數(shù)的精細化更新通常是聯(lián)合電壓和溫度數(shù)據(jù)做一次整體優(yōu)化微調(diào)所有參數(shù)讓兩組曲線的殘差同時降低。3.3 COMSOL中參數(shù)估計的具體操作方法COMSOL從6.x版本開始優(yōu)化模塊里的參數(shù)估計功能已經(jīng)很成熟可以直接把實驗數(shù)據(jù)導(dǎo)入為插值函數(shù)然后在研究設(shè)置里添加參數(shù)估計步驟選擇擬合參數(shù)和目標函數(shù)。擬合參數(shù)的定義需要注意量綱問題。COMSOL的優(yōu)化模塊允許對參數(shù)做對數(shù)變換也就是把你選擇的參數(shù)以log10尺度參與優(yōu)化。這個功能非常實用因為電化學(xué)參數(shù)的量級跨度很大比如擴散系數(shù)可能在1e-15到1e-10之間變化直接在原始數(shù)值尺度上優(yōu)化優(yōu)化器根本無法感知微小變化。取對數(shù)之后參數(shù)變化以數(shù)量級為單位梯度更平滑收斂性會好很多。目標函數(shù)的定義上建議把電壓均方根誤差和溫度均方根誤差分開定義以加權(quán)和的形式合并。我個人經(jīng)驗是電壓權(quán)重先保持為1溫度權(quán)重初始設(shè)為0到0.1等電壓擬合收斂后再逐步提高溫度權(quán)重。如果一開始就溫度電壓等權(quán)重結(jié)果通常兩頭不討好。優(yōu)化算法方面少參數(shù)情況下用Levenberg-Marquardt方法收斂最快它是基于局部敏感梯度的適合候選參數(shù)已經(jīng)離真值不太遠的情況。如果初始值偏差大或者參數(shù)個數(shù)多先用Nelder-Mead或COMSOL里的全局優(yōu)化方法粗搜索一遍再把結(jié)果作為初值切換到Levenberg-Marquardt精修。如果你對自動化流程有更高需求COMSOL 6.4也支持通過Python或MATLAB腳本控制模型修改參數(shù)、運行求解、讀取結(jié)果。我曾經(jīng)用Python寫過一個循環(huán)每次更新一組候選參數(shù)調(diào)用COMSOL求解再把電壓和溫度殘差反饋給外部遺傳算法全程無人值守跑了幾百次迭代效果比內(nèi)置全局優(yōu)化還要靈活。4. 實操參數(shù)參考表與求解器配置4.1 一組可復(fù)用的電化學(xué)與熱參數(shù)基準值這里給出一組我做NCM622/石墨體系軟包電芯時用過的基礎(chǔ)參數(shù)量級上對常見三元體系都有參考價值。注意這些參數(shù)不是萬能值但可以作為你第一輪擬合的起點。電化學(xué)參數(shù)參考范圍參數(shù)負極石墨正極NCM622說明顆粒半徑5-10 μm3-6 μm影響固相擴散路徑長度固相體積分數(shù)0.45-0.550.4-0.5受壓實密度約束最大固相鋰濃度28000-33000 mol/m340000-50000 mol/m3決定容量上限固相擴散系數(shù)1e-15 到 1e-13 量級1e-15 到 1e-14 量級放電末期電壓跌落的主控參數(shù)反應(yīng)速率常數(shù)1e-11 到 1e-10 量級1e-11 到 1e-10 量級單位與交換電流密度設(shè)置相關(guān)孔隙率0.25-0.350.25-0.35影響液相擴散和電導(dǎo)率彎曲因子1.5-31.5-3常用Bruggeman關(guān)系電解液參數(shù)方面典型1M LiPF6體系的參考值包括液相擴散系數(shù)約2e-10到4e-10平方米每秒鋰離子遷移數(shù)約0.3到0.4電解液電導(dǎo)率約0.5到1.5西門子每米并隨鋰離子濃度非線性變化。COMSOL內(nèi)置電解液庫可以自動計算這些值但如果你用的是針對性實驗數(shù)據(jù)也可以手動覆蓋。熱參數(shù)參考范圍體積熱容通常在1.5到3兆焦每立方米開爾文之間鋁集流體和銅集流體的熱容差異不大主要差別在極片和隔膜材料。整體等效熱導(dǎo)率取0.3到1瓦每米開爾文這個值遠低于金屬因為多孔極片的實際導(dǎo)熱路徑受孔隙和接觸電阻限制。對流換熱系數(shù)上面提過自然對流5到15強制風(fēng)冷20到50。4.2 網(wǎng)格、時間步與求解器配置建議一維模型的網(wǎng)格劃分很簡單但很關(guān)鍵。六個域分別設(shè)置最大單元尺寸電極和隔膜區(qū)域用5到10微米量級的網(wǎng)格即可集流體區(qū)域因為只求解電位方程網(wǎng)格可以放粗到幾十微米。顆粒內(nèi)部的徑向離散在COMSOL中是自動完成的一般設(shè)置在10到20個離散點就足夠再多并不會顯著改善精度只會增加自由度。瞬態(tài)求解器的設(shè)置是很多人忽略但影響巨大的環(huán)節(jié)。P2D模型的特征時間尺度跨越很大從電化學(xué)反應(yīng)的毫秒級到熱傳導(dǎo)的百秒級默認的瞬態(tài)求解器如果時間步長控制不當(dāng)要么收斂失敗要么計算時間爆炸。我的習(xí)慣是把求解器設(shè)為BDF最大階數(shù)5相對容差調(diào)到1e-4并設(shè)定最大時間步長為放電總時間的百分之一比如1C放電3600秒最大步長取36秒。需要注意的是在放電初期電壓曲線快速變化階段求解器應(yīng)該能自動加密步長所以不用手動設(shè)置固定的初始步長。不要忽視代數(shù)變量的尺度問題。電壓是伏特量級濃度是幾千摩爾每立方米量級過電位是幾十毫伏量級數(shù)值跨度很大。如果求解器出現(xiàn)收斂困難先檢查各個因變量的尺度是否合理必要時在因變量設(shè)置里手動手工縮放。我遇到過隔膜區(qū)域液相濃度在最開始出現(xiàn)負值的情況就是因為初始濃度和邊界濃度差太大瞬態(tài)求解的振蕩引起的縮小初始步長后問題迎刃而解。模型跑通之后強烈建議先做一個自檢清單放電容量是否和實驗接近、開路電壓是否正確、固相鋰濃度是否始終低于最大濃度、液相濃度是否始終為正、產(chǎn)熱項是否符合物理直覺。這些檢查看起來基礎(chǔ)但能幫你省下大量的擬合調(diào)試時間。5. 常見問題與排查技巧實錄5.1 擬合發(fā)散和參數(shù)跑飛的應(yīng)對策略參數(shù)擬合發(fā)散是每個做電池仿真的人都會遇到的頭疼問題。表現(xiàn)通常是優(yōu)化迭代幾步之后某個參數(shù)跳到了設(shè)定邊界值然后整個仿真收斂失敗優(yōu)化過程中斷。我遇到這種情況的經(jīng)驗是第一優(yōu)先檢查目標函數(shù)的物理尺度。電壓殘差的量級應(yīng)該在毫伏級溫度殘差的量級應(yīng)該在一到兩開爾文左右如果兩者數(shù)值跨度太大優(yōu)化器會本能地優(yōu)先降低數(shù)值更大的一項另一項就失去約束力。加權(quán)系數(shù)在這里要做歸一化處理建議把溫度殘差除以一個參考溫度差電壓殘差除以一個參考電壓差再乘以各自權(quán)重。第二減少同時擬合的參數(shù)數(shù)量。COMSOL參數(shù)估計步驟里有一個靈敏度分析功能可以在正式擬合前跑一遍看每個參數(shù)對目標函數(shù)的影響權(quán)重。那些靈敏度極低的參數(shù)比如某些條件下固相擴散系數(shù)對電壓曲線幾乎沒有影響直接固定為文獻值就好。把它們放進擬合只會增加優(yōu)化器的自由度制造假的相關(guān)性。第三給參數(shù)設(shè)置合理的上下限。這個看似簡單實際上很多人為了不限制模型把邊界設(shè)得很寬結(jié)果優(yōu)化器跑飛。擴散系數(shù)的下限設(shè)置在1e-16量級上限設(shè)置在1e-10量級反應(yīng)速率常數(shù)同理。取對數(shù)尺度后這些邊界大約跨2到3個數(shù)量級足夠模型自由探索但不會飛出物理合理的區(qū)間。5.2 電壓曲線對不上應(yīng)該從哪里開始排查擬合過程中電壓曲線對不上是最常見的問題而且不同的偏差形態(tài)指向的參數(shù)完全不同這個經(jīng)驗非常有用。電壓整體偏高或偏低一個固定值首先檢查初始荷電態(tài)SOC和正負極平衡電位曲線是否對準。低倍率下的端電壓近似等于正極平衡電位減去負極平衡電位并扣除歐姆壓降如果SOC的設(shè)置差了一個百分點電壓就會偏移幾毫伏到幾十毫伏。放電初期電壓迅速下墜過度通常是動力學(xué)參數(shù)偏小導(dǎo)致的極化過大檢查交換電流密度或者反應(yīng)速率常數(shù)的數(shù)量級。這里有一個快速判斷方法看放電開始后幾秒內(nèi)的電壓跳變這個跳變主要由歐姆內(nèi)阻和反應(yīng)的電荷轉(zhuǎn)移過電位貢獻如果跳變量比實驗值大優(yōu)先調(diào)大電導(dǎo)率和反應(yīng)速率常數(shù)。放電末期電壓提前急劇下跌主導(dǎo)參數(shù)是固相擴散系數(shù)或顆粒半徑。電壓曲線的膝部位置直接反映了固相鋰離子濃度是否已經(jīng)接近表面耗盡如果膝部出現(xiàn)得太早說明有效擴散太慢。這里要特別提醒顆粒半徑和擴散系數(shù)的影響方向一致兩者存在強相關(guān)性單憑一條放電曲線很難區(qū)分。解決辦法是換一個不同倍率的數(shù)據(jù)聯(lián)合擬合或者使用不同厚度極片的實驗數(shù)據(jù)。還有一個經(jīng)常被忽略的地方放電C率的定義必須一致。COMSOL里設(shè)置電流或者C率時要注意當(dāng)前是以活性材料容量為基準還是以電池標稱容量為基準C率定義錯誤會讓放電時間整體偏差電壓曲線形狀不變但時間軸對不上這種錯誤最難察覺。5.3 溫度曲線偏低或偏高的典型排查思路溫度擬合偏低的概率比偏高大得多。第一個原因是熱源項漏算很多人在做一維模型時會用電流和端電壓的乘積減去開路電壓和內(nèi)阻發(fā)熱來估計總產(chǎn)熱但P2D模型里電解質(zhì)濃差極化帶來的熱量并不完全包含在這個簡化計算里。最穩(wěn)妥的做法是直接在COMSOL里分別輸出極化熱、可逆熱和歐姆熱然后求和確認總產(chǎn)熱是否與能量平衡吻合。第二個原因是對流換熱系數(shù)設(shè)置過大。我做擬合時曾經(jīng)糾結(jié)于熱容值無論如何調(diào)都不對后來才發(fā)現(xiàn)是沒開強制風(fēng)冷的場景里把h值參考了別的文獻的強制對流數(shù)據(jù)。自然對流加上電池支架的接觸散熱h取8到12是比較合理的第一輪假設(shè)。溫度曲線偏高的情況通常是因為熱容太小或者可逆熱被當(dāng)作純發(fā)熱處理。磷酸鐵鋰體系尤其容易出現(xiàn)這種情況放電過程中的熵?zé)嵩诓糠諷OC區(qū)間吸熱如果不考慮這一項整個溫度曲線會被高估。解決方法是檢查熵?zé)嵯禂?shù)的符號COMSOL支持按SOC區(qū)間給出dU/dT的插值函數(shù)可以直接輸入實驗測得的熵?zé)釘?shù)據(jù)。另外還必須提醒如果熱模型直接把電化學(xué)模塊產(chǎn)生的全部熱作為源項同時又在電化學(xué)模塊里設(shè)置了歐姆熱計算并且熱模型的固體傳熱節(jié)點還自動加了一份電阻熱那么產(chǎn)熱會被重復(fù)計算。這種問題檢查方式很簡單做一個1C放電到50%SOC的瞬態(tài)仿真看總累計產(chǎn)熱是否和總能量損失一致不一致就是重復(fù)計算了。5.4 多目標擬合的權(quán)重調(diào)整與驗證方法同時擬合電壓曲線和溫度曲線本質(zhì)上是一個多目標優(yōu)化問題。用加權(quán)和把兩個目標合在一起時權(quán)重的設(shè)置會直接影響最終參數(shù)解在哪個位置折中。我的建議不是一開始就追求完美的聯(lián)合擬合而是先做純電壓擬合得到一組候選參數(shù)再固定這組參數(shù)單獨看溫度預(yù)測記錄溫度偏差的方向和大小。然后把溫度權(quán)重從0開始逐步增加每輪擬合完都看一下兩個目標的Pareto前沿分布。通常你會發(fā)現(xiàn)電壓殘差稍微增加一點點溫度殘差就能大幅改善這說明當(dāng)前候選參數(shù)處于一個合理區(qū)間如果電壓殘差急劇增加而溫度殘差改善甚微說明權(quán)重已經(jīng)過大了需要回退。最后一步不要省略用一組沒有參與擬合的獨立實驗數(shù)據(jù)做驗證。比如你擬合用的是0.5C和1C放電數(shù)據(jù)那就用0.75C或者2C的數(shù)據(jù)跑一遍仿真看預(yù)測和實驗是否吻合。一個只擬合得好訓(xùn)練數(shù)據(jù)而預(yù)測失效的參數(shù)集意味著過擬合P2D模型的自由度足夠多出現(xiàn)這種情況完全不意外。6. 寫在最后的個人實操體會參數(shù)擬合一圈做下來我最大的體會是擬合的本質(zhì)不是在找一組讓曲線重合的數(shù)字而是在不斷用實驗數(shù)據(jù)修正模型對物理的理解。把參數(shù)拆開、分步擬合、逐步驗證遠比一次性把所有參數(shù)塞給優(yōu)化器靠譜。還有一個經(jīng)驗是前期數(shù)據(jù)質(zhì)量決定了后期擬合效率的上限?；〞r間把OCV測量做準、把倍率數(shù)據(jù)的容量歸一化做好、把溫度測點固定好后面參數(shù)校準省下的時間是以倍計算的。實驗和仿真之間的反復(fù)校準本來就是電池建模的常態(tài)接受這件事之后心態(tài)會穩(wěn)很多。這個一維P2D加熱模型搭好之后能擴展的方向很多。把電化學(xué)產(chǎn)熱映射到三維幾何做熱分布仿真把擬合好的模型參數(shù)用于SOC估算和壽命預(yù)測甚至用它生成虛擬實驗數(shù)據(jù)去訓(xùn)練簡化等效電路模型都是順理成章的下一步。核心是先把這一套電化學(xué)-熱耦合的擬合流程走通后面的路就好走了。