基于电化学-热耦合模型的锂离子电池组件产热分析.pdf

第 12 卷 第 8 期2023 年 8 月Vol.12 No.8Aug.2023储能科学与技术Energy Storage Science and Technology基于电化学-热耦合模型的锂离子电池组件产热分析杨佳兴，张恒运，徐屹东（上海工程技术大学机械与汽车工程学院，上海 201620）摘要：研究电池电化学过程产热对锂离子电池的热管理至关重要。本工作建立了三元NMC锂离子电池的电化学-热耦合模型，首先通过对该电池进行不同倍率的放电与温度实验测试，验证了该模型在电压和温度变化预测准确性。然后针对不同温度下的表现进行模拟仿真研究。在室温下，无论倍率大小，负极产热总是小于正极产热，虽然负极的极化热高于正极，但其可逆吸热较大，导致产热水平低于正极。而随着放电倍率的增加，正极产热所占比例减小，负极所占比例先增加后减小，而集流体产热所占比例持续增加。然而，低温条件下的电池放电表现出与室温情况不同的产热特性，首先，低温导致低倍率负极产热率比例大大增加，负极可逆热为总可逆热的主要贡献热。而高倍率负极产热率减少，正极则呈相反趋势。其次在低温下放电时间随倍率增加呈现不同趋势，高倍率下放电电压快速降低导致放电不完全，在低倍率0.51 C放电运行时出现了电压反弹现象但基本放电完全，这是由于低温限制了负极颗粒内部锂离子及时向外扩散，造成电阻增加与电压快速降低，同时大量产热导致自身温升，从而在低倍率下获得电压反弹并保持持续放电的能力。关键词：锂离子电池；电化学热耦合模型；产热；低温特性；电压反弹doi：10.19799/ki.2095-4239.2023.0082 中图分类号：TM 911 文献标志码：A 文章编号：2095-4239（2023）08-2615-11Heat generation analysis for lithium-ion battery components using electrochemical and thermal coupled modelYANG Jiaxing,ZHANG Hengyun,XU Yidong(School of Mechanical&Automotive Engineering,Shanghai University of Engineering Science,Shanghai 201620,China)Abstract:Understanding the thermal aspects of the battery electrochemical process is essential for managing the heat generated in lithium-ion batteries.Herein,an electrochemical-thermal coupled model was developed for an NMC lithium-ion battery.Experimental tests were conducted to validate the accuracy of the model in predicting voltage and temperature variations at different discharge rates and temperatures.Simulation studies were then conducted to estimate battery performance at different temperatures.At normal temperatures,regardless of the discharge rate,the NE heat generation was always smaller than the PE heat generation.While the NE polarization heat was higher than the PE polarization heat,the reversible heat absorption of NE was larger than that of PE.Consequently,the NE heat generation remained lower than that of PE.As the discharge rate increased,the proportion of PE heat generation decreased,whereas the proportion of NE heat generation increased first 储能测试与评价收稿日期：2023-02-20；修改稿日期：2023-03-14。基金项目：国家自然科学基金（51876113），上海市自然科学基金项目（21ZR1426300）。第一作者：杨佳兴（1996），男，硕士研究生，研究方向为动力电池热管理，E-mail：；通讯作者：张恒运，教授，研究方向为动力电池热管理，E-mail：。引用本文：杨佳兴,张恒运,徐屹东.基于电化学-热耦合模型的锂离子电池组件产热分析J.储能科学与技术,2023,12(8):2615-2625.Citation：YANG Jiaxing,ZHANG Hengyun,XU Yidong.Heat generation analysis for lithium-ion battery components using electrochemical and thermal coupled modelJ.Energy Storage Science and Technology,2023,12(8):2615-2625.2023 年第 12 卷储能科学与技术and then decreased.Simultaneously,the proportion of collector heat generation continued to rise.However,the heat generation characteristics differed when the battery discharge was discharged at subzero temperature compared to the normal temperature case.Firstly,subzero temperatures considerably increased the proportion of NE generation at low discharge rates.The reversible heat from NE became the primary contributor to the total reversible heat.Conversely,at high discharge rates,the NE heat generation rate decreased,while the PE heat generation rate exhibited the opposite trend.Secondly,the discharge time at subzero temperatures showed different trends with increasing discharge rates.At high discharge rates,the voltage rapidly decreased,leading to incomplete discharge and a phenomenon known as voltage rebound.At low discharge rates of 0.51 C discharge operation,the basic discharge process was completed,but the voltage rebound still occurred.The subzero temperature restricted the timely diffusion of lithium ions from the NE particles,leading to increased resistance and a rapid voltage drop.Meanwhile,a large amount of heat generation contributed to the temperature rise in the battery.Therefore,the battery exhibited the ability to rebound and maintain a continuous discharge at low rates.Keywords:lithium-ion battery;electrochemical-thermal coupled model;heat generation;subzero temperature characteristics;voltage rebound锂离子电池由于其高能量密度、快速充电和放电能力、无记忆效应和低自放电率而被认为是最受欢迎的动力源。然而，锂电池的性能与工作条件和环境温度密切相关。温度过低或者过高都对电池的使用寿命和内部组成产生不利影响1。锂离子电池是一个典型的电化学储能装置，包含了各种电化学反应、离子输运以及热传递过程。Bernardi等2和Thomas等3提出了电化学系统及多孔插入电极的通用能量平衡方程，包括不可逆热源项、可逆热源项以及相变和热容量变化项。国内外学者也基于不同产热模型对锂离子电池的电热特性进行了分析4-5，但这些模型在使用过程有所简化，导致其内部各组件产热特性无法评估。Newman等6提出了基于多孔电极理论的准二维模型(P2D模型)，考虑了液态电解液中的浓溶液和固态电极的传输动力学，并利用Fick扩散定律来描述固相材料中锂离子的插入和脱离，自此该模型得到了广泛应用7-11。Lyu等8利用P2D模型研究了高镍NMC储能锂离子电池的热特性。Du等9建立了一个电化学-热耦合模型来考虑锂离子电池内部不可逆热量的产生和演变。张志超等10利用P2D模型研究了层叠式锂离子电池在多尺度的电化学与热行为。需要指出的是，目前的研究大都按照电池各部分热量产生机制进行求解。而在实际应用中，基于电池内部分层组件产热进行区分更加合理，故需要对电池分层组件产热机制进行研究。因此本工作旨在研究NMC锂离子电池分层组件在不同温度、不同倍率下的产热特性。为此，基于P2D模型分析微观层面的电化学浓度和电荷传输并获得电池各层组件的产热变化规律。研究了不同放电倍率下常温与低温产热特性。其中正、负极的极化热，正、负极的可逆热，以及电解液、颗粒的欧姆热，都是通过数值求解确定的。1 电化学-热耦合模型选用额定容量为2.7 Ah、以三元材料为正极材料的N18650CL-29圆柱形锂离子电池进行模型分析和实验测试，其基本参数见表1。电化学-热耦合模型的流程如图1所示，它可以分解为以下几个步骤：拆解真实电池，确定集流体以及电极尺寸，通过建立集流体的等效电阻计算集流体产热量。表1实验中使用的电池的基本参数Table 1Basic parameters of the battery used in the experiment参数标准容量/mAh标称电压/V充电截止电压/V放电截止电压/V数值27003.604.202.502616第 8 期杨佳兴等：基于电化学-热耦合模型的锂离子电池组件产热分析输入电极模型的参数，并根据电极模型计算电极产热。用电池的总产热耦合进三维热模型中计算温升，并将电池温度再度耦合至电极模型中重新计算产热率，实现双向耦合直至迭代收敛。图2是包含负极、隔膜和正极的锂离子电池一维模型，在宏观方向上沿电极厚度方向，在微观方向沿活性材料的固体颗粒的半径方向。1.1电化学模型负极由人造石墨制成，电解液为LiPF6溶液(1 11 EC/EMC/DMC)，正极为LiNi0.5Mn0.3Co0.2O2。P2D模型用表2中的5个控制方程表示，首先对于正负极电极颗粒而言，在电极颗粒内部的物质传输分别用Fick第二定律与欧姆定律表示，Fick第二定律将固相球形颗粒径向上的任意一点锂离子浓度随时间t的变化率转化为方程右侧中该点锂离子浓度随空间的变化率。而用欧姆定律描述了电子在固相颗粒沿电极厚度方向的移动。在电极颗粒内部的电荷守恒用式(3)表示。在电解液中的输运通常用浓溶液理论修正的Nernst-Planck方程来描述，而基于浓溶液理论，给出液相中锂离子的物质守恒方程。等效参数由下式进行描述：eff=iBrug，i=ne，peeff=lBrugDeff=DllBrug(6)图1电化学-热耦合流程Fig.1Electrochemical thermal coupled flow diagram图2电极和颗粒尺度的电化学模型Fig.2Schematic diagram of the electrochemical model in both electrode scale and particle scale表2P2D模型控制方程Table 2P2D model control equations位置正负电极颗粒电解液物质传输电荷守恒物质传输物质守恒控制方程 cs t=1r2 r(Dsr2 cs r)is=-eff sil+is=I，il+is=0il=-eff l+2effRTF(1-t+)(1+dlnfdlncl)lncl cl t=(Deff cl)+(1-t+)ilF(1)(2)(3)(4)(5)边界条件 cs t|r=0=0,-FDs cs r|r=rp=jP,cs|t=0=cs,0is|x=0=is|x=lne+lS+lpe=I,is|x=lne=is|x=lne+lS=0 il=PjP，is=-PjPil|x=0=il|x=lne+lS+lpe=0il|x=lne=il|x=lne+lS=I26172023 年第 12 卷储能科学与技术式中，为固相电导率；为液相电导率；Dl为锂离子液相扩散系数；Brug为Bruggeman系数。电化学反应是锂离子电池产生持续充放电能力的基础，利用Buter-Volmer方程来描述电池内部的电极反应动力学。jP=j0 exp(aFRT)-exp(-cFRT)(7)式中，j0为交换电流密度；a和c为传递系数；为过电位。B-V方程将反应电流密度与过电位联系起来。j0=Fk(clcl，ref)a ccs，surf(cs，max-cs，surf)a(8)式中，k为电极界面反应速率常数；cl,ref为液相锂离子参考浓度；cs,surf为固相颗粒表面锂离子浓度；cs,max为固相颗粒最大锂离子浓度。过电位定义为=s-l-E。s为固相电势，l为液相电势，E为电极固相材料的平衡电势，与电极SOC和温度(T)有关，相应的方程为E=Eref+E T(T-Tref)(9)全电池的熵系数通过测量不同温度下的开路电压得到。采用文献12中负极熵系数，并利用实测的全电池熵系数计算得到正极熵系数13-14，如图3所示。这里电极SOC为锂离子的化学计量浓度，其公式为：x，y=cscs，max(10)1.2产热模型锂离子电池的产热率包括来自负极和正极、集流体与电解液的产热qtot=qne+qpe+qC+qel(11)其中单个正、负电极包含3个主要的热源项，即极化热(qpol)、欧姆热(qOhm)和可逆热(qrev)，由式(12)(14)表示，集流体中的热量由式(15)的欧姆热表示，电解液欧姆热由式(16)计算。产热公式汇总在表3中。1.3模型参数P2D模型参数包括正极和负极颗粒的平均半径、正极和负极的平衡电位、正极和负极的熵系数、电极活性材料的最大锂离子浓度、电极的固相扩散系数、电解液扩散系数、电极的固相电导率、电解液电导率等。本工作中NMC锂电池的参数来自于文献8,15-22，室温下电化学参数见表4，变温动态参数见表5。1.4实验测试验证图4(a)展示了电池热测试实验图。使用K型热电偶和温度采集装置(Agilent 34972A)监测并记录电池温度，时间间隔为3 s。为了更准确地测量电池的放电温度，电池表面去除绝缘膜并焊接热电偶，直接测量电池的金属表面温度。在验证研究图3正负电极熵系数随锂离子化学计量的变化Fig.3NE and PE entropy coefficients with lithium-ion stoichiometry表3产热方程Table 3Heat generation equation位置正负电极颗粒集流体电解液方程qpol=PjPqOhm=eff s s+eff l l+2effRTF(t+-1)(1+dlnfdlncl)lncl lqrev=PjPT E TqC=(RC,Cu+RC,Al)I2qel=eff l l+2effRTF(t+-1)(1+dlnfdlncl)lncl l(12)(13)(14)(15)(16)2618第 8 期杨佳兴等：基于电化学-热耦合模型的锂离子电池组件产热分析中，在室温25，电池外表面的对流传热系数确定为13.4 W/(m2K)23。同时记录电压和电流。根据电池内部集流体结构构建集流体等效电路并得到负极集流体电阻为 9.19 m，正极集流体电阻为3.4 m，该部分产热不能忽略，需要考虑在内。图5为室温(25)下不同放电倍率0.5 C、1 C和3 C下仿真与实验的对比。其中图5(a)为电压对比而图5(b)是温升对比，可以看出该模型在电压和温度上都有很好的一致性。2 结果与讨论2.1室温放电的产热分析本节分析在室温和不同放电倍率下(0.5 C、1 C、2 C、3 C)电池内各组件的产热，不同倍率下的电池平均总产热符合文献24产热与倍率的拟合关系，并着重对0.5 C放电倍率下各组件的瞬时产热进行了分析对比。图6显示了室温下不同放电倍率的负极、电解液、集流体和正极的产热率。由图6(a)可以看出，尽管负极的极化热高于正极，但由于负极的可逆产热存在吸热阶段，因此负极产热量低，而正极产热率大于负极，容易局部升温导致安全问题。而随着放电倍率的增加，正极、负极、电解液与集流体的产热均有所增加。由表 6 可以看出，随着放电倍率的增加，负极产热比例先增加后减小，正极产热的比例在减小而集流体产热的比例在增加，这是由于集流体欧姆产热与电流平方成正表5变温电化学动态参数Table 5Variable temperature electrochemical dynamic parameters参数等效电解液电导率/(S/m)锂离子固相扩散系数Ds/(m2/s)锂离子液相扩散系数Dl/(m2/s)活性系数(f,t+)动态表达式(cref,T)=10-4cref-8.2488+0.053248T-0.00002987T2+0.26235 10-3cref-0.0093063 10-3crefT+0.000008069 10-3crefT2+0.22002 10-6c2ref-0.0001765 10-6c2refT2Ds=Ds,refexp EaRR(1Tref-1T)lgDl(cref,T)=(-4.43-54T-229-0.005cref-0.00022cref)-4(cref,T)=1-0.3631-0.399 0.601-0.24 10-2c0.5ref+0.982 1-0.0052()T-293.15(10-3cref)1.5表4室温下电化学模型参数Table 4Electrochemical model parameters at normal temperature参数厚度/m颗粒半径r/m固相体积分数s电解液体积分数l固相最大锂离子浓度cs,max/(mol/m3)固相初始锂离子浓度cs,0/(mol/m3)电解液参考锂离子浓度cl,0=cref/(mol/m3)传递系数ne,peBruggeman因子Brug锂离子固相扩散活化能EaR/(J/mol)锂离子固相参考扩散系数Ds,ref/(m2/s)等效固相电导率/(S/m)转移数t+电池径向热导率xy/W/(mK)电池轴向热导率z/W/(mK)电池密度/(kg/m3)比定压热容cp/J/(kgK)负极72.570.53860.2807286070.84cs,max,ne10000.51.568025.71.452310131000.3631.6336.9625101028隔膜320.4310001.5正极69.550.63480.3061493880.279cs,max,pe10000.51.59976.8510133.826192023 年第 12 卷储能科学与技术比，导致产热比例快速上升。而从低温下的电池各组件产热比例可以看出，随着放电倍率的增加，正极和负极占比在减小，而集流体产热的比例在增加。2.1.1不可逆热分析图7显示了室温0.5 C下负极、正极的极化热与过电位的比较。可以看出，负极的极化热总是比正极的极化热大。放电结束时，负极的极化热快速增长，而正极的极化热几乎保持不变。这主要是因为在接近放电结束时，大量的锂离子在负极被消耗掉了，导致负极内部严重缺锂的极化状态。负极过电位快速增加作为对比，正极过电位基本不变。欧姆热分为集流体欧姆热、负极和正极颗粒欧姆热以及电解液欧姆热，其中集流体欧姆热是主要贡献热。在电极内部，电解液欧姆热是主要的贡献热。如图8所示，电解液欧姆热大于颗粒欧姆热。随着放电时间的增加，电解液中的热量缓慢减少，固相材料的欧姆热保持不变。正极的欧姆热大于负极的欧姆热，两者受温度的影响较小。2.1.2可逆热分析图9显示了正极与负极在室温与低温时的可逆热比较。在室温下正极的可逆热全程放热，而负极在放电后期才转为放热，而在低温下，0.5 C倍率的负极可逆放热在放电中期超过了正极可逆热，总可逆热的主要贡献来自负极可逆放热。2.2低温放电的产热分析低温限制了锂离子在液体中的移动和在固体中的扩散，使得锂离子在电池内的传输动力不足，导致内阻和产热增加。另外，通常在低温下发现的电压反弹现象25与内部电化学状态的关系也将在本节阐述。图5(a)显示0.5 C、1 C和3 C的放电电压与时间的关系；(b)显示0.5 C、1 C和3 C的放电温度与时间的关系。Fig.5(a)shows 0.5 C,1 C and 3 C discharge voltage versus time;(b)shows 0.5 C,1 C and 3 C discharge temperature versus time图4(a)18650锂离子热测试电池；(b)基于极耳与集流体关系的负极集流体ECM；(c)基于极耳与集流体关系的正极集流体ECMFig.4Schematic of(a)a 18650 Li-ion thermal test battery;(b)NE collector ECM based on lug-collector relationship;(c)PE collector ECM based on lug-collector relationship2620第 8 期杨佳兴等：基于电化学-热耦合模型的锂离子电池组件产热分析表6室温与低温下不同放电倍率各组件产热比例Table 6Heat production ratio of each component with different discharge rates at normal and subzero temperature电池组件负极正极电解液集流体室温0.5 C19.56%59.26%0.95%20.23%1 C29.02%45.96%0.93%24.08%2 C26.3%34.08%0.52%39.11%3 C23.73%28.36%0.58%47.33%低温0.5 C69.34%20.98%0.27%9.41%1 C65.13%20.81%0.92%13.13%2 C36.58%38.73%0.78%23.91%3 C29.96%39.85%0.96%29.23%图6室温下不同放电倍率电池内部组件产热率Fig.6Heat comparison of different heat generation components at 25 图70.5 C下25 时，正负极极化热与过电位随时间变化的曲线Fig.7Plot of polarization heat and overpotential of PE and NE with time at 0.5 C at 25 图8在0.5 C下25 时，集流体欧姆热与不同组件欧姆热随时间变化的曲线Fig.8Collector ohmic heat versus different component ohmic heat at 0.5 C at 25 26212023 年第 12 卷储能科学与技术2.2.1电池内部热源分析图10显示了低温下(-15)不同放电倍率的负极、电解液、集流体和正极产热率。与室温情况相比，低温低倍率(0.51 C)下负极的产热率明显大于正极，而在23 C高倍率，正极产热较负极大。图10(a)为0.5 C放电倍率下的产热，其发生在负极的极化热和可逆热达到了69.34%，而正极的产热量下降到20.98%，其余9.68%的热量发生在电解液和集流体的欧姆热。随着放电倍率的增加，尽管负极极化热一直贡献了最大产热，但在3 C下，集流体产热超过负极极化热。图 11 展示了不同放电倍率下的正负极表面SOC(锂离子的化学计量浓度)的变化。其中小倍率放电时间较长，并且正极表面SOC变化较大，说明此种倍率下放电是完全的。而在2 C与3 C高倍率时，正极表面SOC升高很小，在此倍率下电池放电不完全。2.2.20.5 C放电机理分析图12显示了不同组件在-15 下的产热与电压的关系，以-15 的0.5 C为例。放电3575 s为电压反弹点，同时也对应于最大产热率。该温度下负极的总产热率高于正极，此时负极的极化热远高于正极的极化热，这与室温情况不同。这是因为低温增加了锂离子从负极颗粒向外扩散的阻力，从而增加了产热率。图13显示了0.5 C在-15 下，4个不同时间节点的负极颗粒锂离子浓度。在10 s时负极颗粒内外都图9正负极在0.5 C放电时的可逆热比较Fig.9Comparison of reversible heat for PE and NE图10-15 下不同放电倍率电池内部组件产热Fig.10Comparison of different heat generation components at-15 2622第 8 期杨佳兴等：基于电化学-热耦合模型的锂离子电池组件产热分析有足够的锂离子，但当放电过程进行到3575 s时，负极颗粒表面锂离子浓度降至低点(199 mol/m3)，此刻对应产热曲线的最高点。由于电池产热增加温度升高，加速了锂离子从电极内部向外部的扩散速度，于是负极颗粒表面锂离子浓度先是增加，而在放电后期浓度耗尽(6405 s)，导致放电趋于终止。图12在0.5 C放电下-15 时，每个组件的端电压与产热率的关系Fig.12Terminal voltage versus heat generation for each component at 0.5 C and-15 图11在低温-15 下不同放电倍率正负极表面SOCFig.11Surface SOC of NE and PE at different discharge rates at subzero temperature(-15)图13在0.5 C放电下-15 时锂离子在负极颗粒内的浓度随时间变化的关系Fig.13Lithium-ion concentration within the NE particle as a function of time at 0.5 C at-15 26232023 年第 12 卷储能科学与技术3 结论构建电化学-热耦合模型来研究18650圆柱形锂离子电池在不同工作条件下的产热情况，分析了电池各产热部分在室温和低温下的演变，并以0.5 C为例，分析了低温下放电电压反弹内在机制，归纳出以下结论。（1）分析了室温下电池不同组件的各种产热情况。尽管负极的极化热高于正极，但可逆热吸热导致其产热率低于正极。而随着放电倍率的增加，负极产热在整个电池中所占比例先增大后减小，正极产热占比减小，而集流体产热与电流平方正相关，导致高倍率占比增加。负极极化热在放电末期快速增长的原因是由于缺锂导致负极过电位增加。（2）在-15 的低温产热与室温工况下有明显不同。首先其在低倍率下放电负极的产热量超过了正极，而在高倍率下负极产热低于正极的产热特性。在低倍率下放电时，负极的极化热与可逆热占比较高，而高倍率下放电，正极极化热与可逆热占比较高，而不论何种倍率，集流体产热总是随放电倍率增加。（3）同时在-15 电池0.5 C放电，会出现电压反弹，反弹时间与负极颗粒表面SOC引发的产热极大值和锂离子浓度极小值相对应。电池低锂态运行使得温度升高并加速了锂离子在负极颗粒中的向外运动，导致放电曲线中的电压反弹。符号说明c 锂离子浓度，mol/m3cp 定压比热容，J/(kgK)D 扩散系数，m2/sE 电压，VEaR 扩散系数活化能，J/molF 法拉第常数，C/molhc 对流换热系数，W/(m2K)i 局部电流密度，A/m2j0 交换电流密度，A/m2k 反应速率常数，m2.5/(mol0.5s)L 厚度，mqC 集流体欧姆热，W/m3qel 电解液欧姆热，W/m3qtot 总产热率，W/m3qne,qpe 负极、正极产热率，W/m3R 理想气体常数，J/(molK)rp 颗粒半径，mT 温度，t 时间，st+转移数U 开路电压，V 活性因子希腊符号n,p 负极或正极传递系数 体积分数 液相电导率，S/m 过电位，V 热导率，W/(mK)电池密度，kg/m3 固相电导率，S/m 电势，V角标Brug Bruggeman因子C 集流体eff 有效值el 电解液max 最大值min 最小值ne 负极pe 正极ref 参考值S 隔膜s 固相surf 表面0 初始值缩略词ECM 等效电路模型NE 负极PE 正极参 考 文 献1 FENG X N,OUYANG M G,LIU X,et al.Thermal runaway mechanism of lithium ion battery for electric vehicles:A reviewJ.Energy Storage Materials,2018,10:246-267.2 BERNARDI D,PAWLIKOWSKI E,NEWMAN J.A general energy balance for battery systemsJ.Journal of the Electrochemical Society,1985,132(1):5-12.3 THOMAS K E,NEWMAN J.Thermal modeling of porous insertion electrodesJ.Journal of the Electrochemical Society,2003,150(2):doi:10.1149/1.1531194.4 HU X S,LI S B,PENG HUEI.A comparative study of equivalent circuit models for Li-ion batteriesJ.Journal of Power Sources,2624第 8 期杨佳兴等：基于电化学-热耦合模型的锂离子电池组件产热分析2012,198:359-367.5 BABA N,YOSHIDA H,NAGAOKA M,et al.Numerical simulation of thermal behavior of lithium-ion secondary batteries using the enhanced single particle modelJ.Journal of Power Sources,2014,252:214-228.6 DOYLE M,FULLER T F,NEWMAN J.Modeling of galvanostatic charge and discharge of the lithium/polymer/insertion cellJ.Journal of the Electrochemical Society,1993,140(6):1526-1533.7 李夔宁,谢运成,谢翌,等.基于电化学热耦合模型的富镍锂离子电池产热分析J.储能科学与技术,2021,10(3):1153-1162.LI K N,XIE Y C,XIE Y,et al.Analysis of heat production of nickel-rich lithium-ion battery based on electrochemical thermal coupling modelJ.Energy Storage Science and Technology,2021,10(3):1153-1162.8 LYU P Z,HUO Y T,QU Z G,et al.Investigation on the thermal behavior of Ni-rich NMC lithium-ion battery for energy storageJ.Applied Thermal Engineering,2020,166:doi:10.1016/j.applthermaleng.2019.114749.9 DU S L,LAI Y Q,AI L A,et al.An investigation of irreversible heat generation in lithium ion batteries based on a thermo-electrochemical coupling methodJ.Applied Thermal Engineering,2017,121:501-510.10 张志超,郑莉莉,杜光超,等.基于多尺度锂离子电池电化学及热行为仿真实验研究J.储能科学与技术,2020,9(1):124-130.ZHANG Z C,ZHENG L L,DU G C,et al.Electrochemical and thermal behavior simulation experiments based on multiscale lithium ion batteriesJ.Energy Storage Science and Technology,2020,9(1):124-130.11 韦雪晴,邓海鹏,周宇,等.锂离子电池组的三维电化学-热耦合仿真分析J.储能科学与技术,2022,11(12):3965-3977.WEI X Q,DENG H P,ZHOU Y,et al.Three-dimensional electrochemical-thermal coupling simulation analysis of lithium ion battery packJ.Energy Storage Science and Technology,2022,11(12):3965-3977.12 FANG R Q,GE H,WANG Z H,et al.A two-dimensional heterogeneous model of lithium-ion battery and application on designing